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Estimation of the Relaxation Factor 
for Small Mesh Size 


Recently Frankel [1] and Young [3, 4] have developed a method of successive 
overrelaxation for the solution of the difference equation analogue of the Dirichlet 
problem for a linear elliptic partial differential equation. The purpose of this note 
is to present a simple asymptotic theory of the overrelaxation process which de- 
pends on the assumption that the mesh size h approaches zero, but which de- 
scribes adequately the practical aspects of the earlier work. Thus the results 
discussed are not essentially new, but it is hoped that our approach to them will 
provide an easy insight into the overrelaxation scheme for those interested in its 
application to specific examples. 

We consider the Laplace difference equations for an unknown function u of 
two independent variables in a region D covered by net points spaced / units 
apart. In a standard notation, we use subscripts i and 7 to refer to the location of 
net points in our grid and we use the superscript to indicate steps in the relaxa- 
tion process, so that, in its simplest form, the method of successive overrelaxation 
can be described by the equation 


(1) 4(uzt* r uz ;) = r(upt} ; + up fi, + U5 41,5 + Uj, i441 = 4uz ;), 


where 7 is the relaxation factor. We express r in the form 


2 


(2) 14+ Ch 


for any positive value of the constant C, and we rearrange (1) to obtain 


n n n n n n+l n n+1 n 
ga1,5 + 553-1 + i435 + i541 — 45 iG Mg — Mg H Mi 
h? h? 


n+l n n+1 ” n+1 
+ M55 — Mig — Mj 1 Mi, 5-1 





(3) 


ele 35. 
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Using the familiar idea that the index refers to a new time variable, ¢, and that 
it indicates the location of new net points spaced at time intervals equal to the 
original mesh size h, we recognize that (3) is the difference analogue of the hyper- 
bolic partial differential equation 


(4) Ure + Uyy = Use + Uye + 2C uy. 


Thus for small values of 4 the convergence of the iterative scheme (1) can be 
investigated by an analysis of the decay of time-dependent terms in the solu- 
tion of (4). 
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184. ESTIMATION OF THE RELAXATION FACTOR FOR SMALL MESH SIZE 
The substitution s = ¢ + x/2 + y/2 brings (4) into the canonical form 
(5) Urz + Uyy -T Ftlss _ aC us = 0. 


For a fixed set of boundary conditions, the method of separation of variables 
yields the representation 


(6) «= Uslx, 9) + E [on exp (— fas) + bn exp (—Gus)]Ua(s, 9) 


for the solution u of (5), where Up is the steady-state solution, where the a,, and 
b,, are Fourier coefficients, where 


(7) Pm = 2C — (4C* — 2km?)*, Gm = 2C + (4C* — 2kp*)!, 

and where the U,, and k,,” are the eigenfunctions and eigenvalues of the problem 
(8) AUn + kn?Um = 0 

with homogeneous boundary conditions. The exponent 

(9) b = Re [pi] = Re [2C — (4C? — 2k,*)*] 


corresponding to the lowest eigenvalue k:? governs the rate of decay of the terms 
on the right in (6) with increasing time ¢. Thus the error EZ in the solution of the 
Laplace difference equations obtained through N iterations of the overrelaxation 
scheme (1) has the order of magnitude 


(10) E = O(e-?¥4), 


since it is approximated by the infinite series used in the representation (6) for 
time ¢ = Nh. It follows that we must iterate N = O(1/h) times to arrive at a 
desired degree of accuracy in the solution. 

By (9), the choice of the positive constant C which maximizes p and therefore 
yields the most rapid convergence is clearly C = k,/2', and an underestimate of 
this best value is less damaging than an overestimate (cf. [3]). Indeed, the 
radical in the expression for p is pure imaginary and can be omitted when 
C < k,/2', whereas this same radical causes a sharp decline in the value of p 
when C is taken slightly larger than k,/2*. A simple lower bound on the principal 
frequency k; for the case of Dirichlet boundary conditions is given by Pélya and 
Szegé [2]. If A denotes the area of the region D, they show that 


(11) kiAt > kr, 


where k = 2.405 denotes the first root of the Bessel function of the first kind of 
order zero. Thus a good approximate formula for the relaxation factor r is 


2 2 


(12) Jes 1 + (x/2A)*kh i 1+ 3.014h/Ai 
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Numerical examples worked out by Young [4] indicate that when (12) is used 
instead of the exact best value of r, the additional iterations required to achieve 
equivalent accuracy will usually not exceed 20 per cent of the number needed 
with the optimal r. 

Our approach shows that the method of successive overrelaxation applies to 
the nine-point Laplace difference equation 


(13) tegen, ¢ A Atte goa Ht Atesg 5 As ja Hb Me ji + Mp 


+ igs, 541 HF Ui1, 541 — 20u,; = 0. 


It is easily verified that a relaxation factor r of the form (2) leads in the case of 
(13) to the hyperbolic partial differential equation 


(14) Buzz a Sigg = Que a 3ttyt a 5Cuz, 


similar to (4), and that the best choice for the parameter C here is C = 13+k,/5, 
while the corresponding value for the exponent » governing the rate of converg- 
ence of the iterations is p = 6k,/13+. Thus we see that these quantities are slightly 
larger for the nine-point equation than they are for the five-point equation, a 
result not inconsistent with the numerical data compiled by Young [4], which 
indicated that the optimal value of r decreases and the correspondingly smallest 
number of required iterations changes insignificantly when we pass from the five- 
point to the nine-point interpolation formula. The strong dependence of the latter 
formula for the Laplace operator on the particular partial differential equation 
involved appears to have obscured until now the explanation of this phenomenon. 
At all events, for the nine-point difference equation the asymptotic expression 
for the best value of the relaxation factor r is 


2 2 


= < . 
on 1 + 134kyh/S ~ 1 + 3.074h/A} 





Finally, our introduction of the second order hyperbolic differential equation 
(4) indicates an analogy between the successive overrelaxation scheme and the 
second order Richardson method, which is based on a more readily accessible 
hyperbolic equation (cf. [1]) and converges at a comparable speed. The advantage 
of successive overrelaxation is that it requires retention of less data and is there- 
fore more adaptable for use on high speed computing machines. 
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Uniqueness of the Projective Plane of Order Eight 


1. Introduction. A finite projective plane with m + 1 points on each of its 
lines is said to be of order m. Such planes are known to exist whenever is a 
prime or prime power, there being at least the Desarguesian plane coordinatized 
by the finite field with » elements [8]. For some time it has been known [9] that 
for order 9 there are some non-Desarguesian planes. It is an easy exercise to 
construct the planes of orders 2, 3, 4, and 5 and to verify that they are unique. 
The results of Bruck and Ryser [1] include as a special case the nonexistence of 
planes of order 6 although this special case had been shown earlier by Tarry [7 ]. 
The combined results of Pierce [5] and Hall [2] show the uniqueness of the plane 
of order 7. In this paper we show the uniqueness of the plane of order 8, thus 
showing that 9 is the smallest order of a non-Desarguesian plane. 

The demonstration of the uniqueness of the plane of order 8 is based on Norton’s 
complete list of Latin squares of order 7 [4], there being a total of 147 varieties 
of these, an omission in Norton’s list having been found by Sade [6]. By a theo- 
retical argument it is necessary to consider only 100 of these, and machine 
calculations, described in section 4 of this paper carried out on SWAC at the 
University of California, Los Angeles, went far enough so that it was not difficult 
to complete the work by hand. The completion of the search is given in section 5 
of this paper. 

2. Theoretical basis for calculation. Let A, B, C be the vertices of a triangle 
in a plane of order n. Call AB the line at infinity L,, AC the line x = 0, and BC 
the line y = 0. Label the » — 1 remaining lines through A as x = 1, x = 2, ---, 

= n — 1 in any order and also the m — 1 remaining lines through B as y = 1, 
+++, y =m — 1 in any order. A point P not on L,, will then lie on a unique line 
x = a and a unique line y = b. Then assign to P the coordinates (a, d). 

The m — 1 lines through C = (0, 0) apart from AC and BC will intersect each 
of x = 1, ---,x = m — 1 once and each of y = 1, ---, y = m — 1 once. Such a 
line L will intersect L, in some infinite point and will also contain (0,0) and 
n — 1 points (7, 7) where 7 and 7 take values 1 to m — 1. With L associate the 
As ') it (0,0), (1,01), «++, (» 1, a1) are the finite 
0, a1, «**, On_1 
points of L. The m — 1 different L’s yield permutations: 


permutation ( 


au, 12, oes “Chai 

21, 22, oes Ca e—t 
(2.1) 

Qn-1,1; On-1,2, °***, On—1n-1 


and the array (2.1) will form a Latin square of order m — 1. Since a line through 
C intersects each of y = 1, ---, y = m — 1 once, each row of (2.1) contains each 
of 1, ---,#” — 1 once. Since a line x = i intersects the m — 1 L’s in different points 
the z-th column will contain each of 1, 2, ---, — 1 once. The three constraints, 
row, column, digit correspond respectively to a line through C, a line x = 7 and 
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a line y = j. The interchanging of these three constraints corresponds therefore 
to an interchange of the roles of the vertices A, B, C in the construction of het 
Latin square. A permutation of the rows is merely a permutation of the order in 
which the lines through C are taken. A permutation of the columns corresponds 
to a relabeling of the lines x = c through A and a substitution on the digits corre- 
sponds to a relabeling of the lines y = c through B. Thus, these equivalent forms 
of Latin squares obtained by permuting rows, permuting columns, and substitut- 
ing on digits, amount to no more than relabeling the lines of the three-net of 
pencils through the vertices of a triangle A B C excluding the sides. And an inter- 
change of constraints corresponds to altering the roles of the vertices A, B, C in 
obtaining the square (2.1). It is immediate that any Latin square of order m — 1 
may be used as an array (2.1) corresponding to a three-net of the pencils through 
three points A, B, C excluding the sides. It is not immediate nor even true 
that any Latin square yields a three-net which can be extended to a complete 
plane. (However, for countably infinite planes this and more is true.) [3]. 

Norton [4] has given a list of the Latin squares of order seven to within the 
equivalences of permuting rows, columns, substituting digits and interchanging 
the constraints. Norton lists 146 varieties. Sade [6] found an omission and 
verified that with this additional variety included the list is complete. 

These 147 squares form the starting point in our search for planes of order 8. 
We seek to complete these by adding as many as possible additional lines not 
through A, B, or C. It is however sufficient to treat only 100 of the 147, effecting 
a saving of approximately one-third of the machine time. Norton lists his squares 
according to the number of intercalates which they contain, an intercalate being 
a sub-array of the type: 

a¥--8 
(2.2) 

Ewe ft 
If the columns of the intercalate (2.2) are the i-th and j-th then the four points 
(z, a), (4, 6), (7, a), and (j, 6) are such that no three lie on a line and so they form 
a quadrilateral. The diagonal points of this quadrilateral are the points A, B, C. 
Conversely if A, B, C are the diagonal points of a quadrilateral these four points 
yield an intercalate in the array (2.1). Now in a plane of order 8 there are 
73-72-64-49/1-2-3-4 quadrilaterals and there are 73-72-64/1-2-3 triangles. 
Now for some quadrilaterals the diagonal points may lie on a line. (Indeed in 
the Desarguesian plane of order 8 this is always the case.) Thus the triangles of a 
plane of order 8 are diagonal points of quadrilaterals at most 73-72-64-49/1-2-3-4 
times and as there are 73-72-64/1-2-3 triangles there must be a triangle A, B, C 
which yields a square with at most 49/4 intercalates and thus at most 12 inter- 
calates. Hence in constructing a plane of order 8 we may start from a square of 
order 7 with at most 12 intercalates. Such squares are numbers 1 through 99 of 
Norton’s list and the omission found by Sade. 


3. Completion of squares to the full plane. The one hundred squares of order 7 
used in the search for planes of order 8 could all be normalized so that the first 
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two lines read, adding an initial 0 for C = (0, 0), 


01234567 


G4) 02345671. 


It is sufficient to deal with the finite points (7, 7) 7, 7 = 0,1, ---, 7 constructing 
an affine plane of order 8 since the completion to a projective plane by adding 
infinite points and L,, is trivial. The pencil through (1, 1) will include the lines 
x* = 1,y = 1 and the line 0 1 2 3 4 5 6 7 of (3.1). The remaining lines through 
(1, 1) will go through the points (2, 3), (3, 4), (4, 5), (5, 6), and (6, 7) of the second 
line of (3.1) and there will also be a line through (1, 1) parallel to the second line 
of (3.1). These six lines will thus have the appearance: 


X13XXXXX 


Xi X4XXXX 
(3.2) XU XXS5SXXX 
X1iXXX6XX 
Xi’XXXXTX 


XUXXXXX X. 


Here X’s stand for digits to be filled in. Similarly when the six lines of (3.2) have 
been added to the lines given by one of the 100 squares there will be five more 
lines of the pencil through (2, 3) whose appearance will be: 


XX3X4XXX 
XX3ZXX5XX 
(3.3) XX3XXXG6X 
XX3XXXX 7 
XX3IXXXX X. 


Earlier calculations of the lines of (3.2) were made on the Engineering Re- 
search Associates 1101 Computer. For each of the 100 squares a number of 
possibilities for (3.2) was found, there being between 20 and 25 for each of the 
100 squares on the average. In this search the 1101 was used to perform somewhat 
complicated sorting procedures. Of the 7! = 5040 permutations of the form 
X1X XX X X X there are 1478 which are consistent with (3.1) in that they 
do not agree with either line of (3.1) in as many as two places. These are dis- 
tributed as follows: 


Type AY’tTI RR Ana A a TS 
X1X4XX X X — 140 
X1XXS5XXX-— 141 
X1XXX 6 X X — 140 
XTAa PAA TRESS 

Parallel to second line X 1X X X X X X — 787 


The set of 5040 permutations was easily built up on punched cards and the 
1478 consistent cards were sorted from these. These 1478 permutations were 
taken as input for the ERA 1101 and were handled in two stages. In the first stage 
permutations were eliminated which were inconsistent with any one of the five 
remaining lines of a particular square and the consistent lines were punched as 








oO —_ = VS VN 


re 


the 


age 
ive 
| as 





UNIQUENESS OF THE PROJECTIVE PLANE OF ORDER EIGHT 189 


output on tape. For each of the squares the number of consistent lines was ap- 
proximately 115 to 120. This was a fast code, most of the time being taken up by 
punching outputs, and averaged about three minutes for each square. The second 
stage involved taking the 115 to 120 consistent lines for each square and forming 
the six lines (3.2) from these in every possible way. This was relatively slow taking 
about half an hour for each square and yielding, as remarked above, about twenty 
solutions for each square. In every case tried, except for the Desarguesian plane, 
only three or four further lines could be added. But working by hand it took at 
least an hour to test the construction of further lines given a square and the 
lines of (3.2). 

4. Description of the SWAC code. Making one hundred separate runs SWAC 
took each 7 X 7 square with no more than twelve intercalates and constructed 
the six lines of (3.2) and the five lines of (3.3), punching them as output whenever 
all eleven were found. The method employed was to try, for each of the X’s in 
(3.2) and (3.3), all possibilities 0,1, ---, 7 compatible with previously assigned 
X’s. To specify the criteria of compatibility let P;(j) be the j-th element of the 
i-th row of the 7 X 7 square augmented by an initial column of zeros, and P;’ (j) 
the similar element of (3.2). 

For any (i, j) referring to an X in (3.2) all P,’() with k <i or with k =i 
and / < j will be assumed already assigned. The restrictions on P;’(j) are then 


(a) Pi’(j) #1, 

(b) P/(j) #i+2 if i < 6, 

(c) Pi (j) ¥ Pi (D for any / < j, 

(d) P,’(j) # P,' (7) for any k < i, 

(e) If for any / <j or 1 = 2 or 1 = i+2 we have P,/(l) = P,(J), then 
Pi'(j) # P, (3). 


The first three conditions insure that each row of (3.2) is a permutation. The 
last two prevent any such row from agreeing in two positions with any previous 
row of (3.2) or any row of the 7 X 7 square. 

For an element P,’’ (7) of (3.3) we have conditions like these plus an additional 
condition similar to (e) but with P replaced by P’. In the following discussion we 
shall consider only the assignment of the P,’ (7); the case of the P;’’(j) is more 
complicated due to the extra condition but is basically the same. 

To apply these criteria we construct the set of allowable values of P;(7). If 


saute tts thu tprou<s. 


then any P,’ (j)eS* (the complement of S*) will satisfy (a), (b), and (c). Similarly 
(d) and (e) will be satisfied if P,’(j)eS;, where 


S; = {Pi'(j)|k < i} U {Pas |Pe(2) = 1 or Pe + 2) =i +2 
or P, (J) = P;' (2) for some 7 < 7}. 


(Because of (3.1) the conditions P,(2) = 1 and P(t + 2) = i + 2 imply k = 1 
and k = 2 respectively. This convenient simplification is not essential.) Hence 
S = §* (1) 8; is the set of values of P;' (7) compatible with those already assigned. 

To select a P;'(j) one computes S*, S;, and S. Suppose first that S is not 
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empty. Then P,’(j) is taken to be the smallest element of S, and one moves on to 
the next unfilled position in (3.2). If this position is in the same row S* is replaced 
by S* U {P//(j)}, j replaced by 7 +1 (or 7 + 2 if 7 = i + 1), and the new S; 
calculated. On moving to a new row 7 is increased by 1, 7 put equal to 1, and S* 
set equal to {1,7 + 2} if ¢ < 6 or {1} if ¢ = 6. When (3.2) is completely filled 
one starts a similar process in (3.3); when this is filled one records all the P;’ (7) 
and P;’’(j). 

Suppose now that for some (7, 7) the corresponding S is empty. This means 
that the previously assigned P;’(j)’s must be changed. One therefore ‘“‘back- 
tracks,” returning to the previously considered P;(j) by reducing j or, if necessary 
to move up to the preceding row, by reducing 7 and putting 7 = 8. In the former 
case one gets back the old S* by removing the element P;’(j); in the latter one 
puts S* = {0,1, ---, 7}. The old S; could be re-computed, or it could have been 
saved from the time of its previous use. Since the calculation of S; is somewhat 
complicated, and since there is much back-tracking and advancing within a row 
but not much between rows, it was decided to save the S; within each row but to 
re-compute them on back-tracking to another row. 

After back-tracking and re-computing S one must of course avoid selecting 
from S the same P;’(7) that had been used before. This is accomplished by re- 
moving from S all elements less than or equal to the P,;’(j) that was last used. In 
this way the successive selections of P,’(j) for a given (i, j) are strictly increasing 
until all possibilities are exhausted and one backtracks still farther. 

After a solution is recorded the process is continued by the standard back- 
tracking procedure, just as if one had proceded to a next step and found S = ¢. 
The whole process stops when all possibilities for P,'(1) have been exhausted. 

In applying a routine of this sort it is obviously advantageous to eliminate 
many possibilities early in the process. In particular, in each row of (3.2) it is 
helpful if S; is as large as possible for small j. In the set-up as described 


Si = {P,’(1)|& < ¢} U {0}. 


However if we make a cyclic permutation of columns so that the first element 
of each row becomes the last we get 


Si = {P,'(1)|k < 4} U {1, 2}. 


This slightly better situation was used in the coding. 

In coding this process for SWAC the significant factor was the small (256 
word) high-speed memory (HSM). The program had to be stored on the magnetic 
drum and appropriate parts transferred to the HSM when needed. Fortunately 
it was just possible to squeeze into the HSM the whole procedure for handling 
the P,’(j) and also for the P;’’(7). The delays due to drum transfers entered only 
in the transfer from one of these routines to the other, and in the negligible number 
of input and output processes. Elements of a set were represented by positions 
within a word; i.e., the subset {a1, a2, ---} of {0,1, ---, 7} was represented by 
the number >> 2-*-'. This representation makes it easy to code the set operations 
of union, intersection, and picking the smallest element. 

The general flow of the code may be described by considering fourteen sub- 
routines as follows: 
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Possible Criterion for Successor 
Routine Predecessors Operation Performed Successor Yes No 
) A (G:) (New run) Input given square, cata- — B, B, 
| ogue it by element and 
| position 
B, A, Dy, Fi, Ge Compute possible values Are any values avail- Ga sz 
available for new posi- able? 
| tion 
Ci B, Store first available value —_ D, D; 
; in proper position 
D, Ci Advance position Is first set of lines B. By 
4 complete? 
4 E, B, Test whether backtrack Is a position open in Fi Gi 
7 is possible which to back- 
track? 
. F, E, Backtrack. Go back to —_ B, 2B; 
1 nearest available posi- 
t tion and recompute 
Gi Ey Halt or input new square — (A) (A) 
v B; Di, Dz, F; As in B, As in B, C2 E; 
D C; B, As in Cr — D: D: 
Dz C2 As in D; Is second set of lines e B; 
complete? 
B E; B, Test whether backtrack Is backtrack in last F, Gs 
i would carry back into five lines possible? 
n first six lines formed 
F, E: Backtrack as in F; — B. Bz 
8 G2 BE, Backtrack to next to last — B, B, 
position of first six lines 
™ P D: Punch 11 lines — E, E; 
b. 
e 
is 
nt 
56 
ate 
ly 
ng 
ly 
er 
ns 
by 
ns 
ib- 
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5. Results of calculations. For only one of the hundred 7 X 7 squares was it 
possible to construct the six lines of (3.2) and the five of (3.3). This was Square 
Number 1: 


ooocooco 
SAM WD 
= ITA & bd 
Ne NTA UP Ww 
On KR TO Hn > 
FPwWN FR ATA Ww 
Or DN RK 81 
Aun & WN = 


, 


and four sets of eleven lines were found for this. The computing time for a square 
varied from eight to fifteen minutes, averaging about twelve minutes. 
The four ways of adding the eleven lines to Square Number 1 are: 


— 
NS 


“SNe POA MH 1 WwW PA WM dO 
Ur ONO Re ee 
WOwoWWW UP OO W 
AonerensnN A T10 Fw 
NArFY OR ACH UW 
PHOTO WN A & =O 
Fam aAnnNnNOON UW DN 
on UP NP UWDN OC A 
rA MT FA U1 & dO 
Aur ON KR RE RR ee 
RWW WWA NF SO = WH 
NONF AON Ak w 
CONF ARP WOON UWA = 
= FP wT Db AT QW A SD OC 
UF ANON OC WH WN > 
NN Or UU PWN OO A 


w 
c= 


SPUR DAA WwW Ps bd 
aac rf aT RRR Re ee 
Owwowuana co & = Ww 
NuwmrePsTONTUN OK DA 
Arent nN RP WH TD © 
KF OP UN OrR AWN ™ 
FnNOA Cre FP AIT®WNH SO WU 
SON N AUN OWA WS 
aAarnsT PW DA = bY 
a a ns ee ee ee 
COowewoeawarns Oo Uw 
Ne anacoauwnn ™ 
NIONF RP ATWO UN A 
CONF UAT WHA” 0 
KF UO ArPD Usb FF WwW © 
FPAITON AN OH WAM 


Of these the first two cannot be completed to planes while the third and fourth can 
both be completed uniquely to planes, the completion being the unique Desargue- 
sian plane in both cases. Indeed the second is equivalent to the first so far as 
completion is concerned since if we interchange coordinates, replacing (x, y) by 
0123456 ) the 


(y, x) and then renumber according to the substitution ( 01642753 
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first six lines of the second case become the first six lines of the first case, while 
this change leaves the lines of Square Number 1 unchanged. Thus if either of the 
first two can be completed then both can. But in fact neither can. In the first 
case no line X X X 6 X 7 X X can be added while in the second no line 
XX7XX6X X can be added. 

In the fourth case if we interchange coordinates, replacing (x, y) by (y, x) and 

01234567 

0147362 ae 
Square Number 1 is unaltered and the first six lines of the fourth case become the 
first six lines of the third case, whence these two cases must lead to the same 
completions. The completion is unique and is the Desarguesian plane. In the 
third case this may be written: 


then renumber all coordinates according to the substitution 


ann ort acre 
orrh Us OD & 
naar Oo = Ww 
AnrnNOF UW > 
KB WoOoNnarnu 
ROWAN HE NA 
COPFP AN WU 


IANS WN HK © 


along with six other 8 X 8 squares whose top rows are the remaining lines of 
Square Number 1 and whose columns are the same as those of the 8 X 8 square 
above. This finally shows the uniqueness of the plane of order 8. 

Although it has already been remarked that the first two cases above cannot 
be completed to full planes, it is a curious fact that a great many consistent lines. 
may be added. These are listed here: 


01234567 02345671 
10475326 36051427 
27046135 64273015 
75362410 71526304 
03456712 04567123 
12637405 32176540 
41072653 40321765 
67310524 61705432 
05671234 06712345 
21357046 14306257 
50216473 23570461 
73104625 51643720 
07123456 

31460275 

45736021 

56374102 
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Thus we have half of seven 8 X 8 squares. In projective terms we have a 
total of 45 lines with 9 points on a line which are consistent with themselves but 
which cannot be completed to the 73 lines of a projective plane. 


MARSHALL HALL, Jr. 
Ohio State University 
Columbus, Ohio 


J. DEAN SwiFt 
University of California 
Los Angeles, California 


ROBERT J. WALKER 
Cornell University 
Ithaca, New York 


. bonne gga was supported by the National Science Foundation and the Office of Naval 
esearch. 

1. R. H. Bruck & H. J. Ryser, “The nonexistence of certain finite projective planes,”’ 
Canadian Jn. Math., v. 1, 1949, p. 88-93. 

2. MARSHALL HALL, jr., “Uniqueness of the projective plane with 57 points,” Am. Math. 
Soc., Proc., v. 4, 1953, p. 912-916; Correction to * wees of the projective plane with 57 
points,” Am. Math. Soc., Proc., v. 5, 1954, p. 994-99 

. D. R. Huaues, ‘ ‘Additive and multiplicative ie of the planar ternary rings,” Am. Math. 
me ’Proc., v. 6, 1955, p. 973-980. 

4. H. W. Norton, “The 7 X 7 uares,’ ’ Ann. Eugenics, v. 9, 1939, p. 269-307. 

5. W. A. Prerce, “The impossibility of Fano’s configuration in a projective plane with eight 
oie r line,” Am. Math. Soc., Proc., v. 4, 1953, p. 908-912. 

SRS Sane, “An omission in Norton’s list of 7 X 7 squares,” Annals Math. Statistics, v. 22, 
1951, p. 306-307. 

7.G. Tarry, “Le probléme des 36 officiers,”” Compte Rendu de Il’ Association Frangaise pour 
l’ Avancement de Science Naturel,”’ v. 1, 1900, p. 122- 123; v. 2, 1901, P. 170-203. 

8. O. VEBLEN & W. H. BUSSEY, “Finite projective geometries,’ Am. Math. Soc., Trans., 
v. 7, 1906, p. 241-259. 

9. O. VEBLEN & J. H. MAcLAGAN-WEDDERBURN, ‘‘Non-Desarguesian and non-Pascalian 
geometries,” Am. Math. Soc., Trans., v. 8, 1907, p. 379-388. 


- 


On the Location of Gauss Sums 


We shall understand by a generalized Gauss sum of order k the sum 
p—l 
= > exp (2rim'/p), (p = kf +1, a prime). 
m=0 
This sum can be thought of as the principal root 2 of the reduced period equation 
of degree k for the so-called f-nomial periods, 2; = kn; + 1, where, as usual, 
und 


= x exp (2xig*’t*/p) G = 0, 1, >° a k— 1). 
v=) 


Since the remaining k — 1 roots of the period equation depend on the primitive 
root g, the singling out of one root 2 as the principal root is justifiable. 


For k = 2, it is well known that [1] 


Vp ifp=4n+1 


Ss. = 
() * lind if p = 4m —1. 
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In general it is known that 
— (k — 1)Vp < Si < (k — 1) vp. 


Fork = 3, and 4p 


L? + 27M?, the three roots of the reduced trinomial cubic 


F(z) 


2—3p2—plL=0 (L =1 (mod 3)) 


lie in the intervals 


(—2~p, —vp), (—~p, VP), (vp, 2p). 


Kummer [2] noticed that the principal root zo lies in these intervals with 
frequencies approximately of 1 to 2 to 3 for primes less than 500 and conjectured 
that these ratios hold in general. Further calculations by Goldstine and von 
Neumann [3] for primes up to 10,000 exhibited ratios of the order of 2 to 3 to 4. 
They conclude that ‘these results would seem to indicate a significant departure 
from the conjectured densities and a trend toward randomness.’ On the other 
hand a case has been made by G. Beyer [4] that these are indeed the true ratios 
since 2 + 3 + 4 = 3?, and since these densities hold for composite as well as 
prime moduli. It therefore seemed of interest to continue this investigation on 
the SWAC at the Numerical Analysis Research of the University of California 
at Los Angeles. The results of Goldstine and von Neumann covering the first 
611 primes were extended to the first 1000 primes of the form 6m + 1, with the 
result that the ratios have now become of the order of 3 to 4 to 5, thus continuing 
the trend away from the conjecture and towards equal distribution noted by 
Goldstine and von Neumann. The actual figures will be given below. 

Before doing this we would like to suggest another division of the primes of 
the form 6n + 1 into three classes with respect to the absolute value of S;. Since 


F;(—~3p) = F:(V3p) = —pL, 
there is always one root of F; = 0 which is greater than 3p in absolute value, 


and we already know that there is a root of F; = 0 which is less than yp in abso- 
lute value. Therefore the absolute value of S; lies in one of the three intervals 


(0, Vp), (vp, ¥3p), (V3, 2vp). 


We have actually counted the number of times that the principal root lies in 
these intervals using the tables of Goldstine and von Neumann and our own 
tabulation from the SWAC for primes above 10,000. We find that the number of 
primes in these three categories is approximately the same over the whole range 
covered. It is hoped that these figures will encourage some characterization of 
these classes of primes, which escapes the writer. The actual figures follow, where 
the primes are divided into 10 groups of 100 primes each. 
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TABLE 1. Number of primes p = 1 (mod 6) such that 
Group Si<—Vp —Vp<Si<vVp Si>Vb |Ss]< VP VP <|Ss|< V3p |Ss|> V3P 


I 18 28 54 28 36 36 

II 21 38 41 38 35 27 
Ill 21 33 46 33 34 33 
IV 29 32 39 32 32 36 
Vv 28 29 43 29 35 36 

VI 18 38 44 38 26 36 
Vil 34 26 40 26 31 43 
Vill 24 32 44 32 41 27 
IX 29 34 37 34 35 31 
x 19 31 50 31 38 31 
Total 241 321 438 321 343 336 


A similar problem has been raised in the quartic case for the sum 
p-l 
Si = D exp (2rim'‘/p) 
m=0 


by Hasse [5] in the closing pages of his book on the theory of numbers. The 
discussion of the problem, as given there, depends on much of what is in the book 
and an elementary treatment of it might be in order together with additional 
numerical evidence. 

It was already known to Lebesgue [6] that S,, as well as the companion sum 


p—1 


Si’ = 2 exp (2xirm*/p), 
m=0 


where 7 is any quadratic, but not quartic, residue of , satisfies the quartic 
equation 


(2) [2 + (1 -- (—1)°-942)p P — 4p(2 — a)? = 0, 


where p = a? + 40? and p = a = 1 (mod 4). 
We can write (2) in factored form as follows 


(3) [st — 2N pe + (1 — (-1)0-9"2)p + 2Vpa] 
x [st + 2vpe + (1 — (-1)-"2)p — 2vpa] = 0. 


Since by definition of S2, S, and S,’ and by (1) 
Si + Si = 2S: = 2Vp, 


it follows that S, and S,’ are roots of the first factor of (3) and hence 





(4) Sa; SS, = pt V2V(—1)@-DMp vai v pa. 
We consider with Hasse the Iccation of 


(S — Si’)/4 = N(p — Vpa)/2 forp =8n+1, «= +41 


(S) A= (S — Si’)/(44) = eV(p + Vpa)/2 forp =8n+5, «= +1. 
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Besides the obvious ambiguity in the sign of «, the sign of a is “known” only 
by reference to a table of quadratic partitions [7]. 
Hasse divides all primes of the form 4n + 1 into four categories as follows: 


neti ¥2p/2 <A < Vp (e = +1, |a| = —1 (mod 4)) 
i= for 
8" + 5 0<A < ¥2p/2 (« = +1, |a| = —1 (mod 4)) 


—1, |a| = +1 (mod 4)) 


oe : —2p/2 <A <0 (e 
= r 
” +41, Ja] = +1 (mod 4)) 


Sn+5  Vp/2<a<v 


ss —vVp <A < —y2p/2 (e€ = —1, |a|= —1 (mod 4)) 
5 = for 
8n +5 — 72p/2<A<0 (« = —1, |a| = —1 (mod 4)) 


nda 0<A < ¥2p/2 (« = +1, |a| = +1 (mod 4)) 
Ph 8n +5 os —vVp <A< —Y¥2p/2 (¢ = —1, |a| = +1 (mod 4)). 


These intervals become meaningful if one thinks of them as the projections 
on the x and y axis respectively of the first, third, fifth and seventh octants of 
the circle of radius ¥p. Hasse (there are two errors in Hasse’s table, which he 
attributes to Kaluza. The primes 677 and 877 belong to class p; and not 3, giving 
ratios 21:21:12:26 instead of 21:23:12:24) found that for primes less than 1000 
the four types of primes appear with frequencies which are approximately 2:2:1:2 
for the primes of the form 8” + 1 as well as for primes of the form 8” + 5. He 
conjectured that there are infinitely many primes in each class and that the 
frequencies for primes of the forms 8” + 1 are the same as for the primes of the 
form 8” + 5. 

Another, possibly more simple-minded, way of dividing these primes into four 
categories would be by grouping together all primes, irrespective of their form 
modulo 8 for which A lies in a given interval as follows: 


po for which — Vp < A < —V2p/2 (e = —1, |a| = —(—1)®-* (mod 4)) 
b2 for which —J2p/2 < A <0 (e = —1, |a| = (—1)®-* (mod 4)) 
psfor which 0 <A < ¥2p/2 (e = +1, |a| = (—1)®"* (mod 4)) 
Pefor which 26/2 <A < vp (e = +1, |a| = —(—1)®* (mod 4)). 


(6) 


If we now denote by (x) the number of primes in the class ~, not exceeding x, 
then obviously 


m1(x) + 25(x) = a1(x), the number of primes < x for which |a| = —1 (mod 4) 

3(x) + 47(x) = a3(x), the number of primes < x for which |a|= +1 (mod 4) 

wo(x) + #6(x) = ao(x), the number of primes < x for which |a| = —(—1)®-)4 
(mod 4) 


w2(x) + 2a(x) = a(x), the number of primes < x for which |a| = (—1)@-)4 
(mod 4). 
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It has been shown by Landau [8 ] that the Gaussian primes of the form a +0i, 
when divided into residue classes with respect to any modulus are asymptotically 
equally distributed. The Gaussian primes can be divided into two classes modulo 
4, according as |a|= +1, or —1 (mod 4), @ =1 (mod 4), or according as 
|a| = +(—1)®-“ (mod 4). It therefore follows that 


lim a@3(x)/a1(x) = lim ae(x)/ao(x) = 1. 


zo a 7nw 


Since the ratios observed by Hasse at x = 1000 give 
a1 (1000) /a3(1000) = 3/4, 


they cannot be the limiting ratios. 

It should also be noted that by (6), ao gives the number of primes x, for 
which |A|> ¥26/2, while a2(x) counts those for which |A|< 42/2; there- 
fore |A| is in the limit equally likely to be in the two intervals (0, ¥2p/2) and 
(¥2p/2, Vp). This could possibly be used to support the conjecture that the 
absolute value of S; is also equally divided between three intervals in the limit. 
However, the consideration of the Eisenstein primes of the form a + bw modulo 
three does not seem to shed any light on the cubic problem any more than the 
consideration of the Gaussian primes modulo four helps us to find the limiting 
densities over the four intervals as taken by Hasse or as given in (6). We have 
therefore turned to the SWAC once more in order to calculate A for all primes 
of the form 4” + 1 less than 10,000. 

The time required by the SWAC for all 2549 primes for which A was calcu- 
lated was about eight hours. Towards the end of the run, for primes in the neigh- 
borhood of 10,000 it took almost a minute per prime, but-considering that there 
were about 5000 sines or cosines to compute and add, not to speak of the calcula- 
tion of the appropriate arguments, the time does not seem excessive. Also the 
accuracy of about 6 decimal places in cases which were computed by hand by (4) 
seemed very gratifying. 

The results are as follows: 


TABLE 2 
x m1(x) a3(x) 5 (x) 7 (x) mo(x) w(x) a(x) ae (x) 
1000 21 21 12 26 17 5 25 23 
2000 42 34 27 44 23 26 47 41 
3000 54 55 43 59 46 © 46 62 57 
4000 77 68 49 75 54 55 84 76 
5000 87 86 60 96 67 74 98 90 
6000 105 99 71 108 78 84 112 109 


7000 125 116 79 122 86 94 134 128 
8000 145 126 93 135 101 106 151 141 
9000 163 136 102 153 116 114 163 161 
10000 179 151 113 166 125 131 176 177 


These results do not change appreciably the picture obtained for primes less than 
one thousand. The frequency 7; is still consistently low, but Hasse’s ratios of 
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2:2:1:2 now look more like 3:3:2:3, giving the ratio of 5/6 for a:(x)/as(x) for 
x = 10,000, instead of 3/4 at x = 1000, which is in line with the fact that this 
ratio tends to unity. The figures also seem to indicate that ro(x) and r2(x) are 
very close. The same can be said of 4(x) and 26(x) giving ao/az very close to 
one all along the line. 

The next case, Ss, has not been studied in detail before, although the quintic 
period equation which it satisfies is known [9] (There is a misprint in equation 
(10) of [9]. A correct expression is (8) of the present paper.) in terms of the 
following quadratic partitions: 


16p 


xw 


x? + 50u? + 500? + 125w? 
v— u—4uv, p =x =1 (mod 5) 


(7) 


and is as follows: 


(8) Fs(z) = 2° — 10p2* — Spxat + Sp[p — (x® — 125w*)/4]]s 
+ px — plx® + 625 (u? — v*)w]/8 = 0. 


The roots of this equation lie between —4Vp and 4p, but when divided by Vp 
they no longer lie one each in any five fixed intervals. However, we can show by 
straightforward algebra that there always is a root z such that |z| > 5p, and 
that there are one, two, or three roots in the middle interval |z| < yp. Figures 
obtained for the five roots of (8) for primes less than 10,000 on the SWAC show 
that all of these cases do arise and that for approximately half the primes there 
are two roots in the middle interval, while for about a quarter of the primes there 
are one or three roots in the middle interval. The majority of cases with one root 
in the middle has one root in each interval, but there are exceptions. Nevertheless, 
the frequency of |.S;| in these three intervals is surprisingly uniform, namely 


TABLE 3 
p< ISsi< Vp Vp <|Ss|<~vV5p V5p <|Ss| <4vp Total 
1000 15 10 15 40 
2000 30 19 24 73 
3000 37 33 33 103 
4000 46 43 45 134 
5000 55 50 58 163 
6000 67 60 67 194 
7000 76 71 79 226 
8000 84 76 86 246 
9000 91 88 94 273 
10000 108 97 101 306 


These figures seem to be in line with the conjecture that |S5| is equally likely 
to lie in any one of the three intervals 


(0, vp), (v2, V5p), (v5p, 4vp). 
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The distribution of S; itself over the corresponding 5 intervals gives at 10,000 
the following numbers from left to right 


33, 34, 108, 63, 68 


with ratios of the order of 1:1:3:2:2. 

Another possible division into intervals is provided by the scaled projections 
of the fifth roots of unity on the x-axis. This gives + (1 + ¥5) for division points. 
For p = 10,000 the distribution from left to right over the five intervals is 


12, 51, 122, 97, 24 


or ratios of the order of 1:4:10:8:2. Although the sum of these figures is tempt- 
ingly 25, our experience with the cubic and quartic cases does not permit us to 
indulge in the hypothesis that these are indeed the limiting ratios. It is also hard 
to see how to interpret these figures as a generalization of Kummer’s conjecture. 

We have pursued a step further the conjecture about the equal distribution 
of |.S,| for odd & in the intervals 


(0, vp), (vp, vVkp), (vkp, (k #1 1) Vp) 


by calculating on the SWAC the values of S; for p = 7n + 1 < 5000 with the 
following results: 


TABLE 4 
p< ISi]}<~vp Vp <|S:|<~vV7p 7p <|S;| < 6vVp Total 
1000 6 11 “10 27 
2000 10 21 13 44 
3000 16 32 21 69 
4000 23 38 27 88 
5000 33 43 32 108 


These figures are not quite as convincing as those for the quintic, but since 
nothing is known about the distribution of the roots of the seventh degree period 
equation or about the form of the equation itself it is hoped that the actual tabu- 
lation of the roots as given by the SWAC will be of some use. 

The actual values of these roots together with those of the quintic and the 
principal root of the cubic equation will be deposited in the UMT File of MTAC. 

We wish to express our thanks to the Numerical Analysis Research Group of 
the University of California at Los Angeles for the use of the SWAC and in 
particular to Ruth Horgan and John Selfridge for checking out the various codes 
and operating the machine. 

The code itself is diagrammed below. It presented no special difficulties and 
was put together from a standard cosine subroutine (a five term approximation 
of the 9th degree), which is used on the SWAC and fundamental number-theoretic 
subroutines, which will be described elsewhere. 
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1 
Program for the calculation of S, = 2 exp (2xim*) /p 





1. | Input last prine p= im+1 | 
1 ra 


2. | promo] cS 


3. Is p a prime? 























Yes | wo 
: 
k 
Se h Find m ies » | 


6 a Find — _/P | 


1. | Accumlate zP* cos 20 e,/o| 


. 


9. Is m < p? 











lee 

































































12. Does Ss, check? 











Yes } No 


(sa 








13. Output +P, Sn 


Box 3 uses the trial divisor test for primality. 

Box 5 uses a standard reduction modulo p. 

Box 6 uses the cosine subroutine. 

Box 11 provides an indication on the size of the sum for output. 

Box 12 consisted of either substituting the principal root into the equation, 
or adding all the roots to zero, in case all were computed. In a few cases, where 
the check failed, the calculation was repeated. 
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For roots other than the principal root the program was modified by replacing 
m* (mod p) by p’m* (mod p), (v = 1, 2, ---,& — 1), where p is a primitive root 
of p. In this case the primitive roots subroutine was incorporated into the routine 
and the k sums were computed abreast and added in 12 as a check. 


EmMA LEHMER 
Berkeley, California 
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A Method for Computing Certain 
Inverse Functions 


A method will be demonstrated for computing the inverse of certain functions. 
The method is applicable to the computation of logarithms and inverse trigo- 
nometric functions. It makes use of the binary expansion of real numbers and is, 
therefore, particularly suitable for use in automatic digital computing machines 
which use the binary number system. It is not recommended for hand computing. 

1. The Method. Let f(x) be a function which satisfies the following conditions, 


(i) f(x) is continuous and monotone on an interval (0, a] 
(including a, but not including 0), 
(ii) f(a@/2) is known, 
(iii) f(2x) and f(2x — a) can be computed when f(x) is known. 


(0, a] is taken to mean [a, 0) if @ is negative. Also, the symbol (f(0), f(a) ] 
will be taken to mean (f(0, +), f(@)] (f0,-),f@], Lf@, f(@, +)) or 
[f(a), f(0, —)), whichever is appropriate. 

Examples are: 


(a) f(x) = 2 a=-1 
f(2x) = (f(«))? f(2x + 1) = 2(f(x))* 
(b) f(x) = cosx a=f 


f(2x) = 2(f(%))? —1 f(x — x) = 1 — 2(f(x))? 


Let y « (f(0), f(a)], and let it be required to compute f(y), that is, to find x 
such that f(x) = y. The existence and uniqueness of such an x in (0,a] are 
guaranteed by condition (i). Let w = x/a. Then w is in the interval (0, 1]. It 
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follows that w has an expansion as a binary number, 


w= > u,-27 


n=l 


where u, = 0 or 1 for each m. Where two such expansions exist, choose the one in 
which there are infinitely many u,’s equal to 1. For each non-negative integer k, let 


° 
uy = > Unye° 2, 
n=l 


xX = wa. 
Then x» = x and each x; is in the interval (0, a]. Also, for each k, 
Xe = Fey + FX K41. 

Thus, 

{‘ if x, ¢ (0,a/2] 

Urq = . 
1 if xe (a/2,a}. 

This is equivalent to 
O if f(xx) « (FO), f(@¢/2)] 
1 if f(xx) e (f(¢/2), f@]. 
By condition (iii), it follows that f(x,4:1) can be computed when f(x;) is known. 
Since f(xo) = y is given, an inductive procedure is established for determining 
the u,’s. These determine w, and the required number x is wa. The procedure will 


be illustrated in three examples: 
Example 1. Compute log: (.6). In this example, 


Mey = 


f(x) = 2? a=-1 
f(2x) = (f(x))? f(2x + 1) = 2(f(x))? 
(f(0), f(@/2)] = (1, V-5] (f@/2), f@] = (V5, .5]. 


Since we have to square each f(x,) anyway, we can test whether f(x,) is in 
(0, f(a/2)] or (f(@/2), f(@)] by squaring and comparing with .5. Thus, each 
step consists of squaring f(x,), comparing with .5 and recording m4; = 1 if 
(f(xx))? < .5, and t41 = O if (f(xx))? > .5. Then we record f(xi41) which is 
either (f(x,))? or 2(f(x:))?, whichever is between .5 and 1. 


k f(x) (f(xx))? Mn+ 
0 6 36 1 
1 72 5184 0 
2 5184 .26873 856 1 
3 .53747 712 .28888 16545 23494 4 1 
4 .57776 33090 46988 vee 


Thus w 


.1011 --- (binary notation) and log: .6 = aw = —.1011 --- (binary) 
—(11/146+¢) 0<« < 1/16. 
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Example 2. Compute arc cos 0.5. In this example, 


f(x) = cos x a=f 
f(2x) = 2(f(x))? — 1 f(2x — r) = 1 — 2(f(x))? 
(f(0), f(@/2)] = (1, 0] (f(a/2), f@] = (0, —1] 


At each step we record m4: = 0 if f(x) > 0, wey: = 1 if f(x) <0. Then 
f (%e41) = & (2(f(xx))? — 1), the sign + is that of f(x;). 


k Ff (xx) Uk+1 2(f(xx))? — 1 
0 5 0 —.5 
1 —.5 1 —.5 
2 s 0 —.5 
3 —.5 1 —.5 


w = .0101 --- (binary), 
w= 1/3 
arc cos 0.5 = 2/3. 


Example 3. Compute arc cos 0.2. 


k f (%x) Unt 2(f(xe))? — 1 
0 3 0 — .92 
1 —.92 1 .6928 
2 — .6928 1 — .04005 632 
3 .04005 632 0 . tee 

w = .0110 --- (binary) 


w= 3/8+¢,0<e< 1/16 
arc cos 0.2 = wr = 3x/8 + ex, 0 < ew < 0.1964. 


2. Error Analysis. Condition (iii) is seldom realized in practice, for required 
quantities are usually ‘computable’ only approximately. Even so simple a 
process as squaring cannot be repeated indefinitely without eventually ‘‘rounding 
off.”” A more practical condition than (iii) is 


(iiia) there exist small positive numbers 6 and e, such that when f(x) is given, 
(x € (0, a.]) numbers can be computed which differ from f(2x) and f(2x — a), 
respectively, by less than «, and such that for any pair x, x’ « (0,a], 
| f(x) — f(x’)| < € implies |x — x’| < 6. 


Hereafter « will designate the maximum roundoff error. Suppose that the 
calculation of x proceeds as indicated above, except that at each stage the true 
value of f(x.) is replaced by an approximation y,, computed in accordance with 
(iiia). Then for each n, 


lyn — f(2f*(ya-1) — Gttn)| < € 
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and this implies 


\f* (yn) — 2f7*(n-1) + au,| < 8. 
Now, let 


25> x 
and for each positive integer N, let 
N 
ty = @ D> Un: 2 + 2-* f(y). 
n=1 


That is, 2y consists of the first N binary digits of w plus an error term. Then, 
for each n, 


Zn — 22-1 = 2-* (f- (yn) 2 2f (ya) + aU») 
| Zn — Zn—1| < 2%. 
Thus for every V 


|tvy — x| =|2v — 20| 
N N 
=|>° (fn — 2n-1)| < D> 2-5 < 5. 
n=l n=l 


Thus each z, differs from the required x by less than 6. If the process could be 
carried on indefinitely, therefore, roundoff errors would introduce a total error 
no greater than 6. Since 6 is a function of the computation facilities available and 
independent of the number of steps executed, it will be designated as the maximum 
residual error. It is interesting to note that the maximum residual error is nu- 
merically equal to the maximum error that could result from rounding the input 
number, y. 

If the process is terminated after NV steps, the result is an underestimate and 
should, therefore, be ‘‘rounded up” by adding 2-*—'a. The computed value for 
x is, then, 

N 
att > u,-2"}+ 2-74 


n=) 
This differs from zy by 
2-"(f"(yw) — a/2) 


which is not greater in absolute value than 2-*—a. Therefore, 2-*-"|a| is desig- 
nated the maximum termination error. The total error is obviously less than the 
sum of the maximum residual error and the maximum termination error. 

In Example 1 above (log: y), the maximum residual error is 1/2(2* — 1) 
where ¢ is the maximum roundoff error. It is less than 


€ log.2 
1/2 1 a log.2 


The maximum termination error is 2~*—, 
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In Examples 2 and 3 above, arccosy, the maximum residual error is 
arc cos (1 — «), which is approximately ¥2e, and the maximum termination error 
is r/2%+, 

3. Programming for Machine Computation. A typical machine program for 
the computation of f-'(x) by the method of section 1 would be as follows: 


(It is assumed that f(x) is an imcreasing function of x. If it is a decreasing 
function, then the relation > should be replaced by < in Step 1.) The symbol 
C,, means “‘contents of cell m.”’ 


Cell Number Initial Contents Representing 
1 y Yn 
2 1/2 2-1 
3 0 > u;° 2-‘ 
i=1 
+ > aed 2" 
5 f7(@/2) f@/2) 
Program: 


Step 1. If Ci > Cs, proceed to Step 3. If not, proceed to Step 2. 
Step 2. Compute f(2f-'(C:)) and store in cell 1. Proceed to Step 5. 
Step 3. Compute f(2f-'(C:) — a) and store in cell 1. Proceed to Step 4. 
Step 4. Store C2 + C; in cell 3. Proceed to Step 5. 
Step 5. Store 1/2C: in cell 2. Proceed to Step 6. 
Step 6. If C2 < C., proceed to Step 7. If not, proceed to Step 1. 
Step 7. Store C2 + C; in cell 3. Store a-C; in cell 3. 
The computed value of f(y) is Cs. 


Of course, the program outlined above can be improved when special prop- 
erties of f(x) are known. If, for example, f(x) is monotone in the interval (0, 2a] 
then comparing of f(x) and f(a/2) is equivalent to comparing f(2x) and f(a). 
(This condition is satisfied by 27, but not by cos x.) If in addition, f(x — a) can 
be computed when f(x) is known, the program outlined above can be modified by 
placing f-'(a) rather than f-'(a/2) in cell 5, and altering steps 1-3 to read: 


Step 1. Compute f(2f-'(C,)) and store in cell 1. : 
Step 2. If C: > Cs, proceed to Step 3. If not, proceed to Step 5. 
Step 3. Compute f(f-'(Ci) — a) and store in cell 1. Proceed to Step 4. 


The modified program has several advantages. For one thing, Step 3 of the 
modified program is likely to require fewer orders than Step 3 of the original 
program. In the case of the function f(x) = 2*, for example, Step 3 of the modified 
program is multiplication by 2, whereas in the original program, it is squaring 
and multiplying by 2. Another advantage in the case f(x) = 27 is that f(a) = 1/2 
is more easily and accurately evaluated than f(a/2) = 1/2. 

It may be possible to take advantage of special features of the particular 
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machine to be used, by combining certain steps. A program for log: x has been 
developed for the CRC-102A Computer, which has only four commands in the 
repeating loop, and only three additional cells referred to in the repeating loop. 
The repetitive part of the program is shown below, coded in octal, with explana- 
tory notes at the right. 


Cell Number (Initial) Contents 
2000 25 2007 2007 2007 Square y 
2001 27 2006 2005 2006 Shift logical 
2002 31 2006 2007 2006 Scale factor 
2003 37 2006 2003 2000 Test overflow 
2004 34 3000 2100 —* Exit 
2005 00 0000 0000 0001 Shift control 
2006 37. 7777 7777 7776 See text 
2007 ¥ Yn 


* Address of next command in main program. 


The number in cell 2006 is 


ans -4 tr 2) —1. 
i=1 


4. Other Inverse Trigonometric Functions. Although the method given above 
is applicable directly to the computation of other inverse trigonometric functions, 
the functions f(2f-'(y) — ua) are in general more complicated than + (2y’? — 1). 
It is advantageous, therefore, to evaluate other inverse trigonometric functions 
by converting them to inverse cosines, as follows: 


arc sin y = 2/2 — arc cosy 


i— 2 
+ 1/2 arc cos ( + =) 


(sign is that of y) 





arc tan y 





i-y 
arc cot y = 2/2 + 1/2 are cos (1 ) 


(sign opposite that of y) 


arc sec y = arc cos (1/y) 


arc csc y = 5are cos (1/¥). 


The application to arc sin y provides an interesting algorithm for computing 
arc sin y, which may be useful even in hand computing. The algorithm is as 
follows : 
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1. Iterate the function 2y? — 1, starting with the number y whose arc sin is 
required. 
2. Record the signs of the iterates in order. 


3. Accumulate the signs; that is, record the ‘“‘partial products” of the signs 
in order. 


4. Write descending powers of 2 between the signs accumulated. 
5. Multiply the series obtained by 2/2. 


Example: Compute arc sin ¥.75 


1. 75, 5, -5, —.5, —.5 
2 + + _ es _ 
\ a + = + _ eee 
4. +1/2 8/6 -1/8. 41/16 =1/32 --. = 3/3 
5 a/2-(2/3) = x/3 
arc sin ¥.75 = 2/3. 


5. Comparison with Other Methods. The usefulness of the method described 
above as compared with other methods depends, of course, on the function to be 
evaluated and on the features of the machine to be used. 

The fact that each iteration yields exactly one binary bit may be an advantage 
or a disadvantage; a method where error decreases faster than 2-* will converge 
with fewer iterations than this one. On the other hand, one iteration of this 
method may consist of fewer commands than an iteration of another method. 
The logarithm program for the CRC-102A described above has 4 commands per 
iteration as compared with 14 commands per iteration in another program for 
logarithm on the same machine. The program for arcsin by the method described 
above has 9 commands per iteration as compared to 22 for another arcsin pro- 
gram. The fact that each iteration yields exactly one binary bit also simplifies © 
error analysis, for the number of digits of accuracy is exactly one less than the 
number of digits computed. 


D. R. Morrison 
SANDIA Corp. 


Albuquerque, New Mexico 


A Method for Solving Algebraic Equations 
using an Automatic Computer 


Introduction. Many methods have been developed for solving algebraic equa- 
tions and several of these have been used with automatic computers [1, 2]. Those 
methods which are most suitable for use with automatic computers are ones 
which apply to a wide class of equations and which are relatively rapid when the 
degree of the equation is large. The method described here has been constructed 
with these considerations in mind and has been programmed for the ILLIAc 
computer at the University of Illinois. 
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A process is to be constructed to find m solutions to the general algebraic 
equations of mth degree 


(1) f(x) = age* + ax*" + --- +a, =0 


where the coefficients do, a1, ---,@, are complex numbers and a») # 0. Each root 
is found by an iterative procedure. Successive iterations toward a particular root 
are obtained by finding the nearer root of a quadratic whose curve passes through 
the last three points. The quadratic will in general have complex coefficients and 
complex roots. This solution is accomplished by a variation of the standard 
quadratic formula. Although the method derived here is rather complicated, no 
evaluation of derivatives of f(x) and only one evaluation of the polynomial f(x) 
is required per iteration. If the degree of the equation is large a greater amount 
of time is spent evaluating the function than is spent in the remainder of the 
process. Thus, the time spent per iteration is less with this process than with 
iterative schemes which require the calculation of derivatives, whenever the degree 
of the equation is large. 
The Lagrange interpolation formula will yield a quadratic 


(2) Li(f(x)) = box* + dix + be 


whose curve passes through the last three points (x;, f(x,)), (xis, f(xi-1)), 
(xs-2, f(x:-2)) where the coefficients bo, 1, b2 satisfy the system of equations 
box? + bx; + b2 = f(x), 

(3) box? + bins + bd: = Sf (xi-1), 

box?s_2 + dixie + be = f(xi-2). 
A somewhat more convenient representation for this quadratic is obtained by 
introducing the new quantities h = x — x;, hy = x; — Xi-1, Ais = Xi-1 — Xi-2, 
h = h/hy, 4 = hi/hy_1, and 6; = 1 + Ay. The Lagrange interpolation formula may 
now be written as the following quadratic in X. 
(4) Li(f(x)) = NOLS (x-2)AZ = Sf (xi-dard; + f(x)as] 

+ ABALS (ei2)AZ — f(xir)b2 + f(x) Ac + 6) ] + f(x). 
A single iterative step is obtained by letting x;,: be a value of x which makes 
L;(f(x)) vanish. We may solve the quadratic equation in \ obtained by setting 


expression (4) equal to zero. We then obtain A = Agi = (X41 — %5)/ (Hi — Xi-1) 
by using the inverse of the standard quadratic formula 





(5) . =4 = 2f (x8; 
oF git Vg? — Af (x)bALSf (xH-2)ai — f(x): + f(x] 


where g; = f(x:-2)A2 — f(xi-s)b2 + f(xi) (Ai + 6). From Aiyz1 we may obtain 
Meigs = Aagahs and Xi41 = xX; + Ais. This x:41 represents a zero of the quadratic 
described above. The sign in the denominator of (5) is always taken so as to make 
the denominator have the greater magnitude. This makes A;,; and h;4; each be 












| 
; 
| 
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that one of the two possible choices having the smaller magnitude so that x;41 
is the root which is closer to x;. 


A convenient starting method for this process uses artificial starting values 
at x» = —1, x; = 1, and x, = 0: 


Qn — An-1 + Qn-2 is used for Ff (xo), 


(6) Gn + Gn_1 + Gn-2 is used for f(x:), 
an is used for f(x). 
thus making A. = —4 and hk, = —1. This choice of starting values makes 


Lo(f(x)) = dn + Gn—1¥ + Ano? 


which approximates to f(x) in the neighborhood of x = 0. The advantage of 
this starting process is that it requires no special evaluations of the polynomial 
f(x) and is therefore rapid. 

Convergence of the Process. A final value of the root x; is taken when 
|x; — x:-1|/|x;| becornes less than some preassigned number. Such a convergence 
criterion is consistent with the use of floating point arithmetic in the calculation. 
As a result of this criterion we see that convergence occurs if x;1 = x;. This 
means that before convergence no two consecutive iterative results will be equal. 
Furthermore, if x; = x;2 we have 6; = aberrant = 0 so by (5) Aw = O and 

t—1 ~ Ai-2 
Xiz1 = x; also giving convergence unless x; = 0. Thus in normal operation of the 
process we see that x;, x;_1, and x;_2 are distinct. 

As each root is found it may be divided into the polynomial f(x) thus reducing 
the degree of the polynomial by one. The algorithm for this reduction is the 
commonly used one : 


(7) a,’ = 704-1 + ai, (i = 0, i, 2, 2 -) 


where a,’ is the new coefficient to replace a; and r is the root which has just been 
found. We make a_; = 0. Errors introduced by this process will be reduced if the 
roots are eliminated in order of increasing magnitude. By always starting at the 
point x = 0 one will tend to find roots in roughly this order. 

No general proof of convergence in the large has been obtained for this 
process, but convergence can be shown to occur whenever the process leads one 
sufficiently close to a single or double root. 

In order to facilitate the study of convergence let us assume that x;4: = 0. 
This loses no generality since a simple shift of origin is always possible within the 
system. At point x;,; we also have L;(xi;1) = 0 so that 


(8) Li(%ig1) = Dov? inn + Di%i41 + be = D2 = 0. 


Now each of the functions appearing on the right hand sides of equations (3) 
may be expanded about x;,: = 0 as 


(9) f(%i-m) = z= x*; mf (0)/k! 
k=0 
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and 5, may be written as 
(10) bz = D buf®(0)/k! = 0 
ko 


where by, is obtained by solving the system of equations 


bux? + bux: + bu = x* 
(11) Doxa + Dux + bu = x*i1 
Dogx*i-2 + Duxi-e + bx = x*i2. 


This system of equations may be solved by elimination provided x;, x;1, and 
*X;-2 are distinct and we obtain boo = 1, b21 = bee = 0 for the first three solutions. 
When k > 3 the solution becomes somewhat more difficult but may be carried 
out by eliminating by, between the first two equations to give 


k—2 
borx Xi-1 — by = XXi-1 ) xfxi-2-?. 
p=0 


A similar elimination may be made between another pair and the result combined 
with the above equations to give 


(12) be = xeirti2 YL xfxf_4xj-. 
ptqts—k—3 


The sum is to be taken over all non-negative integral p, g, and s satisfying 
p+q+s =k — 3. All of these results may be inserted in (8) and (10) giving 


(13) f(O) = —xpxi1%i-2 (uo + > x aia!) 


k=4 pt+qts=—k—3 


Up to this point no approximations have been made and no limits have been 
taken. Equation (13) expresses the same relationship contained in (5). We now 
assume that the points x;, x;_1, X;_2 lie in the neighborhood of a root r. Thus if we 
let 641 = Xe — 7, KH HT, G1 KH Xi-1-—-7 and €;_2 = xj-2 — r the magni- 
tudes of the last three quantities are all assumed to be less than some upper 
bound én 


(14) les] < €m, |€i-1] << €m, [| €s-2] < €m- 
We also now make the tentative assumption 
(15) | €c41| < Em» 


and we shall seek to justify this assumption later. If (14) and (15) are inserted 
in (13) and the functions are expanded about r we obtain 


(16) ecnaf"(r) + sdf" (r) + yt’ (”) 
= (€ = €:41) (€-1 - €441) (€-2 a 41) "(r) + 0(e&,)- 
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If r is a single root we see that ¢:,1 = 0(é,) so that we may write 


(17) €i41 = —€s€;_-16:-2 ao + O(€,). 


A solution ¢;,1 to this equation does exist if ¢,, is sufficiently small and will satisfy 
(15). This solution will also satisfy L;(f(xi,:)) = 0, and hence we are justified 
in assumption (15) and hence (17) for at least one of the two x;,: for which 
L;(f(xi+1)) = 0 holds. We now wish to show that (17) holds for that x;,; which 
is actually chosen by the process described in connection with equation (5). If 
both x;41 satisfy (17), the proof need not be given. If, however, one does and one 
does not, we must make some further analysis. It was pointed out that the process 
chooses the point x;,: which is nearer to x;. The point given by equation (17) 
must satisfy 


(18) | Xe41 — x; | = | esa ~ e;| < 2am 


This must therefore also hold for the x;,: which is chosen by the process in (5) 
and hence |¢41| < 3e, for this case. But |€:41| < 3€m, is adequate to give (17) for 
sufficiently small «,, and we may therefore assume (17) for the ¢:,: obtained in 
the process. 

A general limiting formula for the ¢; in the neighborhood of a root may be 
obtained from equation (17). If logarithms are taken on both sides we obtain 


T(r ) 
6f'(r) 


Neglecting the terms 0(e,) we may solve (19) as a difference equation using 
standard techniques and obtain 


(19) log 4.1 = log é+ log 6-1 + log €;-2 + log (- FO) + O(ém ). 





ie (r) 
6f" o) 


Where the constants ¢1, cz, and c; are determined by the starting values and the 
three orders of convergence 1, m2, and m; are roots of the characteristic equation 


(20) log €; = cymy! + cym2i + cum; — 35 log ( - 


(21) m* = m? +m +1. 
The roots are 
m= 1.84, 
M2,m; = —.420 + .606:. 
Since the last two roots have magnitude less than 1 their effect will die out and 


the order of the process is given by m. After these approximations become valid 
we have from (20) 


(22) : €41 = Ke) - 
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where 





«-(-55) - 


In the case of a double root a similar argument exists. Equation (17) is then 
replaced by 


(23) du = FO + od) 


— €;€;-1€;-2 3f"(r) 


and the characteristic equation becomes 


(24) 2m* = m? +m +1. 
It has roots 
m= 1.23, 
Me, Mm; = — .367 + .5201, 


and again the order of convergence is given by m:. We therefore have in the limit 
(25) Conn = Ke} 8 


with 





r-(-§8)" 


Convergence of the Generalized Process. One might imagine a generalized 
process in which an a degree Lagrange interpolation polynomial L;(x, a) is used 
rather than the quadratic of equation (2). This presumes that some new method 
for obtaining the nearest root of this polynomial is to be used. Since the direct 
method corresponding to equation (5) would no longer be practical, one would 
probably use some iterative method for solution of the equation L;(x;4:, a) = 0. 

We now wish to investigate the convergence rate for such a process. 

A general set of equations corresponding to equations (11) may be formed. 
They are 


(26) > baxtai = xt_,, s =0,1,---,a. 
i=0 


When k > athe quantity b., may be obtained by elimination. (This direct method 
for obtaining b., was pointed out to the author by Mr. W. Scott Bartky.) Let us 
eliminate b.-1,, between the first equation and each succeeding equation giving 


(27) XS djs (x FP! — xf) — dan (xs — xe.) = xeei.(xf* — f=), 
1=0 


s=1,2,---,a. 











214 SOLVING ALGEBRAIC EQUATIONS BY AN AUTOMATIC COMPUTER 


Each equation may be divided by (x; — x;_.) giving 


a—2 a—j—2 ; k-—2 
(28) yA DieX Xie > i — bar = XKi-s > xix? s = 1,2, ---,a. 
i=0 1=0 1=0 
We next eliminate b,_2, ---, bo, in a similar manner until the result 


(29) bar = (—1)*x ii +++ Xia D ehh, - ++ aft, 
potpit---+p.=k- (+1) 


is obtained. In this expression the summation is made over all terms for which 


the exponents Po, pi, ---, Pa are non-negative integers and }> p; = k — (a + 1). 
i=0 
We also see directly from (26) that bu. = 0 if k < a, except that dao = 1. 
We may therefore obtain a generalization of equation (17) 


fiom (r) 


(30) €ip1 = (—1)* Hees --- “CE DYO 


+ 0(¢**) 


and obtain for the characteristic equation corresponding to (21) 
(31) met = m+ mei+.-.- +1 


This equation has one root m, which lies between 1 and 2 on the real axis and 
which approaches 2 with increasing a. The remaining roots lie within the unit 
circle and therefore represent perturbations which die out. The order of con- 
vergence of the process to single roots is therefore given by m. Since this can never 
reach 2 we conclude that there is little to be gained in speed of convergence by 
letting a exceed 2. 

One should not ignore the possibility of letting a = 1. In this case the formula 
corresponding to (5) is greatly simplified since a linear equation rather than a 
quadratic now must be solved. This choice, however, suffers from a disadvantage 
if all the coefficients of the original equation are real. If one starts from a real point 
xo then all successive iterative results x; will also be real and hence only real roots 
will be found. 

Tests of the Method. The process with a = 2 as outlined in the preceding 
sections was altered slightly in practice. Whenever the new value of the function 
f (41) is calculated the quantity | f(x,1)|/|f(x.)| is formed. If this latter quan- 
tity exceeds 10 the quantity A: is halved and Aj41%i41, and f(x;,1) are recomputed 
accordingly. With this revision the process has produced convergence in all the 
cases tested. 

Another alteration was made to handle the case in which the denominator 
of (5) is zero. This occurs whenever f(x;) = f(xi-1) = f(xi-2) and in such cases 
the arbitrary value \;,: = 1 is chosen since (5) may no longer be used. 

The process (a = 2) was tested for equations of varying degree. Fourteen 
equations were solved starting with degree 10 and progressing in steps of ten 
through degree 140. Each equation was formed by choosing random points as 








—~ ~~ FP 


war we? 


~ VS 
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roots within the square having vertices + 1 + i. Polynomials were then formed 
from these roots. The solutions to these polynomials were then compared with the 
original random numbers which were used to generate the polynomial. Results 
are summarized in the table. 


Time Taken by 
Accuracy of Last I:z1ac for Com- 


Degree of Accuracy of Least Root to be plete Solution 
Equation Accurate Root Found in Minutes 
10 107 10° 1 
20 10-* 10-* 2 
30 10-5 10-* 5 
40 10- 10° 6 
50 10-* 10-5 10 
60 10-5 1077 12 
70 10- 10-5 17 
80 10- 10-5 20 
90 10-1 10-5 20 
100 — 10-* 33 
110 — — 42 
120 — — 43 
130 — 10-* 48 
140 — — 60 


Dashes in the table indicate that some roots were too inaccurate to be identi- 
fied. In all equations some roots appeared which were correct to 10-* or better. 
The solutions to the 100th degree equation, which had some unidentifiable roots, 
were used to generate a polynomial. All coefficients of this polynomial agreed with 
the coefficients of the original polynomial to at least 6 decimal places. This result 
indicates that the obtaining of accurate values of the roots of the equations of 
higher degree was precluded by the limited accuracy of the coefficients and 
independent of the method of solution. 

The equation x"* — 1 = 0 was solved as an example of a special type of 
equation whose solution could be easily checked. The maximum error occurring 
in any root was of order 10-7 and the time for solution was 70 minutes. 

No equation whose solution has been attempted has failed to yield convergence 
although as indicated in the table, the solutions of equations of large degree may 
be greatly in error. We conclude that convergence in the large does occur in most 
practical cases in spite of the fact that convergence has only been proved for 
single and double roots when the process has brought one to the neighborhood 
of a root. 


Davin E. MULLER 
University of Illinois 
Urbana, Illinois 
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Everett’s Formula for Bivariate Interpolation and 
Throwback of Fourth Differences 


It is well known [1], [2], [6], [7], [8], [9], that it is possible in the case of 
univariate tables for use with Everett’s formula, to eliminate columns of higher 
order differences with practically no loss of accuracy by modification of one or 
more lower order differences through a process known as throwback. That the 
same thing is possible with bivariate tables (and, presumably, with other multi- 
variate tables) seems not to have been recorded in print. 

Everett’s formula for bivariate interpolation, as far as fourth order differences, 
can be written as follows, using symbolism similar to that of [2] (see also [3], 
p. 8-9 and [4], p. viii—xiv). 

(1) Upp = u(xo + ph, yo + pk) 

gq'to,0 + gp'to,1 + pq'tr,0 + pp'ts,1 

— E2(q){q’'5."uo,0 + p’b°uo,1} — E2(p){q/5.2t1,0 + p'5.2t1, 1} 
— E2(q’){q6y?u0,0 + pb,7u1,0} — E2(p’){g5,?uo,1 + 8,721, 1} 
+ E4(q) {q'b24t0,0 + p’bs*to,1} + Ea(p) {q/52°u1,0 + p'5.°u1,1} 
+ E4(q'){Q5y*uo,0 + pby4te1,0} + Ea(p’) {q6y4uo,1 + p5,4u1, 1} 
+ E2(q)E2(q’)6.76,?u0,0 + E2(q)E2(p’)6.76,7uo, 1 

+ Es(p)E2(q’)676,7u1,0 + E2(p)E2(p’)6.76,7u;, 1 


where 
pb+q=1 and ~’+q=1 
and where 
E3(x) = — “ 1) 
and 


ra = ("1") 


We write the above as follows: 


(2) Upp = gq'Uo,0 + gp'Uo,1 + pg'tr0 + pp'ur1 

— Ex(q)q'[6.?uo,0 — cbz*uo,o — d6,75,7u0, 0] 

— E2(q’)q[6,?u0,0 — c’d,4u0,9 — d’5,5,?u0, 0] 

— Ex(q)p’[6.?u0,1 — cds4ue,1 — d6,°6,2u0, 1] 

— Ex(p’)q[6,*u0,1 — c'6,4u,1 — d’5,%,2uo, 1] 

— Es(p)q’(6.?u1,0 — cbs*t1,9 — d5.78,2u1, 0] 

— Ex(q’)p[6,7u1,0 — cby'u1,9 — d’5,76,2u1, 0] 

— Ex(p)p'[82u1,1 — cbatus,1 — 8.28201, 1] 

— Ex(p’)p[67u11 — c’84u1,1 — d'5,76,2u1,1] + R, 











EVERETT’S FORMULA FOR BIVARIATE INTERPOLATION 217 


where c, c’, d, and d’ are constants to be selected so that | R| is “small” for all 
p, p’ on ranges 0 < » < 1,0 < p’ < 1, and where 
(3) R = (Eu) — cE2(q) lq’b4uo,0 + [Ea(q’) — c'E2(q’) ]95,‘uo,0 
+ E2(q)[E2(q’) — dq’ — d’q76,7u0,0 
+ (Es(q) — cE2(q) 1p’6.4uo,1 + [E«(6’) — c’E2(p’) ]96,‘uo, 1 
+ E2(q)(E2(p’) — dp’ — d’q\76,7uo,1 
+ [E«(p) a. cE2(p) ]q’6.*u1,0 + [E.(q’) “? c’E2(q’) \pi,‘u1,0 
+ E:(p)(E2(q’) — dq’ — d’pi76,7uro 
+ [Ea(p) — cE2(p) ]p’b.4u1,1 + [Ea(p’) — c'E2(p’) lpi,'ur,1 
+ Ex(p)LE2(p’) — dp’ — d’p \i76,7uy 1. 
In equation (2) above, it is implied that the entire term in (1) involving 6,“u, ; 
is to be thrown back onto 6,°u; ;, and that the entire term in (1) involving 6,‘u; ; is 
to be thrown back onto 6,7; ;, while the single term in (1) involving 676,’u; ; may 
be thrown back onto both 6,7u; ; and 5,2u; ;. Let us define R, as that part of R 


involving 6,‘u; ;, R, as that part of R involving 6,‘u; ;, and R,, as that part involv- 
ing 675,7;, rn Thus 


(4) R=R.+R, + Ry. 


Since 5,4; ;, 5,/u;,; and 6,°6,?u; ; are essentially independent of each other, it 
seems natural to select.c and c’ so that |R.| and | R,| are small, and take care of 
R,, later. This has the further advantage that it will take care of the special case 
where 6,76,2u is zero or nearly so. We do this by the method of [1] and [2]. In 


[2] it is shown (in effect) that if c = co = (3 + ¥2)/24 ~ 0.1839256, then the 
absolute value of 


(5) R' = [E,(p) — cE2(p) Jiu: + [Es(q) — cE2(q) ]b*uo 


is less than L/2000 provided |é*u| < L and |5'u| < 0.07L. In equation (3) we 
have four groups of terms, two each in R, and R,, which resemble R’. For example, 
R, can be written as 


(6) Re = q {LEs(q) — cE2(g) J6.*uo,0 + [Ea(p) — cEa(p) b.‘u1, 0} 
+ p'{LEs@ — cE2(q) Ws'uo,1 + [Es(p) — cE2(p) 6.411}. 


The terms in braces are precisely analogous to the univariate expression R’. 
Thus, if we select c = Co, it follows immediately (since g’ + p’ = 1) that 


(7) |Rz| < L/2000 provided |éAu|<L and |6,5«| < 0.07L. 
If we also select c’ = ¢o, we get, in a similar fashion 


(8) | R,|< M/2000 provided |8,‘u|< M and |,5u| < 0.07M. 


Now R,, consists of four terms, of which a typical one is 


(9) $(g, 9')6.76,7u0, 0, 
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where 
¢(g,9') = E2(g)LE2(q’) — dq’ — d’q]. 


The natural method of attack on the problem of minimizing the maximum 
absolute contribution of R., is to proceed in a manner similar to that of [2], viz., 
to write 5,76,2u,1 = 6,76,?u;, 0 + 6.76,'u;,, 7 = 0, 1. After these substitutions we 
can write R,, as follows: 


(10) Rey = Gi,76,7uo, 0 + H5,76,7u,, 9 + Ab,76,5u0,4 + Bb 76,u1,3, 


where 
A = E;(q)[E2(p’) — dp’ — d’q], 
B = E,(p)[Ex(p’) — dp’ — a’p], 
G = E,(q)LE2(q’) + Ex(p’) — d — 2d’q], 
H = E,(p)(E2(q’) + E2(p’) — d — 2d’p]. 


If we replace 6,76,?u1, 9 by 6.75,?u0,0 + 5,°5,2u;,0 we can write above equation as 
(11) Ray = G*5,%,2t40,0 + H.%5,25,0 + Ab,%5,%u0,4 + Bd.%,2u1,4 
where 


Gt =G+H = [E.(7) + £2(¢’) — d][E:@ + E2(6)] 
— 2d’'[gE2(q) + pE2(p)]. 


We assume, for the moment, that the three terms in (11) involving fifth order 
mixed differences are “small’’ in comparison with the fourth order term; we 
examine the problem of selecting d and d’ so as to minimize the maximum abso- 
lute deviation from zero of the fourth order term in the region D:0 < p <1, 
O<?p' <1. 

Note that G*(, ’) is symmetric with respect to p = 1/2 and with respect to 
p’ = 1/2. Therefore, we may restrict attention to the region D*: 0 < p < 1/2, 
0<?p' < 1/2. 

Evidently we may assume both d > O and d’ > 0, since any positive extremum 
of G* for d < 0 could be made smaller by taking d > 0; the same is true of any 
positive extremum for d’ < 0. 


We have 
aG*/dap = (1 — 2p){LE2(p’) + Ex(q’) — d]/2 — 2d’[2p* — 2p + 1]/3} 
= — 4d’(1 — 2p)(b — pi) (> — D2)/3 
where 
pi = 1/2 — Vg — 1/4, p2 = 1/2 + vg — 1/4, 
g = [Ex(¢’) + Ex(o) - a} 
Also, 


aG*/ap’ = (1 — 2p’)[Ea(p) + E2(q) 1/2. 
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The extrema of G* either lie on the boundary of D* or satisfy the conditions 
8G*/dp = aG*/dp’ = 0. The latter conditions imply that 


(a) p=1/2, pi or pe 


and 
(b) p’ =1/2 or p=0. 


Since G*(0, p’) = 0, the only non-zero extrema of G* occur on one of the other 
three boundaries of D*. 

Along p’ = 0, g = — 3d/8d’ < 0; hence both p; and p are imaginary here, 
and therefore 0G*/dp < 0 on 0 < p < 1/2, = 0 at p = 1/2. Thus, G*(1/2, 0) 
= — (d + d’)/8 < 0 is the minimum value of G* on boundary p’ = 0 of D*. 

Obviously, 8G*/dp’ > 0 on the interior of D* and along p = 1/2, and 
8G*/dp’ = 0 in D* only along p = 0 and p’ = 1/2. Hence, we can infer that 
G*(1/2, 0) is the absolute minimum of G* in D* and the only possible maximum 
is along p’ = 1/2. We therefore next investigate G* along p’ = 1/2. 

Along p’ = 1/2, g = 3(1/8 — d)/8d’ and g — 1/4 = (3 — 24d — 16d’) /64d’. 
It can be shown, by examination of 8G*/dp, that 


a. if 3 — 24d — 16d’ < 0, G*(1/2,1/2) is a (negative) minimum and the 
absolute maximum of G* on D* is therefore 0. 
b. if 3 — 24d — 16d’ > 0, 
and (1) 3 — 24d — 32d’ < 0, G*(1/2,1/2) is a maximum, absolute if 
positive, 
and (2) 3 — 24d — 32d’ > 0, G*(1/2,1/2) > 0 and the absolute maxi- 
mum of G* on D*. 


In all cases, G*(1/2, 1/2) = 1/64 — (d + d’)/8. 

Since we wish to minimize the maximum absolute deviation of G* from 0, since 
G*(1/2, 1/2) — G*(1/2, 0) = 1/64, and since both G*(1/2, 1/2) and G*(1/2, 0) 
are monotonic increasing functions of (d + d’), we can evidently solve this 
minimax problem by making G*(1/2, 1/2) > 0 and =|G*(1/2, 0)|. From this it 
follows that d + d’ = 1/16 and maximum absolute deviation of G* from 0 is 
1/128. Note that the condition d + d’ = 1/16 is consistent with case (b) (2) above. 

Having determined that we should select d + d’ = 1/16, whence |G*| < 1/128, 
we now consider the selection of d (and of d’ = 1/16 — d) and its effect on fifth 
order terms in R,,. 

An investigation of the coefficients A, B, and H of the fifth order terms in 
R., by methods similar to the investigation of G* shows that (when d’ = 1/16 — d) 


(1) Maximum |H| is a monotonic decreasing function of d, 0 < d < 1/16. 

(2) Maximum |A| and maximum |B| are monotonic increasing functions of 
d on above range, and have same values, although attained for different 
values of (p, p’). 


Thus, if |5.°6,24| << N, we can say that 


(12) |Rey|< N/128 + max |H| max |6,°6,2u|-+ 2 max |A| max |6,76,*u|. 
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Unless we make some assumption about the relative magnitudes of 6,°5,24 and 
6,75,*u, there is no rational basis for selection of the value of d other than a selection 
based on the fact that max|H| is monotonic increasing and max|A| is monotonic 
decreasing. For example, if we assume that max|é,*s,2u| = 2 max|676,5u| = P, 
then we would be confronted with a problem of minimizing P[max|H|+ max|A|]. 
If this were the case, it turns out we should choose d close to 1/32 (midpoint of 
interval 0 < d < 1/16) to minimize [max|H|+ max|A|]. Since this choice at 
least has the virtue of symmetry (if d = 1/32, then d’ = 1/32), we make it. 
With the above choice, we find that 


max|H| < 0.00451 
and 
max|A| < 0.00323. 


However, maxima above are attained for different (p, p’) from the pairs for 
which |G*| has its maximum, viz., (1/2, 1/2) and (1/2, 0). The maximum values 
of |H| and |A| for these pairs of arguments are respectively 1/256 ~ .00391 
and 1/512 = .00195. 

In order to determine a reasonable set of restrictions on |6,°4,?u| and |5,°6,42|, 
we remember that we already have assumed the fifth order terms ‘‘small’’ in 
comparison with the fourth order term whose maximum absolute value is V/128. 
If, e.g., N < 60 so that V/128 < .46875, and this term constitutes most of | R.,|, 
then in (12) we may allow max|676,24| = max|é,*6,'u| < 4 and have |R.,| < 1/2. 
Thus, if |5,°6,20| << N and |6,%6,2u|, |68,5| << N/15, then |R.,| < N/120. 

Obviously, other sets of restrictions are possible; e.g., we might allow the two 
fifth order differences different maxima, or we might allow them to be larger 
fractions of total | R.,|. A similar result would be: if |6,6u|< N and |6,°6,?u|, 
|5.°5,2u| < .28N, then | R.,| < N/100. 

However, it seems natural to select a restriction set corresponding roughly to 
that chosen for |6.4~| and |.5w| in [2], viz., L and 0.07Z (thus, 1/15 ~ 0.07). 

In summary, then, we can say that with the choice c = c’ = ¢) and d = d’ 
= 1/32 for the “‘throwback’”’ constants, 


(13) |R| < L/2000 + M/2000 + N/120 


provided 
\8.4u| < L, |8.5u| < 0.07L, 
|5,4u| < M, |5,5«| < 0.07M, 
|8,76,2u|< N, |6,%5,3u| < N/15, 
and 
|5.°6,2u| < N/15. 


In some instances, it may be desirable to have available the general result 
corresponding to (13), which does not necessarily assume that the fifth order 
differences are ‘‘small’’ with respect to the fourth order terms. 

In order to see how we get such a result, consider, as in [2], that we have 
written equation (5) as follows: 


(5’) R’ = Giuy + [Ea(p) — cE2(p) ]6*uy 
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where 


G = p(p’ — 1)(p — 2)/24 + cp(p — 1)/2. 


With the choice c = co (which minimizes maximum |G]|), we find that 
|G| < .000447 and | E4(p) — cE2(p)| < .000792. Consequently, if |5*wo| = 1000 
and |6°u,| < 70, we would have | R’| < .503, which is approximately 1/2. Indeed, 
if we make use of the fact that max|G| and max| E4(p) — cE2(p)| occur for differ- 
ent values of p, and that at p = .88265, |G| = max, whereas | E,(p) — cE2(p)| 
= .000744 instead of .000792, we find that, probably, | R’| < .499. From this 
arbitrary selection of max |é‘u| = 1000 comes, then, the result stated in connection 
with equation (5). Evidently, the corresponding general (and conservative) result 
is that if |8*w| << Land |é*u| < KL, then 


(5’’) | R’| < L(.000447 + .000792K). 
Obviously, this implies that 

(7) |R.| < L(.000447 + .000792K,) if |é4u|<L and |é6,5u|< Kil. 
Similarly, we have 

(8’) |R.| < M(.000447 + .000792K:.) if |éu|< M and |6,5u|< KM. 


(In the above, the only restrictions on the K; are that they be >0.) 
From the discussion in connection with equation (12) above, we see that 
we have 


(12’) |Rey| < N[1/128 + .00451K 32 + .00646K:3] 
if 
| 526,74 | < N, | 5.°5,21 | < KN, 


and 
|6.°8,5u| < KesN, 


where, again, the only condition on the K;,; is that they be positive. 
Consequently, we have the following general result to replace (13): 


(13’) | R| < L(.000447 + .000792K:) + M(.000447 + .000792K:) 
+ N(1/128 + .00451Ks2 + .00646K:s) 


under the hypotheses governing (7’), (8’), (12’). 

The advantages of the use of modified second differences as in equation (2) 
are obvious—a published table need contain no fourth order differences and 
formula (2) requires only 24 multiplications to evaluate 12 terms as opposed to 
the 48 multiplications required to evaluate 24 terms in formula (1). 

It should be noted that the fraction 1/120 is disappointingly large in com- 
parison with 1/2000. However, there really seems to be nothing one can do about 
this situation. Even so, we have “‘gained”’ in some sense, since the condition for 
negligibility of terms involving fourth order (pure and mixed) differences is that 


(14) L/42 + M/42 + N/64 < 1/2 (units of last decimal place tabulated) 
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whereas the condition for ‘“‘throwback ability’’ of the same differences is 
(15) L/2000 + M/2000 + N/120 < 1/2, 


in the presence of auxiliary conditions in (13). 

As an example of the use of throwback, we may consider its application to the 
tables in [5]. Cursory examination of [5] indicates that in a large proportion of 
the table, inequality (15) (or its more general version corresponding to (13’)) is 
satisfied, so that modified second differences could be made to work. As a specific 
illustration, we take the numerical example on p. xxiv—xxv, viz., computation of 
I(4.025, 7.05). We use formula (2) as modified by introduction of the throwback 
constants ¢, c’, d, d’ selected as above, viz., 


(2") Upp ~ Qq'to,0 + gp'uo1 + Pg'ur0 + Pp'ui1 
— E:(q)q’M.?uo,o0 — E2(q)p’MPuo.1 — E2(p)q’M.2u1,0 — E2(p)p’MPu11 
— Es(q')qM,7uo,0 — E2(p’)qgMuo1 — Es(q')pM,7u1,0 — E2(p’)p>M,u11 


where M,?u, ; is the modified second difference of u, with respect to the variable z, 
evaluated at “point” i, 7. 

At the point in question, (4.0, 7.0) ~ uo, o, I(4.1, 7.0) ~ u1, 0, I(4.0, 7.2) ~ uo, 1, 
and [(4.1, 7.2) ~ 1,1. Thus, the full array of differences for use in the unmodified 
formula would be: 








82 by? 
54 by! 
5,6,? 
(0, 0) (0, 1) 7 
—17520 —1747 —17985 —1745 
63 0 80 0 
16 17 
(1, 0) (1, 1) ' 
—16319 —1643 —16802 —1649 
30 0 43 0 
10 13 














This is replaced by the following set of modified differences : 








M2 M, 
(0, 0) (0, 1) 

—17532 —1748 —18000 —1746 
(1, 0) (1, 1) 

—16325 —1643 —16810 01649 
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Here we have p = 0.25, p’ = 0.25. Using (2’), together with a table of E;(x), 
as in [2], we get J(4.025, 7.05) ~ .878, 5410, which agrees with result obtained 
by Pearson. 

In principle, throwback of other kinds in the bivariate Everett formula is 
possible. There is, also, no a priori reason why it could not be accomplished in 
the trivariate Everett formula, which is stated in [4], page x. 


Tuomas H. SouTHARD 
University of California 
Los aie California 
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On the Treatment of Monte Carlo Methods in Text Books 


Though Monte Carlo methods have not yet reached the text book stage, they 
are gradually being afforded a brief mention in new books on numerical analysis. 
This is probably stimulated by a desire for completeness but the brevity of treat- 
ment often results in a very narrow picture of Monte Carlo. For any given problem 
there may be several methods of solution by Monte Carlo and there are general 
principles which should guide one in choosing between these. A description of a 
single method, leaving the reader with the impression that this is the sole method, 
is therefore a mistake. Care must also be taken to see that, though a method may 
be illustrated on a simple example for which Monte Carlo would not normally 
be used, it is easily extensible to problems for which Monte Carlo would be needed 
and provides a reasonably practical method of solving these. 

In two otherwise excellent books by Householder [2] and Kopal [4] a tech- 
nique is given for evaluating integrals by Monte Carlo. It is the only example of 
Monte Carlo in either book and is rather an unfortunate one; for it is not easily 
applicable in practice to multiple integrals though it is largely for these that 
Monte Carlo can be useful; in addition, it is always less efficient, and often much 
less so, than another well-known simple technique (defined as crude Monte Carlo 
by Hammersley and Morton [1]). Their comparative efficiency is also an ad- 
mirable illustration of a general precept to be used when choosing even the 
simplest Monte Carlo methods. 
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Consider the estimation of the integral 


la f " $(x)dx 


where f(x) is a prescribed function whose minimum is made zero by a change of 
origin and whose maximum is a positive constant M. In the method described in 
[2] and [4], we draw a rectangle of height M and unit breadth about the curve 
y = f(x) and choose m points (&;, Mn), i = 1, 2, ---, m, at random in this rec- 
tangle; &; and »; are independent random numbers uniformly distributed in the 
interval (0, 1). We count a score #;; = ¢,/£;, 9;) for each point, t:; being M if the 
point lies between y = f(x) and the x-axis and zero otherwise. Then the mean 
score t; = (1/m) >-?_; #1; is an unbiased estimator of J and i,/M is binomially 
distributed with parameter J/M, i.e., 


e(i:) = I and vari; = I(M — I)/n. 


The crude Monte Carlo technique, first described by Kahn [3] in 1949, 
consists of using #2 = f(£) as an estimator of J. Thus, for the 2 random numbers 


€:, we have a mean score é, = (1/n) } f(é,) for which clearly 


e(iz) = I and vari, = (1/n) {f f(x)? dx — r| = (S — I*)/n. 


It is clear that, when dealing with multiple integrals, the bounds of the 
integrand which it is necessary to know for the first method may be very difficult 
to obtain even when finite; and that any overestimaté of them can lead to a 
considerable increase in var t,;. Moreover, since f(x) > 0, S*< IM so that 
var te < var t;, the equality occurring only when f(x) is a copataitt the impor- 
tance of this difference is enhanced by the fact that it is the ratio of the variances 
which measures the relative efficiency of the two methods. For example, when 
f(x) =x, IM = 1/2and S = 1/3 but var t; = 1/4 = 3 var te, so that the second 
method is three times as efficient as the first. Finally, it should be noted that 


1 
t= jf ti(é, »)dn, i.e., M times the probability that, for given £, the point 
0 


(é, Mn) will lie below y = f(x): thus the second estimator is simply related to the 
first, being its average over 7, and its higher efficiency is due to following the 
principle that random processes should be replaced by analysis wherever possible, 
as pointed out in [1]. 

K. W. Morton 
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On the Simultaneous Numerical Inversion of a Matrix 
and All its Leading Submatrices 


This note considers the case in which a square matrix A may be factorized 
(1) A = X'DY, 


where X and Y are square matrices whose diagonal elements (where i = 7) are 
unity, and whose subdiagonal elements (where i > 7) are zero, whilst D is a 
diagonal matrix, that is, one whose off-diagonal elements (where i + 7) are zero 
[i]. Where this factorization exists, it is unique [1]. Where it does not exist, 
there exists another factorization of the form (1) in which the matrices X and Y 
differ from the foregoing specification by having their rows in some way. inter- 
changed [2]. 

Let the matrices appearing in relation (1) all be partitioned after the rth row 
and column, the resultant leading submatrices of order r by r being denoted by 
A,, X,', D,, Y,, where r = 1, 2, ---,m — 1. Then it is easily seen that 


a => YD. Y,. 
If D, is non-singular, so is A,, and hence 
(2) At = Y-°D“"X,'-. 


This holds also for r = m, when suffixes may be omitted. It is, in fact, equiva- 
lent to taking r = m and replacing the elements in the last  — r columns and 
rows of all matrices by zero. This suggests that if a typical element, say },;, of 
A-— = B is obtained as the scalar product of the ith row of Y"D~— and the jth 
row of X-', the corresponding element b,j.) (where 7, 7 < r) of A; = B,, say, is 
obtained as the rth partial sum of the m products of the row-by-row multiplication. 
Thus, with a desk calculating machine, on which products can be accumulated, 
the inverses B;, Bz, ---, of A1, As, --- are obtained simply as the partial sums 
of the row-by-row multiplications needed anyhow for the computation of the 
inverse B of A. 

As the factorization is unique any process yielding X’, D, and Y (or Y“D“ 
and X’-", etc.) may be used. The following arithmetic check is suggested. Relation 
(1) gives 

X’14 = DY. 


Comparison of diagonal elements on either side shows that the scalar product 
of the rth column of X,— into the rth column of A, is equal to d,, the rth diagonal 
element of D,. This relation may be used to check each column of X~ fresh after 
computation. Similarly it is found that the scalar product of the rth column of 
Y,—' into the rth row of A, is also equal to d,. 

H. NAGLER 


25 Helenslea Avenue 
London N.W.11 
England 
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Moti. A. ig mee “Rounding-off errors in matrix processes,” Quart. Jour. Mech. and Appl. 
v. Pp 

2. R.A. Frazer, W. jy. Duncan, & A. R. Cotrar, Elementary gg and Some Applications 
to Dynamics and Differential Equations, Cambridge Univ. Press, 1938, p. 99 


A Note on Large Linear Systems 


Towards the end of 1955 the author of the present note was charged by 
Svenska Aeroplan Aktiebolaget (SAAB Aircraft Company), with setting up a 
program for the automatic computer BESK for the solution of large symmetric 
systems of linear algebraic equations. The reason was that the solution of 4 
systems with 214 unknowns was required. The systems had been set up in con- 
junction with certain calculations for airplane structural analysis. Approximately 
40% of the coefficients in the matrices of these systems differed from zero. This 
fact was taken into account in setting up the program. Nevertheless, the whole 
system could not be accommodated at one time in the magnetic drum memory of 
8192 words associated with BESK, and a splicing method was developed which 
does not involve excessive loss of time in connection with the input and output of 
intermediate results. The four systems were solved simultaneously by the BESK. 
The method used for solving the equations was Gauss’ method of elimination. 
The calculating time was about 2 hours. If p% of the coefficients of the systems 
differ from zero, and the order of the systems is m, the calculating time will be 
proportional to p*n'. 

By formation of residuals and one iteration a solution was obtained which, 
when inserted into the original systems, proved to satisfy the equations with 
deviations less than 10-*%, i.e., if the computed solution was Z; and the system 
was AX = B, the “error vector’ E = > a;;4; — 6; satisfied |EZ| < 10-*|B|. The 


J 
systems were badly conditioned since, if the matrix was (ax), the ‘condition 
number”’ was: 


|determinant (a,x) | 
24, 214 
I (= Jaa! ) 
kat \i=t 


The matrix is therefore comparable with the m X m matrix 


R(A) = 


= 2-886 





5 —4 1 

—4 6 —4 1 . 
1 —4 6 —4 i , 
: 1 4 6 —4 1 
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studied by Todd, for which R(A,) ~ c-n*-2- for a certain constant c. No direct 
comparison of the results of our experiments with those made by Todd [1] (for 
a 49 X 49 matrix, without iteration) is possible, since his results are expressed in 
terms of the norm of the error matrix. 

The calculating program thus described operates at a fixed binary point and 
uses 40 binary digits. A program operating at a floating binary point and including 
32 binary digits in its digital portion and 8 binary digits in its scalar exponent 
has also been developed. This program works at an approximately 30% slower 
rate and gave, when applied to the above-mentioned systems, a slightly better 
solution. 

Hans RIESEL 
Essingeringen 9 
Stockholm 
Sweden 


1. J. Topp, “‘The condition of certain matrices II,” Arch. d. Math., v. 5, 1954, p. 249-257 
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65[A ].—Rosert L. Causey, Decimal to Octal and Octal to Decimal Conversion 
Tables, U. S. Naval Air Missile Test Center, Point Mugu, California, 1952, 
30 p., 26.5 cm. Deposited in the UMT Fine. 


These are radix tables with three figure groups; when using the decimal to 
octal table the tabular entries are added together using octal arithmetic. The 
decimal to octal table has 11 place octal values of m X 10-* for n = 0(1)999 and 
K = 3, 6, 9. The octal to decimal table has 10 place decimal values of  X 8-* 
for m = 0(1)511 and K = 3, 6, 9, 12, together with an argument column giving 
n X 8-* in octal. The tables were printed by photo-offset from output of an IBM 
tabulator. This machine has the unpleasant characteristic of leaving-more space 
between adjacent columns than rows. The tables were developed for use with 
the RAYDAC. 


J. L. SELFRIDGE 
University of California 
Los Angeles, California 


66[D, E].—Harvarp University, Computation Laboratory, Annals, Tables of 
the Function arc sin z, v. 40, Harvard University Press, Cambridge, Mass., 
1956, xxxviii + 586 p., 27 cm. Price $12.50. 


This table lists values of w = arc sin z, both the function w and the argument 
z being expressed in cartesian form. Symmetries make it unnecessary to consider 
values of z belonging to points outside the first quadrant. If w = u + iv and 
z = x + ty, both uw and v/are tabulated to 6 decimals for various values of x 
and y. The results are arranged in 11 tables, in every one of which the lattice is 
square, that is, the intervals in x and y are equal. The intervals 6 = Ax = Ay in 
Tables I to XI are respectively 0.002, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2, 5. 
In a table of a function of a complex variable, the earlier intervals are strikingly 
small; previous tables of sin z and arc sin z have rarely used intervals finer than 
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0.02. The regions in which the various intervals apply are clearly set out in detail 
on the contents page, and need not be reproduced here. It seems sufficient to 
state that 


(i) the smaliest intervals relate to the regions which keep closest to the 
branch point at z = 1; 

(ii) the region covered by Table I (6 = 0.002) is the rectangle 0.916 < x 
< 1.090, 0 < y < 0.140; 

(iii) the other regions are not single rectangles, but are each composed of two 
or three rectangles, the maximum width of the composite region vary- 
ing in x from 776 to 1398, and in y from 956 to 1058; 

(iv) the whole region covered in the volume is the square 


0<x<475, 0<y< 475, 


minus that part of the square for which both x > 195 and y > 350. 


In general, interpolation is by Taylor’s series, using the first two derivatives ; 
the quantities actually tabulated are not these derivatives, f’(z) and f’’(z), but 
the ‘‘pseudo-differences” 5f’(z) and &f’’(z)/2, of which the real and imaginary 
parts are given. (The reviewer regrets the use of the symbol 6 to denote width of 
interval; he would have preferred to call it almost anything else, for example h.) 
Where use of two derivatives does not suffice, reference is given to the appropriate 
procedure—usually the employment of an earlier (finer mesh) table, but in the 
case of a region of Table I near to z = 1, the employment of an auxiliary function 
tabulated on pages xxxii and xxxiii. 

It is interesting to see what an extensive table of a function of a complex 
variable can be produced in a fairly routine fashion with modern equipment. 
The present table was produced by the Computation Laboratory for the Aero- 
nautical Research Laboratory of the United States Air Force, and its extent was 
presumably fixed by the requirements of that institution. Previous tables of both 
the direct and the inverse sine of complex argument have been mostly the work 
of electrical engineers. As the authors point out, the real and imaginary parts of 
the direct function sin z are given so simply by the expressions sin x cosh y and 
cos x sinh y that tabulation seems unnecessary, though the importance of the 
function has caused a number of tables to be made. In contrast, the computation 
of the inverse function arc sin z from ordinary tables of functions of one real vari- 
able is rather tedious, yet, as the authors rightly state, not much tabulation of 
the inverse function has been undertaken ; useful tables listing a cartesian function 
of a polar argument were computed in the 1920’s, the moving spirit being F. Emde. 
Concerning the references on page xxiv, it may be added that an abridgement of 
Hawelka’s table has been printed by Boll [1], and that the small table given in 
the second edition of Jahnke and Emde was later transferred to Emde’s volume 
of tables of elementary functions [2]. As the present tables list a cartesian func- 
tion of a cartesian argument, they break new ground as regards form, besides 
being of record extent. They are manifestly to be welcomed for their own sake. 
They could also provide material if it should ever be desired to produce a smaller 
and cheaper table, such as would satisfy the needs of many. The reviewer ven- 
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tures this suggestion in case any persons should feel that a little acorn planted by 
Emde has grown into a splendid oak which is too big for their garden. 


A. F. 


. Marcett B Numériques Universelles des Laboratoires et Bureaux d'Etudes, 
Dunod, Paris, 1947p. 132 on (see MTAC, v. 2, 1947, pes | 
E, Tafeln elementarer Funktionen, Leipzig and i 1940, second edition, Leipzig, 
1945, p nS (see Mrac. v. 1, 1945, p. 384 and v. 4, ane 79). There ‘are American reprints of 
betty dhogeunell edition of Jahnke and Emde and the first edition of Emde. 


67[F].—R. J. Porter, “On irregular negative determinants of exponent 9m,” 
MTAC, v. 10, 1956, p. 22-25. 


The determinant D = 5b? — ac of the quadratic form ax* + 2bxy + cy* is 
called irregular if the principal genus is non-cyclic. The number of classes divided 
by the highest period is called the exponent of D. The first 58 values of —D with 
exponents 9n are listed, giving the complete table up to 150,000. Tables are given 
which reveal special properties and a method for squaring a class of forms is out- 
lined. The discriminant —3299 seems to play an exceptional role in various 
respects. 


O. Taussky-Topp 
National Bureau of Standards 


Washington, D. C. 
This review was prepared by O. Taussky-Todd for Mathematical Reviews. 


68[F, K, W ].—Gésta NeEovivs, “Artificial traffic trials using digital computers,” 


Chapter 2. Generating pseudo-random numbers, Ericcson Technics 2, 1955, 
p. 279-291. 


Five methods of generating pseudo-random numbers to radix 16 are described. 
The first method produces the Fibonacci series 1, 1, 2, 3, 5, 8, 13, --- reduced 
modulo 2*° and the others are rather complicated modifications of this method 
which, however, are also convenient for the electronic computer BESK. For a 
description of the mathematical properties of Fibonacci sequences modularly 
reduced, see Duparc, Lekkerkerker, and Peremans [1]. The four modified 
Fibonacci sequences, but not the unmodified sequence, were found to satisfy the 
four tests for randomness described by M. G. Kendall and B. Babington Smith 
[27]. It was in a sense fortunate that the radix used was large (16) instead of 2; 
otherwise the serial test would certainly have led to the rejection of all the 
methods, unless this test had been interpreted in accordance with the reviewer's 
modification [3] of the analysis in [2]. 


I. J. Goop 
25 Scott House 
Princess Elizabeth Way 
Cheltenham, England 
wali H. J. A. Duparc, C. G. LEKKERKERKER, & W. PEREMANS, Reduced sequences of —— 
aw, out fs oabesanton numbers, Math. Centrum, ‘Amsterdam, Report ZW 1953-002 (also Report 


. M. G. Kenpatt & B. Basrncton Suita, Tables of Random Sampling Numbers, Tracts for 
Computers XXIV, Cambri diy Aap Press, 1951. 


. I. J. Goon, * "The or sampling numbers and other tests for randomness " Camb. 
Phil "Soc., Proc., 29, 1953, p. 26-284. 
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69[G, H, X].—Marx Lorxm, “A set of test matrices,” MTAC, v. 9, 1955, 
p. 153-161. 


The finite Hilbert matrix B, = ((¢ + k — 1)), i, k = 1,2, ---, 2 is often 
used for testing numerical methods in matrix computations. In order to have test 
data available for a non-symmetric matrix a study and tabulations are carried 
out for A, = (@#) with ay = 1,k =1,---,m and ag = (¢+k2 — 1) for 
4=2,---,n,k =1,---,m. It is known that the det B, is the reciprocal of a 
positive integer and it is now shown that det A, is (—1)*~ times the reciprocal 
of a positive integer. The elements of A, are integers, which is also known to be 
true for B,—. While all the characteristic roots of B, are positive the tabulations 
indicate that for A,, apart from the dominant characteristic root, all others have 
negative signs. Estimates for the characteristic root and vectors are studied as 
well as the M- and P-condition numbers (see John Todd, “‘The condition of a 
certain matrix,” Cambridge Phil. Soc., Proc., v. 46, 1950, p. 116-118). The paper 
includes tables for = 1 to 10 of det A,, the M- and P-condition numbers, the 
characteristic roots of largest and smallest absolute value and the corresponding 
vectors; Ag“! and Hs = (A¢)7A¢ are also given. Results are given exactly, or to 
about 7S. 


O. TaussKy-Topp 
National Bureau of Standards 


Washington, D. C. 
This review was prepared by O. Taussky-Todd for Mathematical Reviews. 


70[G, K].—F. N. Davip & M. G. KENDALL, “Tables of symmetric functions, 
Part V,”’ Biometrika, v. 42, 1955, p. 223-243. 


These tables complete the services by the authors {1] in which each of the 
4 kinds of fundamental symmetric functions, monomial, one-part, unitary, and 
homogeneous product sums are expressed in terms of each of the others for weights 
up to and including 12. The present tables are for the one-part, unitary pair. The 
table for the one-part, homogeneous product sum pair would be the same except 
for obvious changes in sign. 

AE oA 
. F. N. Davip & M. G. KENDALL, ‘‘Tables of symmetric functions Part I ; Parts II, and III’’; 


Part IV, Biometrika, v. 36, 1949, p. 431-449; v. 38, 1951, p. 435-462; v. 40, 1953, p. 429-446. 
[MTAC, v. 4, p. 146, v. 6, p. 224-225; see also corrigendum, v. 8, p. 188, v. 8, p. 150.] 


71[1, L, P].—Ross E. Graves, ‘Design and Analysis of an Optimum Comb 
Filter Based upon an Information Storage Device,” GERA-20, Goodyear 
Aircraft Corp., Litchfield Park, Arizona, 1953, ii + 65 p.; “Miscellaneous 
Considerations in Comb Filter Design,’” GERA-24, Goodyear Aircraft Corp., 
Litchfield Park, Arizona, 1954, v + 74 p., photostat copies in notebook 
binders, 28 cm. Both copies deposited in the UMT FILe. 


These papers contain, among other things, several devices for computation of 
the coefficients in the (finite) Fourier series for the function T,_1(Z» cos u), where 
Tx(x) is the Chebyshev polynomial of degree K, defined by cos (K arc cos x). 

GERA-20 has two recursion formulas for computing the above coefficients; 
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numerical examples of their use for » = 40 relations between the coefficients and 
the hypergeometric function, plus relations between the coefficients and the 
Jacobi polynomial, are given in the Addenda. 

GERA-?24 includes “An Exact Method for Computation of the Tchebycheff 
Coefficients.”” The coefficients referred to are proportional to the above coefficients 
(in GERA-20), the so-called weights for an n-signal Chebyshev filter. The formula 
given is stated to be well suited for computation. 


Tuomas H. SoOUTHARD 
University of California 
Los Angeles, California 


72[1, L, P].—IsaBELLE ArsHam, Chebyshev Coefficients for Chebyshev Polynomials 
of Orders 12 and 24 under the General Linear Transformation (U), Report 
No. TR-326, Diamond Ordnance Fuze Laboratories, Ordnance Corps, De- 
partment of the Army, Washington, D. C., 1956, 27 p., 26.6 cm. 


Consider the equation 7,(ax + 6) = },.0 bmZm(x), where T,,(x) is the 
Chebyshev polynomial of degree m. 

This paper lists the coefficients b,,(m = 0,1, ---, ) as functions of a and B, 
for the cases m = 12 and 24. 

The special case in which 6 = 0 is related to the problem considered in 
GERA-20 and GERA-24 (see preceding review 71). The },, are needed in such 
situations as antenna design. 


Tuomas H. SouTHARD 
University of California 
Los Angeles, California 


73[K].—NBS Applied Mathematics Series, No. 21, A Guide to Tables of the 
Normal Probability Integral, U. S. Gov. Printing Office, Washington, D. C., 
1952, iv + 16 p., 26.0 cm. Price $0.15. 


This guide is not meant to be either an exhaustive or historical compilation 
but rather to assist a researcher or student to find the normal probability table 
he needs in standard texts and other sources likely to be available to him. The 
Guide is in two main portions. In the first, part I lists tables of segments of area 
under the normal curve corresponding to the abscissa x of five kinds: from —x 
to x, from 0 to x, exterior of the interval, (—x, x), from x to ©, and from — 
to x. Part II contains list of tables inverse to those in part I, i.e., for given areas 
of any of the five kinds, the corresponding x is tabulated. Each set of listings is 
headed by a sketch of the area in question, the formula for it, and a key showing 
a few illustrative tabled values. The actual listings contain for each table: (i) the 
number of decimals or significant figures, (ii) intervals and ranges of the argu- 
ment, (iii) indication of the differences, derivatives, or proportional parts given, 
and (iv) authorships, date, and page or table number. Part III gives descriptions 
of appropriate interpolation methods and Part IV gives relationships between 
some types of normal tables and other tabulated functions. The second main 
portion is the bibliography of tables arranged by author and year for reference 
from column (iv) of the listings in parts I and II. The arrangement and notation 
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follows that in Chapter 15 of A. Fletcher, J. C. P. Miller, & L. Rosenhead, An 
Index of Mathematical Tables, McGraw-Hill Book Co., New York, 1946 from which 
much of the content is taken. This Guide will be very useful and should adequately 
meet the ordinary needs of the researcher or student. 


C. c. C. 


74(K ].—Biometrika, Editors, Editorial, ‘‘The normal probability function: 
Tables of certain area-ordinate ratios and of their reciprocals,’ Biometrika, 
v. 42, 1955, p. 217-222. 


In distribution problems such as that of correlation between two continuous 
variables when one is dichotomized or in the case of percentile points one fre- 
quently is interested in the ratio between the area under the normal distribution 
to the left or right of a given deviate and the corresponding ordinate. The tables 
given in this editorial give values of the four ratios P/Z, Q/Z, Z/P and Z/Q for 
X = 0(.01)3. Here 


1 
Z= Z(X) = Fale 
T 


z 


1 oo 
P= P(x) =f" Terran, Q=1-P = f poeta 
a Te 


The ratios P/Z, Q/Z and Z/Q are tabled to 5D and those for Z/P to 5S. 

These tables form useful companions for Table II, Tables for Statisticians and 
Biometricians, Part II, which give the same four ratios for an argument of P 
instead of X. Both tables will be included in the second volume of Biometrika 
Tables for Statisticians. ‘ 


C. F. Kossack 
Purdue University 


Lafayette, Ind. 


75(K ].—D. Tr1cHroEw, “Numerical Analysis Research unpublished statistical 
tables,” Amer. Stat. Assn., Jn., v. 50, 1955, p. 550-556: 


This note describes tables connected with probability and statistics which 
were computed at what is now Numerical Analysis Research, University of 
California at Los Angeles. It is likely that many of these tables will not be pub- 
lished, and most of them exist on punched cards. The tables are listed and indi- 
vidually summarized under 5 categories: 


I. Tables associated with the normal distribution. 
Most of these tables are of functions useful in connection with computing 
moments of order statistics. 
II. Tables associated with the gamma distribution. 
These are tables used in the computation of probability points of x?. 
III. Tables associated with the #-distribution. 
These are tables used in computation of probability points of ¢. 
IV. Tables for selecting samples from certain distributions 
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V. Miscellaneous tables. 

1. Multinomial coefficients 

2. Coefficients for curve fittings by Chebyshev polynomials 

3. Tables for Probit Analysis with Poisson error models. 

M. G. NATRELLA 

Statistical Engineering Laboratory 
National Bureau of Standards 
Washington 25, D. C. 


76[K ].—A. A. Anis, ‘The variance of the maximum of partial sums of a finite 
number of independent normal variates,’ Biometrika, v. 42, 1955, p. 96-101. 
Let %1,%2,°++,%, be independent standard normal variates, let S; = x; 

+ x2+ --- +x, and let U, = max {5S:, S2, --- S,}. It was shown by Anis and 

Lloyd [1] that E(U,) = *z} 1-4/2 which is known to be [2Vn — 1 + ¢(4) 

+ 0(1/Vn) 1/2. In the present paper it is shown that E(U,?) = (m + 1)/2 


+ 522? ¥ jaifs(r — s + 1)}*/2s. There is a table of E(U,), E(U,2) and «(U,) 
to 4D for m = 3(1)25. 


Anis approximates the double sum by a double integral in order to write 
E(U,2) ~n -— Vn(2 + V2)/" to terms of order Vn. This formula, together 
with that for E(U,) mentioned above, implies o(U,) ~ V(x — 2)/xVn — 
[2 + V2 + 2¢(4)]/2Vx(e — 2) which does not agree well with the tabled values. 
Using another method, the reviewer finds V(m — 1)/1+ Vin — 2)/2 + --- 
+ Vi/n —1~an/2+ (a) Vn with the aid of which the double sum may be 
evaluated. We get E(U,”) ~n + 2¢ (4) Vn/ «x and hence find constant term 0 in 
the expression for o(U,). This is in agreement with the tabled values, which 
exceed 0.60281Vn by about 0.478/Vn. 





J. L. Hopcgs, Jr. 
Univ. of California 


Berkeley, California 


1. A. A. Ants & E. H. Lioyp, “On the range of partial sums of a finite number of independent 
normal variates”, Biometrika, v. 40, 1953, p. 35-42. 


77[K ].—Josern Berkson, “Estimate of the integrated normal curve by mini- 
mum normit chi-square with particular reference to bio-assay,”’ Amer. Stat. 
Assn., Jn., v. 50, 1955, p. 529-549. 


This article is the latest of a series of discussions of the relative merits of 
various methods of estimating the parameters (a and @) of the straight line 
transform of the integrated normal curve, used to a great extent in bio-assay. 
The work “‘normit” is coined to represent a normal deviate, i.e., 


vy; = normit P; = a + B;, 


1 at+fzi 
Pp = fra, 
V2" J—0 


and x; is a function of the dosage producing P; kills. 


where 
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The author favors the minimum normit x? method over that of maximum 
likelihood to estimate a and 8 from a set of x; and p; (estimate of P;). He sub- 
stantiates this preference with the results of a sampling experiment consisting of 
600 samples for each of four sets of P’s, each set consisting of P; (low dose), 
P: (mid dose) and P; (high dose). 

He considers both estimation of a for 6 known and estimation of both a 
and £. In all cases the mean square error for the minimum normit x? method was 
less than for the maximum likelihood method. 

Detailed computing procedures are presented for the minimum normit x? 
method plus a numerical example. As an aid to computation, the following tables 
have been prepared : 


(1) Normits (v) to 5D for p = 0(.001)1.000. 

(2) Normit weights (w = 2*/pg), used in estimating and to 4D for 
p = 0(.001)1.000, where z is the ordinate of the normal curve for a given 
pandg=1-—p. 

(3) Antinormits (f) to 5D for » = 0(.001)2.499, used to estimate P for a 
v estimated from the regression line. 


R. L. ANDERSON 
North Carolina State College 
Raleigh, N. C. 


78(K ].—Artuur Grav & HERBERT SOLOMON, “Distribution of quadratic forms 
and some applications,” Ann. Math. Stat., v. 26, 1955, p. 464-477. 


Tables I and II of this paper are probabilities P(Q, < #) for k = 2, 3 and 
t = .1(.1)1(.5)2(1)5 to 4S where Q, = S}_, ax? and x; are normally and inde- 
pendently distributed with zero mean and unit variance,->> a; = 1, a; > 0. The 
probabilities for k = 2 are given for a; = .5(.1).9, .95, .99, 1 and for k = 3 for 
(a1, @2) = (1/3, 1/3), (.4, .3), (.4, 4), (.5, .3), (.6,.2), (5, .4), (.6, .3), (.7, .2), 
(.8, .1). Various approximations are investigated for accuracy. Use of the tables 
is illustrated with several applications. 


W. J. Drxon 
University of California 
Los Angeles, California 


79[K ].—G. S. WEtsH, “‘A tabular method of obtaining tetrachoric r with median- 
cut variables,” Psychometrika, v. 20, 1955, p. 83-85. 


A table is given for obtaining tetrachoric r when it is possible to make cuts at 
the medians. The computing chart [1] is used to set up a table which can be used 
if the proportion in the plus-plus cell is known. The table gives rie, to 3D, for 
proportions to 3D. 

J. R. VATNSDAL 


State College of Washington 
Pullman, Wash. 


L. CHesere, M. Sarrer, & L. L. THurstone, Computing Diagrams for the Tetrachoric 
Cubtaation Coefficient, University Chicago Bookstore, Chicago, 1933. 
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80[K].—S. R. BroapBent, “Quantum hypotheses,” Biometrika, v. 42, i955, 
p. 45-57. 


The author considers populations consisting of normally distributed com- 
ponents, particularly if the means of the components differ by a constant amount 
or quantum, 25. The hypothesis that the means of the components are equally 
spaced is called the quantum hypothesis. The problems treated are the estimation 
of the quantum which determines the spacing of the modes, the estimate of the 
scatter within each subdivision, and a test of the quantum hypothesis. Actually 
the author tests the hypothesis of a rectangular distribution in a subdivision, 
against the alternative of the quantum hypothesis. This test makes use of the 
statistic s*/* where s*? is the lumped variance. A table of the critical values of 
s*/# is given for sample sizes of nm = 20(5)100(50)1000 to 4S for the 5%, 1%, and 
.1% probability levels. 


L. A. AROIAN 
Weapon Systems Laboratories 


Hughes Aircraft Co. 
Culver City, California 


81[K ].—GaretH HorsneLt, “The determination of single-sample schemes for 
percentage defectives,’’ Applied Stat., v. 3, 1954, p. 150-158. 


In sampling a lot of unknown quality for attributes where each item may be 
classified as ‘“‘good” or ‘“‘defective,’’ it is possible to design a single sampling in- 
spection plan which will accept practically all (say a = 95%) of the lots which 
are of fraction defective 1, for example, and which will reject practically all lots 
(say 1 — 8 = 90%) which are of fraction defective p2(f2 > p:). The single 
sampling inspection plan for attributes will be given by N, the number of items 
to be drawn at random for the sample, and K the maximum number of allowable 
defective items in the sample for which the lot will be accepted. In the case of 
sampling binomial populations, this problem has been solved by Grubbs [1] for 
the case a = .95 and 6 = .10. 

The present work by Horsnell covers the case of sampling Poisson populations, 
i.e., percentage defective less than about 10, and gives some very useful tables for 
several acceptance and rejection probability levels. In Horsnell’s notation, P(a) 
is the percentage defective for which the chance of acceptance of the lot is a and 
P(8) is the percentage defective for which the chance of acceptance of the lot is 8. 
[P(8) > P(a) and a > £.] Horsnell gives tables for the Poisson case for a = .99 
and a = .95 and 6 = .1, .05, .01 which are tabulated usually to 4S in the form 
K = 1(1)20 vs. NP(.99), NP(.95), P(.10)/P(.99), P(.05)/P(.99), P(.01)/P(.99), 
P(.10)/P(.95), P(.05)/P(.95), P(.01)/P(.95) and NP(.50). 

Thus, one takes the ratio of the objectionable percentage defective P(8) to 
the acceptable percentage defective P(a) for a lot, looks up this ratio in Horsnell’s 
table for the selected a,and 6 and finds the required acceptance number K im- 
mediately. The sample size N is then found by dividing the corresponding value 
of NP(a) to the K found by the acceptable percentage defective P(a). 

The tables of the article can be used to obtain sufficient points to plot the 
operating characteristic of the sampling inspection or the chance that a lot of a 
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given percentage defective will pass the inspection scheme—and also for finding 
confidence intervals for the true percentage defective of the lot for an observed 
number of defectives in a sample of size N drawn at random from the lot. 


F. E. GruBBs 
Ballistic Research Laboratories 
Aberdeen, Maryland 


1. F. E. Grusss, “On designing single sampling inspection plans,” Ann. Math. Stat., v. 20, 
1942, p. 242-246. 


82[K].—S. H. Brooks, “The estimation of an optimum subsampling number,” 
Amer. Stat. Assn., Jn., v. 50, 1955, p. 398-415. 


This paper discusses the use of pilot samples to estimate the optimum sub- 
sample size in two-stage sampling. For given ratios between the cost of sampling 
and population variance in primary and secondary units, and assuming infinite 
primary and secondary populations, the pilot design of minimum cost such that 
the average efficiency of the final design relative to the use of the true optimum 
will be 90%, is given in Table I for cost ratios <1, 2, 4, 8, 16, 32, 64, 100 and 
variance ratios <4, 1(1)4, 6, 8, 12, 16, 24, 32, 48, 64. For cost ratios, 1, 2, 4, 8, 
16, 32, 64, 100 and subsample size, 1(1)10, 12, 16, 25, 50, 100, Table III gives the 
upper and lower limits to 1D of the range of values of the variance ratio for which 
the efficiency will exceed 90%. 

These tables and the associated theory should be valuable in connection with 
other uses of the basic two-factor nested design. 


C. A. BENNETT 
Hanford Atomic Products Operation 
General Electric Company 
Richland, Washington 


83[K].—SoLomon Ku upack, ‘‘The Poisson distribution and its differences up 
to the 6th order,” 1947, vi + 100 p., 26.7 cm. Introduction in typescript ; body 
of table lithoprinted. No copies available but plates are available. Address 
inquiries to Professor Solomon Kullback, George Washington University, 
Washington, D. C. 


These tables give values of the Poisson frequency function, (e~*a*)/x! to 7D 
for a = 0(.001).01(.01).4 and to 6D for a = .5(.1)15(1)100 which are taken from 
Molina [1]. In addition the first six ascending differences are given to the same 
D for the same values of a. They were computed to facilitate the use of the Gram- 
Charlier type B series which is briefly discussed and illustrated in the introduction. 

c. ¢..& 
ouin E. C. Motina, Poisson’s Exponential Binomial Limit, D. van Nostrand Co., New York, 


84[K].—M. A. Grrsuicx, H. Rusin, & R. Sitcreaves, “Estimates of bounded 
relative error in particle counting,” Ann. Math. Stat., v. 26, 1955, p. 276-285. 


This paper gives a method of estimating the parameter A in a Poisson distribu- 
tion by gradually expanding the area (or time) observed until a fixed number M 
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of events occur. If ay denotes the necessary area, and y a constant, then confidence 
limits for \ are b/[(1 + y)au] and b/[(1 — y)ay_] with confidence level 


ers xM—le-2 d 
a= ae 
wary (M — 1)! 


A preassigned bound on the percent error in estimating A is given by 1007. For 
M(1i — ¥’ 1 
fixed M and y, the value b = b* = ( y Y) log ; + maximizes a, and, for 
pag 

fixed -, M may be chosen as the least integer so that by using 5*, a is just larger 
than a preassigned level of confidence. 

Table I, using 6 = b*, gives values of a to 4D for M = 2(2)40, y = .01, .05, 
.10, .20. Most of these tabled a turn out to be small. The largest, for M = 40, 
are .0503, .2479, .4733, .7988 corresponding to the four values of y. An approxi- 





7 
mation formula is vau (2 + Zr) = Z, where 2, is the 100(1 + a)/2 per- 


centile of the unit normal distribution. A listing fora = .90, .95, .99 and y = .05, 
.10 shows this approximation as good for these cases as a more complicated one 
which in turn checks with the tabled values for M = 40. 


FRANK MASSEY 
University of California 


Los Angeles, California 


85[K ].—Lro Katz, “Probability of indecomposability of a random mapping 
function,” Ann. Math. Stat., v. 26, 1955, p. 512-517. 


Tables are given to 5D for the expression, 


(N — 1)!¥=** mi 
(N — 1)! z, oo 


for 1 = 0, 1; N = 2(1)20(2)40(5)100. A formula is given for large N. These 
quantities are the probabilities referred to in the title for 7 = 0 corresponding to 


a general case and i = 1 to the “hollow” case defined in the article. They are 
computed from values of P given in Molina’s tables [1] by the formula, 


(N — 1)!e%P(N; N —i — 1)/(N — 1)’. 
W. J. Drxon 
University of California 
Los Angeles, California 


1. E. C. Moura, Poisson's Exponential Binomial Limit, D. van Nostrand Co., New York, 1945. 
86[K ].—P. G. Moors, “The properties of the mean square successive difference 


in samples from various populations,” Amer. Stat. Assn., Jn., v. 50, 1955, 
p. 434-456. 


Three estimators of the parent population variance (¢*), based on samples of n, 
are compared: the usual sample variance (s*) and two definitions based on suc- 
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cessive differences, 


n—1 


= > (xi41 — x4)?/(m — 1), 


7 = > (x25 - X2-1)?/2m, 
i=1 


where 7? is only used for n = 2m. 

The first four moments of each estimator are derived. The relative efficiencies 
as estimators of o? (for fixed mean) are tabulated to 2D for = 5(5)25, © and 
various values of the kurtosis parameter, B2(= y/o‘). 

The bias of the estimators is studied when the population mean shifts; only 
s? is seriously biased. Also the variance of s? increases materially, the amount 
being inversely proportional to #2. 

Since & has good properties for shifting populations, its distribution is studied. 
This distribution is best approximated by a Type VI Pearson curve but is close 
enough to a Type III (x?) for most purposes. Four types of parent populations 
are considered to represent different values of 6: and the skewness parameter 
B1(= us/o*): Normal (6: = 0, 82 = 3), Rectangular (6; = 0, 8. = 1.8), Double 
Exponential (B1 = 0, B2 = 6), and Type Ill (Bi = 1,82 = 6;6:1 = 2, Bs = 4). 
Exact and Type III and Type VI approximate distributions of 5*/¢? are considered 
for 2 = 10(10)50, 75, 100. The following tabulations are given: 


(1) Values of 8; and £2 to 4D 
(2) Prob [(6/o*) < £] for selected — to 5D for normal parent only 
(3) Upper 5% points for 5*/c? to 2D for normal parent only. 


Values of 8: and 8; of s? and log s* are also tabulated fer normal (to 2D) and 
Type III parent populations (to 4D). The log transformation is also advo- 
cated for 6. 

Four examples are given to illustrate the usefulness of 5? as an estimator of o°. 


R. L. ANDERSON 
North Carolina State College 
Raleigh, N. C. 


87[K ].—E. S. Pace, ‘Control chart with warning lines,”’ Biometrika, v. 42, 1955, 
p. 243-257. 


First we consider four of Page’s five rules. y is the “‘ideal’”’ process mean, 
o the (known) process standard deviation, N the (regular-interval) sample size, 
# the sample mean. 

Rule I. Choose k, n, N. Take samples of size N. Take action if any point (2) 
falls outside the action lines (u + 610/ VN) or if any & of the last points falls 
outside the warning lines (u + B20/ VN, B: < 81). Now consider a score x; assigned 
to the ith sample as follows: x; = —a (if the sample point is within the warning 
lines), = (if the sample point is between the warning and action lines), =c (if 
the sample point is outside the action lines) (¢ > 0,c > 6 > 0). This scoring 
scheme leads naturally to Rule II. 
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Rule II. Take action after the mth sample if any of the inequalities, 
DX xXn-1 > hyy1, 7 = 0,1, 2, --- s, are satisfied, s and h, suitable constants. 
i=0 

Limiting s only by the number of samples taken since action was last taken 
but requiring 4, = h for all r, the following two restrictions on Rule I are readily 
obtained (see text and references) : (1) = 2 for any n, or (2) k = n. We are led to 

Rule III. Choose n, N. Take samples of N. Take action if 2 points in any 
sequence of n fall between the warning and action lines or if any point falls outside 
the action lines. 

Rule IV. Choose n, N. Take samples of N. Take action if m consecutive points 
fall between the warning and action lines or if any point falls outside the ac- 
tion lines. 

For the case of a constant percentage of inspected output, Page suggests that 
the average run length (the average number of units inspected before the scheme 
requires action when the parameters remain constant) is a good criterion for 
selecting an inspection scheme. Page derives and tables average run length to 0D 
applied to the mean of a normal population, known variance; for Rule IV as 
follows (the mean yp’ = uw + Ag): 


(1) N=5, n = 3,4, By = 3.00(.125)3.25, Bz = 1.00(.25)2.00, A = .00(.02)1.8 
(2) N = 10, n = 3,4, Bi = 2.875(.125)3.25, Bz = 1.00(.25)2.00, \ = .00(.02)1.0 
(3) NW = 15, m = 3,4, By = 2.75(.125)3.125, Bz = 1.00(.25)2.00, A = .00(.02)1.0 
(4) N = 20, m = 3,4, Bi = 2.75(.125)3.125, Bz = 1.00(.25)2.00, A = .00(.02)1.0. 


Finally Rule V is: Choose and N and take action if (i) any point falls outside 
the action lines, or (ii) ” consecutive points fall outside the warning lines, or (iii) 
two out of any set of m consecutive points fall outside opposite warning lines. 
For Rule V tables show the average run lengths to 0D for B; = 2.875(.125)3.25 
and B; = 1.5(.25)2 for the two cases: wu’ = » + Ao, o’ = o with A = 0(.2)1 and 
uw’ = p, o’ = ko with k = 1(.25)2.25. 


HAROLD FREEMAN 
Massachusetts Institute of Technology 


Cambridge, Massachusetts 


88[K].—Fe.rx Ronce, “Die Verhdltnisschatznung (ratio estimate) nach der 
Methode des ‘Veranderungsfaktors’ und der ‘additiven Veranderungsgrosse,’ ”’ 
Mitteilungsblatt fiir Math. Stat., v. 6, 1954, p. 221-232. 


In this paper the author introduces an additive adjustment method for esti- 
mating the total of some population from a sample when information is available 
from some previous period. Thus Y,* = XN(g — Z), where Y;* is the estimate 
for the second period, X is the total for the population at the first period, N is 
the number in the population, 7 is the mean of the sample for the second period, 
and Z is the mean of the same sample for the second period. 

A comparison of the efficiency of this estimate with the standard ratio esti- 
mates is made. This relative efficiency is found to depend upon the correlation, r, 
between observations made on the two different occasions, the ratio, s, of the 
standard deviations of individual observations, and the ratio, f, of the two means. 
For small s and small f or for large s and large f, the additive adjustment method 
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is found to be more efficient than the ratio estimate. Tables of the relative effi- 
ciency to 4D for r = .6(.1).9, .95, S = .4(.2)2, and f = .6(.2)2.4 are given. 


C. F. Kossack 
Purdue University 
Lafayette, Ind. 


89[K ].—J. E. Jackson & E. L. Ross, “Extended tables for use with the ‘G’ test 
for means,’’ Amer. Stat. Assn., Jn., v. 50, 1955, p. 416-433. 


Define Gi, G: as follows: G: = |X — u|/R where X is the mean of a sample of 
size mn, and R is the average range of m random non-overlapping subgroups of n; 

= |X, — X:|/R where the samples from which X;, X: are obtained are made 
up of m, mz random non-overlapping subgroups respectively, of size m each and 
R is the average range of the m; + mz subgroups. 

Table I contains percentage points of the distribution of G; for significance 
levels a = .10, .05 and .01; = 2(1)15(2D); m = 1(1)15(2D). Table II contains 
percentage points of the distribution of G: for percentage points a = .10, .05 and 
.01; 2 = 2(1)15(2D) ; mi and mz = 1(1)15(2D). 

Both tables were computed directly from Lord’s tables of the percentage 
points of u, [1] using indicated transformations. For certain values of m, not 
included in the original tables, corresponding percentage points for u’s were 
obtained using the interpolation formula suggested by Lord. 

The indicated tests are of the two-tail variety. 

P. N. SOMERVILLE 
Statistical Engineering Laboratory 
National Bureau of Standards 
Washington 25, D. C. 


1. E. Lorn, “The use of range in place of standard deviation in the ¢ test," Biometrika, v. 34, 
1947, p. 41-67. - 


90[K ].—H. F. Drneman, “A computing chart for the point biserial correlation 
coefficient,’’ Psychometrika, v. 19, 1954, p. 257-259. 


A computing chart is given for quick estimation of a point biserial correlation 
coefficient when a normally distributed continuous variable is artificially di- 
chotomized at the median. Use is made of [1] and [2] to make a chart which is 
claimed to have an accuracy to 1D. 


J. R. VATNSDAL 
State College of Washington 
a Wash. 


CuestrE, M. Sarrir, & L. L. THURSTONE, Computing Diagrams for the Tetrachoric 
Comeianion Cosficient, University of rs Bookstore, Chicago, 1933. 
2. W.B EL, N. C. Perry, & J. P. Gurtrorp, “The estimation of a point biserial coeffi- 
cient of correlation from a phi coefficient,” Brit. Jn. Psych., Stat. Sec. 1952, v. 5, p. 139-150. 


91[K].—R. F. Tate, ‘‘The theory of correlation between two continuous vari- 
ables when one is dichotomized,” Biometrika, v. 42, 1955, p. 205-216. 


Karl Pearson’s biserial correlation coefficient r* is shown to be consistent and 
asymptotically normal as an estimate of p, the bivariate normal correlation. If w 


is the (standardized) point of dichotomy, and p(w) = f ‘ (24)-te-*/*dt, the dis- 
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tribution of r* depends upon p and ~. Soper [1] obtained the asymptotic variance, 


AV(*) = A+ o| et 4 eons mae s/2| +7 


r 


which is a function of p and p only, A being the normal density with argument w. 
Table 2 gives {nA V(r*)}!, the square root of the expression in curly brackets 
above, to 3D (3 and 4S) for + p = 0(.10)1 and p or (1 — p) = .05(.05).5. 


Leo Katz 
Michigan State University 


East Lansing, Michigan 


1. H. E. Soper, “On the probable error for the biserial expression for the correlation coeffi- 
cient,”” Biometrika, v. 10, 1913, p. 384-390. 


92[K].—J. L. Hopces, Jr., “A bivariate sign test,”” Ann. Math. Stat., v. 26, 
1955, p. 523-527. 


Consider n pairs of bivariate vectors (x;, y;) and (x,’, y,’), i = 1, ---,#. The 
2n vectors of this set are assumed to be statistically independent. The null hy- 
pothesis is that the bivariate distribution for (x;, y;) is identical with that of 
(x;’, ys’) for « = 1,---,m. The alternative hypothesis of interest is that the 
probability distribution for (x,’, y;’) has been shifted relative to that for (x;, y,) 
and that this unknown shift direction is roughly the same for all z. For given shift 
direction, the problem reduces to the sign test by projecting the vectors of 
differences (x; — x,’, yi — y,/) onto the given direction and counting the number S 
of these vectors which have a specified one of the two possible direction senses. 
Consider all possible shift directions and let M be the maximum of the correspond- 
ing S values. The null hypothesis is rejected if M is too large. A procedure is given 
for determining M for all possible cases. The quantity K = m — M is used for 
the derivations and tabulation. For k < 2/3, 


Pr (K <k) = in — 24) (") Jo, 


The table in this article lists values of Pr (K < k) to 5D for m = 1(1)30 and 
k < n/3. 


J. E. WatsH 
Military Operations Research Division 
Lockheed Aircraft Corporation 
Burbank, California 


93[K].—D. R. Cox & A. Stuart, ‘“‘Some quick sign tests for trend in location 
and dispersion,’ Biometrika, v. 42, 1955, p. 80-95. 


In a linear regression model with y; = a + Az + «;, (¢ = 1,2, ---, N), the «; 
are independent standardized normal variates. The null hypothesis, A = 0, is to 
be tested against alternatives A > 0, using a (quick) test based on comparisons 
of independent pairs of the observations. It is shown that the best weighted sign 


test is given by S; = pm LA (N — 2K + 1)hx, w—-x+1, and the best unweighted sign 
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2N 
test is given by S; = D4, he, 7 + K where, for i < j, hij = 1Lify; > yjand 


hi; = 0 otherwise. Since the asymptotic efficiency of S; relative to S; is 96%, 
the authors recommend use of S;. It is also established that the asymptotic rela- 
tive efficiency of S; compared to the best parametric test based on 6, the sample 
regression coefficient, is 83%. 

Table 3 gives the exact power of S; to 3D for sample sizes N = 15(15)135, 
significance level the largest value <.05, and alternatives given by p = .50, .49, 


—V2NA ' 
.45(.05).05, where p = ® Pe ee the unit normal c.d.f. The values of » are 


identified as the true significance level a, and the rows a, b, ---, 7, give powers 
for the various A’s used. Table 4 gives the exact power of the 5 test to the same 
accuracy for the same alternatives. Table 3 also shows, in parentheses, the values 
of the corresponding standardized regression coefficient, A, to 4D. The smallest 
ratio of powers tabulated is 50.3%, for N = 30, p = .25 and A = .0477. 


Leo Katz 
Michigan State University 
East Lansing, Michigan 


04[K ].—Evetyn Fix & J. L. Hopcgs, Jr., “Significance probabilities of the 
Wilcoxon test,” Ann. Math. Stat., v. 26, 1955, p. 301-312. 


The problem considered is that of obtaining significance probabilities for 
the Wilcoxon unpaired two-sample test [1]. Let m <n be positive integers 
while R; < R. < --- < R, represent a random sample of size m drawn with- 
out replacement from the first m + m positive integers. Let S; = R; — i and 
U = S,+ --- + S, while r(u, m,n) = Pr (U < x). The value of x(u, m, n) can 


be expressed in the form A (u, m, n) / ( oe > ") , where 


A(u,m,n) = > (—1)*Ai[u — kn — 4k(kR +1), m — k] 


and the summation is over all non-negative values of the integer k for which 
u — kn — 4k(k + 1) is non-negative. Table I contains exact values of Ao(u, m) 
for m < 12 and u < 100 while Table II contains exact values of A2(u, m) for 
m < 11 and u < 75. In general, A, (u, m) for any k can be found from 


A,(u,m) = Y Ax_i(u — rk, m) 


where summation is over all non-negative integers r such that u — rk is non- 
negative. To extend the range of validity past m < 12 and u < 100, the Edge- 
worth approximation to (u, m, m) is given to terms of order 1/m?. This approxi- 
mation appears to be accurate to about 4D for m > 12. 


J. E. WaAtsH 
Military Operations Research Division 
Lockheed Aircraft Corporation 
Burbank, California 


1. F. Wiicoxon, “Individual comparisons by ranking methods,’ Biometrica, v. 1, 1945, p. 
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95[K].—Des Raj, “Relative efficiency of gauging and exact measurement in 


estimating the proportion of a population between given limits,” Sankhya, 
v. 15, 1955, p. 191-196. 


Suppose that we want to estimate the proportion, P, of a population falling 
in the interval (a, 6). The author compares the large sample relative efficiency of 
two methods for estimating this proportion ; one method consisting of measuring 
each item in the sample, the other method consisting of using go-not go gauges 
set at the values a and b. Tables for the proportion P and the efficiency are given 
for the Type III, Cauchy, and normal distributions, for the case where a and b 
are symmetrically located about the mean of the distribution, and for the one- 
sided case where b = «. In the normal case efficiencies are tabulated to 3D for 
those values a = uw — do for which A = 0(.05)3.00. 


BENJAMIN EPSTEIN 
Stanford University 


Stanford, California 


96(K].—R. E. GrrEENwoop, “Coupon collector’s test for random digits,” 
MTAC, v. 9, 1955, p. 1-S. 


In an ordered set of digits beginning with a specified position one may deter- 
mine the length of the sequence of digits required in order to have the complete 
set 0,1, ---,9 appear in it. The author proposes that the comparison of the 
observed lengths of non-overlapping sequences with the theoretical distribution 
of such lengths, on the assumption that the probability that any one of the ten 
values be assumed by any digit is 1/10, be used as a test of randomness in 
the ordered set. If », is the probability that a sequence have the length n, 

—1 
?. = — m = 10,11, ---. In order to tabulate values of ,, the author ex- 
tended the table of differences of zero of Fisher and Yates [1 ]. His table is claimed 
to give exact values of p, for m = 10(1)35 and values to 20D for.m = 36(1)75. 
A check through m = 26 disclosed one error, apparently typographical. For 
n = 13, p, should read 0.0080 8315 2 instead of 0.0080 9315 2. This correction 
was noted with an obvious misprint, MTAC, v. 9, p. 224. This coupon collector’s 
test is illustrated on the 2035D approximation to # and the 2500S approximation 
to e. A x? test, using grouped frequencies, gave no indication of non-randomness 
in these two interesting ordered sets. 
. 4.2. 


1. R.A. FisHer & F. Yates, Statistical Tables for Biological, Agricultural and Medical Research, 
3rd edition, London, 1948, Table XXII. 


97[K ].—G. W. Tompson, “Bounds for the ratio of range to standard deviation,” 
Biometrika, v. 42, 1955, p. 268-269. 


This table gives upper and lower bounds (distribution free) for the ratio of 
the range w to standard deviation estimate s, both from the same sample of size n. 
It can be shown that the upper bound of w/s is V2(m—1) and the lower bound 
of w/s is 2Vm — 1/m for n even and 2Vn/n + 1 for n odd. 

A table is given for m = 3(1)20(10)100, 150, 200, 500, 1000, to 3D. For samples 
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of size 3 from a normal population, a table of upper and lower percentage points 
of w/s is given for percentages 0, 0.5, 1.0, 2.5, 5.0, and 10.0, to 5D. 
These can be used for routine checks of computation of s. 


J. R. VATNSDAL 
State College of Washington 
Pullman, Washington 


o8[K ].—T. J. Terpstra, “A non-parametric test for the problem of k samples,” 
K. Ned. Akad. v. Wetensch., Proc., v. 57, s. A., 1954, p. 505-512. 


Let X1, ---,X be independent random variables and let Hy be the hy- 
pothesis that all the X; have the same continuous distribution function. Let 
xi,n(h = 1, ---, ) bea set of 2; observations on X;(i = 1, ---, k). For i ¥ j, let 
Ui; + 3nm; equal the number of times that x;,, > xj,m(h < n;,m < nj), while 
U;; = 0. Finally, 

2. 
o=-6y FHL 12 rt 


ig mm n+i1 Fn; 








(X Ui;), 


where 2 = >> n;. Let n; = aym + 0(Vm), where a; and m are positive integers, 


a; being independent of m. It is shown that, under Hp and for large m, the sta- 
tistic Q has approximately a x?-distribution with k(& — 1)/2d.f. Further, in Table 
1 the exact distribution of Q to 4D is given for k = 3, mi < mz < m3 < 3, while 
Table 2 gives the corresponding upper and lower .10, .05, and .01 significance 
levels also to 4D. In Table 3, for »: = mz. = n; = 3, the exact distribution of Q 
is compared with the limiting x;’-distribution, indicating that for n; > 3 the 
x?-approximation is quite satisfactory. 

J. H. B. KemPpeRMAN 
Purdue University 
Lafayette, Ind. 


99[K ].—M. Zia up-Din, “Expression of the K-statistics ky and kio in terms of 
power sums and sample moments,” Ann. Math. Stat., v. 25, 1954, p. 800-803. 


To the k-statistics of Fisher expressed in terms of the power sums of the 
observations in a sample of m which had been given through weight 8 by Dressel 
[1], the author adds ky and kyo. 

Ce GC. 


1. Paut L. DressEL, ‘Statistical Semi-invariants and their estimates with particular emphasis 
on their relation to algebraic invariants,” Ann. Math. Stat., v. 11, 1940, p. 33-57. 


100[K, W].—D. R. Cox, ‘A table for predicting the production from a group 
under the care of one operative,” Roy. Stat. Soc., Jn., ser. B, v. 16, 1954, 
p. 285-287. 


For the case of N machines under the care of a single operator, assuming 
random occurrence of stoppages and an exponential distribution of clearing times, 
the expected fraction of the operator’s time utilized is tabulated to 3D for 
x = .005(.001).01(.01).12 and Nx = .1(.1)1.9 (also Nx = .85(0.1)1.15 for 
x < 0.04), where x is the product of the mean clearing time and mean ratio of 
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-occurrence of stops. This table essentially extends Tables 3 and 4 of Benson and 
Cox [1] to larger NV. A simple relation enables one to compute the expected pro- 
duction rate per machine. 


C. A. BENNETT 
Hanford Atomic Products Operation 


General Electric Company 
Richland, Washington 


1. F. Benson & D. R. Cox, “The productivity of machine requiring attention at random 
intervals,” Roy. Stat. Soc., Jn., ser B, v. 13, 1951, p. 65-82. 


101[L ].—NATIONAL PuysicaL LaBoratory, Tables of Weber Parabolic Cylinder 
Functions. Computed by Scientific Computing Service Limited, Mathematical 
Introduction by J. C. P. Miller, Editor. H. M. Stationery Office, London, 
1955, 233 p., 28 cm. Available from British Information Services, 30 Rocke- 
feller Plaza, New York 20, New York. Price $11.34. 


This book offers much more than the title seems to promise. More than one 
third of the book is taken up by a mathematical introduction containing all the 
information needed, by a computer, on Weber’s differential equation and its 
various solutions. The remaining part of the book contains the numerical tables. 

Weber’s differential equation is written in the alternative forms 


eon 2 a 
(1) qa +a)y =0 

pa = one 
(2) a3 + (3x? — a)y = 0. 


In terms of the parabolic cylinder function D,(z), and of the function 


4ie*I'(3 + ia) 
I'(} — ia) 


relevant solutions of (1), for real a, may be defined as 


E(a, x) = I D_ia+(xe~**"") 


U(a, x) = D_.+(x), U (a, x) = tan xraU(a, x) + sec raU(a, —x) 
Via, x) ad U(a, x)/T(3 — a), V(a, x) = U(a, x)/T(3 =? a), 


and relevant solutions of (2), for real a, as 


W(a, x) = Bt Re E(a,x), W(a, —x) = k*Im EG, x), 
where 
k = Vi + & — ev, 


These solutions have been selected as being most suitable for numerical tabulation. 

The mathematical introduction by J. C. P. Miller contains power series ex- 
pansions; asymptotic expansions (ji) for large |x|, (ii) for large |a|, (iii) for |a| 
and |x| both large; integral representations; relations with Bessel functions and 
other special functions; a very extensive collection of formulas; a description of 
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the tables together with an account of their preparation and instructions for 
interpolation; and a bibliography of previous tables and researches connected 
with parabolic cylinder functions. Both differential equations, (1) and (2), are 
included in this discussion. In the course of the work certain slowly varying 
auxiliary functions are introduced for computational purposes. 

The numerical part contains tables relating to solutions of (2). Certain 
numerical tables of solutions of (1), especially for integer values of 2a, are avail- 
able but on the whole a systematical tabulation of solutions of (1) has not as yet 
been undertaken. 

Almost all computations were carried out by Scientific Computing Service, 
Ltd., and were supervised by J. C. P. Miller who was supported in this task by 
the late L. J. Comrie, and by C. W. Jones. The Mathematics Division of the 
National Physical Laboratory cooperated in the production of the volume. 

Table I gives W(a, x), W(a, —x) and reduced derivatives for a = —10(1)10, 
x = 0(.1)10. W is generally given to 8S, and as many reduced derivatives are 
given as can be used effectively in interpolation in the x direction. The a-interval 
is too large to permit interpolation in the a direction. 

Table II gives auxiliary functions in the non-oscillatory region. log W(a, +x) 
and the first derivatives of these functions are given to 8S with second and fourth 
(or modified second) central differences for a = 1(1)10, x = 0(.1)2Va. 

Table III gives auxiliary functions in the oscillatory region: these are defined 
by the equations 
E(a, x) = Fe*x, dE(a,x)/dx = — Ge, 


F, G, x, ¥ real. 8S values of F, G, x, ¥, with differences as in Table II, are given 
for a = —10(1)2, x = 0(.1)10 and a = 3(1)10, x = 2Va(.1)10. 

Tables IV to VI give auxiliary constants: e+**, k+#, 1 — k+! for a = 0(.1)10; 
e-** for a = 0(.001).1; W(a, 0) and W’(a, 0) for a = 0(.02)1(.1)5(.2)20; real and 
imaginary parts of log '(m/4 + ia/2) for m = 0(1)3, a = 0(.1)5(.2)20. 

This book is a model of presentation. It contains all the material needed by 
the computer, and its arrangement is admirable. The printing is excellent. 


A. ErpgELyYI 
California Institute of Technology 


Pasadena, California 


102[L ].—M. ScouLer & H. GEBELEIN, (a) Acht- und neunstellige Tabellen zu den 
elliptischen Funktionen (Eight and Nine Place Tables of Elliptical Functions), 
Springer, Berlin, 1955, xxiv + 296 p., 26.7 cm. Price DM 58, (b) Fiinfstellige 
Tabellen zu den elliptischen Funktionen (Five Place Tables of Elliptical Func- 
tions), same publishers, 1955, xi + 114 p., 23.8 cm. Price DM 29.60. 


Both these volumes contain tables of Jacobian elliptic functions and theta 
functions, or related auxiliary functions, on a new system. The novelty lies chiefly 
in the arguments. As modular argument, instead of the usual modular angle (here 
called ©) or of k? = sin? ©, use is made of Jacobi’s nome g (or of g* or g*). As the 
other argument, instead of the usual argument u of the Jacobian elliptic functions 
or of the argument x = ru/2K often used for the theta functions, the authors 
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make use of z = cos 2x = cos (ru/K). The aim of the new system is ease of 
bivariate interpolation. This aim has been notably achieved, and the tables 
promise to be useful when interpolation in both arguments is unavoidable; such 
interpolation is sometimes avoidable when © or &? is a natural argument in the 
problem to which elliptic functions are being applied. 

The authors consider that the advantages of their system would be lost if 
ordinary bivariate interpolation were used. Their bivariate tables are meant to 
be interpolated first in one variable and then in the other, either way round 
according to circumstances. Consequently these tables are given twice over, first 
arranged with g (or g* or g*) as top argument and z as left argument, with differ- 
ences in 2, and then vice versa. 

In the detailed descriptions which follow, the first-mentioned argument in 
bivariate tables is the top argument, and differences relate to the second-men- 
tioned (left) argument. All logarithms are to the base 10. The functions G, H, 
G, H used by the authors are related to ordinary theta functions by the equations 


#1(x) = 2gtsinxG #8,(x) = A 
G=1-¢6 A =1+@H. 
The detailed contents of (a) are: 


Table I. G(q*, z) to 9D with A? for g* = .001(.001).1, z = —1(.05) + 1. Also 
q to 7D and ® to 0’.0001, both without differences, for the same g‘*. 

Table II. G(q*, 2) to 9D with A for z = —1(.05) + 1, g* = 0(.001).1. Also 8 
and g again. 

Table III. H(g'*, z) to 9D with A? for g* = .002(.002).176, = —1(.05) + 1. 
Also g to 7D and @ to 0’.0001, both without differences for the same gq’. 

Table IV. H(q', z) to 9D with A for z = —1(.05) + 1, g* = 0(.002).176. Also 
© and 2 again. 

Table V. log (snu/sinx), log (cn u/cosx), logdmu to 8D with A for 
q = .01(.01).55, z = —1(.05) + 1. Also —log cos 8 to 8D, © to 0’.01, and K, 
K/E to 8D, all without differences, for the same gq. 

Table VI. log (sn u/sin x), log (cn u/cosx), logdnu to 8D with A for 
z= —1(.05) +1, g = 0(.01).55. Also —log cos 8 to 8D and 9 to 0’.01, both 
without differences, and K, K/E to 8D with A, all for the same g. 

Table VII. 1/(1 — q), K, K/E, all to 8D with A for —log k’ = —log cos 8 
= 0(.005)3. The values of K and E for —log k’ > 0.5 were obtained from a table 
by E. L. Kaplan [1]. 


The detailed contents of (b) are: 


Table I. log (sn u/sin x), log (cn u/cos x), log dn u to 5D with A or A? (the sec- 
ond difference being modified when italicized) for g = .01(.01).5,z = —1(.1) +1. 
Also @ to 0’.01 without differences for the same gq. 

Table IT. log (sn u/sin x), log (cn u/cos x), log dn u to 5D with A or A? for 

= —1(.1) +1, g = 0(.01).5. Also © to 0’.01 and —log cos 0, K, K/E to 5D, 
all without differences for the same g. 

Table III. G(q, 2), H (q, 2) to 5D with A or A? for g = .91(.01).5, 2 = —1(.1) +1. 

Also © to 0’.01 without differences for the same gq. 
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Table IV. G(q, 2), H(g, 2) to SD with A? for s = —1(.1) + 1, ¢ = 0(.01).5. 

Table V.1/(1 — ¢g), K, K/E to 5D with A, and @ to 0’.01 without differences, 
all for —log k’ = —log cos @ = 0(.01)2.5. 

Table VI. Coefficients of second differences for Everett interpolation of both 
function and first derivative, both to 4D with A for fraction of interval = 0(.01)1. 


Both volumes have explanatory text in German and English with worked 
examples and illustrative diagrams. 


A. F. 


1. E. L. Karwan, “Auxiliary table of commune elliptic integrals,’’ J. Math. Phys., v. 25, 1946, 
p. 26-36, RMT 311, MTAC v. 2, 1946-7, p. 127-128. 


103[L].—F. G. Tricom1, “Valori numerici di funzioni ortogonali di Laguerre,”’ 
Atti. R. Acc. Sc. Torino, v. 90, 1955-1956, p. 1-8. 


This is a revision and extension of the author’s earlier table [1] which was 
reproduced in Jahnke-Emde [2 ]. It gives 1, (x) = e-*L,(x) = e-* > es — 
wo 1S =m! 


to 6S for m = 0(1)10, x = 0.1(.1)1(.25)6(1)14(2)34. The range was chosen so as 
to include all zeros of all the functions. The present table was prepared on a desk 
calculator. The computations were carried out to 8S using the recurrence relation, 
(nm + 1)la41 = (2n + 1 — x)l, — ni,_1. This method is known to be dangerous 
due to the building up of rounding-off errors, so that thorough checking is essential. 
n 2 
Consider £,(x) = n!L,(x) = > (") (n — m)!(—x)™. This has integral 
m=0 





values for integral arguments, and these can be calculated exactly by the use of 
the recurrence relation 


Lai = (2n +1i- x)Ln — n'L_-1. 


From these the corresponding values of /,,(x) were obtained by multiplication by 
e—*/n!. The other values tabulated were obtained by use of an addition formula [3] 


I(x + h) = [In(h) — bai (h) Vox) + [ln—1(h) — Ina) Vi(x) + --- 
+ [hi(A) — lo(A) Vai (x) + lo(h)la(x), 


and a bisection formula [3] 


ln(3x) = e2-" [ + (") h(x) + (*) Ia(x) +--+ + (") h(e)|. 


The table was compared with that of L. J. Slater [4], where the arguments 
coincided. After multiplication by the appropriate exponential there was agree- 
ment, up to 1 or 2 units in the last place, except for 6 values. Dr. Slater informed 
Professor Tricomi that his values were correct; they are 


Le (5) = —2.090 278, Lr (5) = 0.325 397, 
Ls (3.5) = —1.252 437, Ly (3.5) = —0.584 567, 
Ly(3.5) = 0.221 114, L(3) = —0.700 022. 
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Dr. Slater informed the reviewer that errors were possible whenever a different 
routine was used for the calculation when scaling difficulties arose: this routine 
came into action at x > 4.1, 3.4, 2.8, 2.4, 2.1, 1.9 for nm = 5, 6, 7, 8, 9, 10 
respectively. 

5. Ti 


1. F. G. Tricomr, ‘‘Generalizzazione di una formula asintotica sui polinomi di e sue 
someone Accad. delle Scienze di Torino, Cl. d. sci. fis., mat., e nat., Ati, v. 76, 1941, p. 

2. E. JAHNKE & F. Empe, Tafeln hiherer Funktionen, ed. 4, Teubner, Leipzig, - p. 32-33. 

3. F. G. Tricoma, Vorlesung en tiber Orthogonalreihen, Springer, Berlin, 1955, p. 218. 

4. Lucy J. Stater, “A cot table of the Laguerre polynomials,” Inst. ‘lec. E 
Monograph 136, 1955. To be repu blished in Inst. Elec. Engineers, Proc., Part C. [See TAC. 
v. 10, 1956, p. 174-175, which pra sarre further bibliographical material. ] 


104[L ].—NBS Applied Mathematics Series, No. 37, Tables of Functions and of 
Zeros of Functions, U. S. Gov. Printing Office, Washington, D. C., 1954, 
ix + 211 p., 30 cm. Price $2.25. 


This is the first of a series of volumes, to be called Collected Short Tables of the 
National Bureau of Standards Computation Laboratory, and to contain some tables 
previously printed elsewhere as well as certain tables, hitherto unpublished, which 
have been found useful in the operation of the Computation Laboratory. 

In the present volume 18 tables are given. Of these, 15 have either been pub- 
lished or reviewed in MTAC; the remaining three have not been previously 
published. All tables are reproduced by a photo-offset process, and in the case of 
published tables the explanatory text is also reproduced. In addition, there is an 
Introduction by A. N. Lowan to all 18 tables, giving, in some cases, additional 
references. 

1. On the function H(m, a, x). See MTAC 1, 1944, RMT 135, p. 156. 


2. Table of the integrals f * Jo(t)dt and f * Yo(t)dt. See MTAC 1, 1944, 
0 0 


RMT 133, p. 154. 

Table of Jio(x) and related functions. See MTAC 1, 1944, RMT 134, p. 155. 
Table of f,(x). See MTAC 1, 1945, RMT 181, p. 363-364. 

The function E,(x). See MTAC 2, 1947, RMT 392, p. 272. 

Table of the Struve functions. See MTAC 2, 1947, RMT 387, p. 268-269. 
Table of Fourier coefficients. See MTAC 1, 1944, RMT 157, p. 192-193. 
Table of sines and cosines. Gives sin x and cos x for x = 100(1)1000 to 8D 
(guaranteed to within a unit in the last place). This supplements the tables in 
NBS AMS 43. 

9. Radix tables for finding logarithms to 25 decimal places, by Herbert E. 
Salzer. Gives 28D values of log,x for x = 2(1)10, and log, (i + 10) for 
m = 1(1)14, together with instructions for using this brief table for the computa- 
tion of logarithms. 

10. Table of x*/n! for x = .01(.01)1.99 to 13D for m taken up to the point 
where x*/n! vanishes to 13D, and also for x = 1(1)10, m = 1(1)40 to 8S. (Both 
tables are guaranteed to within one unit in the last place.) 

11. Zeros of Legendre polynomials. See MTAC 1, 1943, RMT 92, p. 51-52. 

12. Zeros and weight factors of Laguerre polynomials. See MTAC 4, 1950, 
RMT 737, p. 86. 


PNAVPoe 
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13. A short table of series . . . . See MTAC 2, 1946, RMT 263, p. 37-38. 
14. Zeros of certain Bessel functions. See MTAC 1, 1945, p. 353-354. 

15. More zeros of certain Bessel functions. See MTAC 2, 1946, p. 118-119. 
16. Zeros of the derivative of Bessel functions. See MTAC 7, 1953, p. 69-71. 
17. Complex zeros. See MTAC 3, 1949, p. 351-352. 

18. Roots of sin z = z. See MTAC 1, 1944, Review 50, p. 141. 


A. Erpéty1 
California Institute of Technology 


Pasadena, California 


105[L, M, S].—A. AsHour & A. Sasri, “Tabulation of the function y(@) 

= ine — nm .” MTAC, v. 10, 1956, p. 57-65. 

This article contains a table of the function }F., sin n6@/n?, @ = 0°(10’)180°, 
6D, with first differences. It is prefaced by an interesting introduction which 
describes the functional properties and the method of computation. This series is 
of particular interest from the point of view of numerical analysis. The reader 
might contemplate the difficulties he must face in evaluating this series to the 
accuracy given in the table by straightforward summation on the most modern 
high-speed computer. The authors avoided these difficulties by working with the 
integral representation 





v(0) = y(0, 6) = - f log (2 sin 32)dt 


to first cbtain the values from @ = $a to by numerical quadrature. Thus the 
logarithmic singularity at ¢ = 0 was avoided by making use of the relation 
¥(0, $¢ + 0) = — ¥(3a + 0, x). In the range ¢ = $a to ¢ = x the integrand is 
continuous and, working with an interval of 10’, the trapezoidal rule, corrected 
to take account of the first derivative, gave satisfactory results. The values of 
the integrand logio sin $¢ and its derivative were available from existing tables. 
The values from @ = 0 to @ = 4” were obtained with the functional equation 
¥(0) = ¥(r — 6) + 4y(26). The function is tabulated for 6 = 0°(10’)180° to six 
places, since they were all readily obtainable as a result of the method of compu- 
tation. Although the table is not linearly interpolable near the beginning of the 
table, first differences are given throughout. Since 


v(4r — 6) = Glog, 2 — > (—)"" 3 (228 — ee sade > 
=, (2m)! ™ 2n(2n +1)’ 





we have ¥(0) ~ (x — 6) log, 2 near @ = =x and so the tables could have been at 
intervals of 20’ rather than 10’ from 150°-180° and still satisfy the criterion for 
linear interpolation. As a check on the tables the sixth differences were computed 
for successive groups of seven entries using the last value of one group as the first 
value of the next. The results indicate that the values are probably correct to 
within a unit in the last place given. The given first differences were checked by 
summation and it was found that the first difference for @ = 179° 40’ should be 
— 2017, not —2016. 
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The tables could have been reduced in size considerably if standard practice 
were followed. For exampie in the neighborhood @ = 0° the authors give the 
representation 

8 6 log|@ : "3B ae 
vv) = —- og | Fore (2n)! *” On(2n +1) 





The presence of the logarithmic singularity makes interpolation unfeasible near 
6 = 0°. However, if the auxiliary function J = ¥(6) + @log|@| were given we 
would obtain (using the tabular entries) 


Y 

0° .000000 
1° 017453 
2° .034908 
3° 052362 
4° .069818 
°° .087276 


The resulting table is linear and the 31 entries given in the table could be replaced 
by the six entries above. The values of ¥(@) could be given to make a total of 12 
entries instead of 31. The values of @ log|@| can be obtained with the aid of a table 
of logarithms. 

For the remainder of the table if a criterion for interpolation with modified 
second differences using Everett’s formula had been adopted the following 
intervals would have been sufficient 


6 Interval No. of arguments 
5°-10° 30’ 10 
10°-30° te 20 
30° io y ad 30 
90°-180° 5° 18 


Thus the total number of arguments necessary to provide the same informa- 
tion would be 84, with 90 functional values and their corresponding modified 
second differences. The resulting table would require less than one page rather 
than the six pages as given. 

It may be argued by some that the ideas mentioned above are given from the 
point of view of the table maker and table users would prefer the table in its 
present form, long as it is. However, the appearance of this table and one or two 
others in earlier issues indicates a change in editorial policy of MTAC. Previously 
only short tables of unusual interest were given in MTAC. It is the opinion of the 
reviewer that this policy should be continued and consequently every device 
leading to economization in the size of tables appearing in MTAC should be 
employed. This will avoid the possibility that there be a slack in articles on 
numerical analysis, bibliographic information, reviews, etc. It is clear that the 
editors are trying to fill a need in providing a place where tables can be published— 
a need which was indicated at the Conference on Mathematical Tables held at 
MIT in September, 1954 (see MTAC 9, 1955, p. 42-44). However, it is question- 
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able whether MTAC should be used for this purpose—it would seem that tables 
might be better published in a supplementary issue or perhaps a new journal be 
started for this purpose. 

This function had been tabulated for much coarser interval by T. Clausen [1 ] 
and Clausen’s table is reprinted in [2]. 


A. ABRAMOWITZ 
National Bureau of Standards 


ste gt Laboratory 


T. Crausen, “Uber die Zerlegung reeller gebrochener Funktionen,” J. f. d. reine u. angew. 
Math x4 8, 1832, p. 298-300. 
. F.W . NEWMAN, The Higher Trigonometry, McMillan and Bowes, Cambridge, 1892. 


106[L ].—J. R. Puiuie, “Numerical solution of equations of the diffusion type 
with diffusivity concentration-dependent,’’ Commonwealth Scientific and 
Industrial Research Organization, Faraday Soc., Trans., No. 391, v. 51, part 7, 
1955, p. 885-892. 


This table contains a table of the function A (x) for x = 0(.2)1(.5)4(1)10, 3D. 
The function A(x) = ierfc (x/erfc x). The calculation of the value of the func- 
tion was based on an alternate exponential of the function and the use of tables 
in [1]. An asymptotic expansion for A (x) is included. 

Then a tabulated function is introduced in connection with an iterative 
method for the solution of a diffusion equation of the type 


s-2(02) 
. i... Bs 


-_ J. Comrie, Chambers’s Six-Figure Mathematical Tables, v. 2, Chambers, London, 1959, 
p. 518. 


where D is a function of 0. 


107[M1.—A. Erpéiy1, W. Macnus, F. OBERHETTINGER, & F. G. TrICoMI, 
Tables of Integral Transforms, McGraw-Hill Book Co., Inc., New York, 
Toronto, London, v. 1, 1954, xx + 391 p. Price $7.50; v. 2, 1954, xvi + 451 p. 
Price $8.00. 


This work is unique in more than one respect. It is the first extensive collection 
of functional relations involving integrals of all those functions which may be 
called hypergeometric in the widest possible sense. There are integrals containing 
Gauss’ series, its familiar special cases (such as Legendre functions, Jacobi poly- 
nomials and the complete elliptic integrals), its limiting cases (such as Bessel 
and Whittaker functions), and many of its generalizations (among them Meijer’s 
G-function and Appell’s and Horn’s hypergeometric functions of two variables). 
There are also integrals involving Theta functions. There are no integrals with 
Mathieu or Spheroidal wave functions, and no systematic attempt is made to 
list integrals involving elementary functions only. (For the latter the authors 
rightly refer to the recent excellent tables by W. Meyer zur Capellen and by 
W. Grébner and N. Hofreiter.) 
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Unique, too, is the organization of the tables. If possible, the integrals are 
listed as integral transforms. This means that they are written in the form 


b 
£0) = [°K ns@ex, 


where the kernel K and the limits a and 6 are standardized. The following trans- 
forms have been considered (in order to adhere to conventional notations, other 
letters than x and y are occasionally used) : 





No. of 
Name of transform K(x, y) a b pages 

Fourier cosine cos (xy) 0 2 56 
Fourier sine sin (xy) 0 oo 54 
Exponential Fourier e~ itv — oc 8 
Laplace ev 0 2 98 
Inverse Laplace e*v c—10 c¢ 1t@ 75 
Mellin xv 0 oo 33 
Inverse Mellin y* c—is« c¢ 1 26 
Hankel J, (xy) (xy)? 0 oo 88 
Y-transform Y, (xy) (xy)! 0 2 25 
K-transform K, (xy) (xy)* 0 2 29 
H-transform H, (xy) (xy)! 0 © 16 
Kontorovich-Lebedev Kiz(y) 0 20 3 
Riemann-Liouville fractional 1 

integral T(u) O97: ° 7 ” 
Weyl fractional integral oe (x — y)- y 2 12 

T(x) 

> 1 
Stieltjes G+) 0 oo 23 
Hilbert ; — 2 © 20 

ah 


For each transform the functions f(x) are arranged according to a certain hier- 
archy of functions which is explained in the preface to volume I. These integral 
transforms fill the whole of volume I and about three-fifths of volume II. The 
remaining part of volume II is taken up by miscellaneous definite integrals which 
cannot be conveniently written as an integral transform of one of the above types. 

Only for a few of the better-known transforms (the first five in the above list) 
the authors have been able to avail themselves of significant previous collections 
of integrals. Most of the other material is either new or has been extracted from 
the periodical literature. Professor Oberhettinger contributed most of the miscel- 
laneous integrals in volume II. 

The following remarks concerning accuracy and completeness of the tables 
should be regarded in the light of the unquestionable merit and great competence 
of the work under review. Unfortunately there exists for tables of formulas no 
simple device for detecting errors such as differencing. As a consequence, few 
such tables can match the high standards of accuracy which are now maintained 
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by the best numerical tables, and also the table under review contains a number 
of misprints and minor errors. The reviewer has been informed that an incomplete 
list of corrigenda to volume I has been prepared and is available from the pub- 
lisher. A more extensive list, covering both volumes, is in preparation and will 
be available at the same place. Turning to completeness, it cannot of course 
reasonably be expected that all integrals which have ever been evaluated can be 
found in the tables. One would expect, however, that for those integrals which 
are listed the results are given in complete form. There are some exceptions to this, 
particularly in the tables of Hankel transforms, where some of the given results 
do not cover the entire range of y. 

The editors of MTAC have agreed to publish desirable future addenda to the 
tables. To this the reviewer would like to add a suggestion of his own. Very 
generally speaking, integral transforms are frequently met in problems with 
infinite domains. Problems with finite domains often are solved conveniently by 
series expansions. This suggests the usefulness of tables of series expansions 
involving higher transcendental functions. These tables could be arranged in 
much the same way as the present Tables of Integral Transforms and would form, 
in this reviewer’s opinion, a welcome companion volume to the present set. 


P. HENRICI 
University of California 
Los Angeles, California 


108[S, X].—BERNARD FRIEDMAN, Principles and Techniques of Applied Mathe- 
matics, John Wiley & Sons, Inc., New York, 1956, ix + 315 p. Price $8.00. 


Professor Friedman has written a most excellent and useful book which will 
go far in bridging the gap between the pure and applied mathematics of the 
spectral theory of ordinary and partial differential equations. In any book of an 
introductory nature of this size it is impossible of course to cover the range of 
applied mathematics in anything like the sense it is covered in Frank-von Mises, 
and so Friedman has wisely concentrated on providing the mathematical tools 
appropriate for scattering problems in acoustics, electromagnetic theory and in 
quantum mechanics. In the course of doing this he has at the same time covered 
a great deal of common ground essential for all of applied mathematics. 

The book consists of five chapters entitled: I. Linear Spaces; II. Spectral 
Theory of Operators; III. Green’s Functions; IV. Eigenvalue Problems of Ordi- 
nary Differential Equations; V. Partial Differential Equations. The pedagogy is 
excellent throughout. Friedman proceeds from the simplest example of a given 
phenomenon to a sufficiently general example of its occurrence to provide complete 
understanding of it. The reader is further helped by the fact that the more diffi- 
cult mathematical theorems are proved in a series of appendices at the end of the 
appropriate chapter. Thus the main thread of any given argument is more readily 
followed since a better separation of principles from detaiis is achieved. 

Novel features not usually found in previous books on these topics include: 
the use of the bra and ket notation of Dirac; the introduction and systematic use 
and justification of the distributions of L. Schwartz; the extensive use of contour 
integration in the manner of the Watson transformation; the use of the concept 
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of the “approximate spectrum” for the determination of the continuous spectrum 
of a given differential operator; the treatment of separation of variables by using 
the inverse of two commutative operators ; the method of Kato for improving the 
results of a Rayleigh-Ritz eigenvalue estimation ; the introduction of the concept 
of impedance after Schelkunoff for the usual Sturm-Liouville problem and its 
relations to the concepts of transmission and reflection coefficients ; the systematic 
treatment of Sturm-Liouville boundary value problems when the boundary con- 
ditions are either inhomogeneous or when the eigenvalue enters explicitly in them 
and many others too numerous to be enumerated here. 

While the reviewer would have preferred a more extensive discussion of varia- 
tional methods and of the meaning of the notions of limit point and limit circle 
for differential equations and their relationship to the usual abstract Hilbert 
space formulation of the spectral theory, he believes that Professor Friedman has 
written a thoroughly enjoyable and useful book which all active applied mathe- 
maticians should certainly have as part of their library, and which would certainly 
represent at the same time a suitable text-book for either a semester or a year 
course in applied mathematics. (By way of incidental criticism, the concepts of 
covariant and contravariant vectors are never defined although on page 22 their 
distinction is likened to that existing between basis and reciprocal basis vectors 
and on page 26 to that between row and column vectors.) 


C. L. Dorn 
University of Michigan 
Ann Arbor, Michigan 


109[S].—H. Preuss, Integraliafeln zur Quantenchemie, Springer-Verlag, West 

Berlin, Germany, 1956, iv + 162 p., 27 cm. Price DM 39. 

In the calculation of the binding energy and electronic structure of molecules 
and solids by quantum mechanics, it is necessary to calculate the integral of 
products of wave functions, each about one or another atomic center, multiplied 
by various direct or inverse powers of the distance 7, from one of the centers to 
an element of integration or the distance 7:2. between two elements of integration. 
For first-order approximations only one- and two-center integrals are needed, 
the effect on a bond between two atoms because of the presence of a third is 
neglected. 

The tables under review give numerical values for the simpler pertinent one- 
and two-center integrals, for atomic wave functions appropriate for K and L 
shells (i.e., for constituent atoms up to neon in the periodic table) for a fairly 
wide range of the parameters involved. In addition, a number of subsidiary func- 
tions are tabulated, in terms of which the other pertinent functions may be com- 
puted. Thus approximate calculations of the structure and energy of such mole- 
cules as HzO, NHs, etc., may be worked out by the use of these tables plus a 
desk computer. 

In the functions themselves, five significant figures are given, in the main; 
for the auxiliary functions six significant figures are given. This is adequate, 
since first-order calculations of binding energy often differ from measured values 
by factors of 2 or 3. In order to give a wide enough range of the parameters and 
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keep the tables reasonable in size, the intervals between tabulated values of the 
parameters are large enough so that interpolation would be difficult in many cases. 

These limitations are not critical, however, if the limitations of the first-order 
computations, for which the tables are designed, are kept in mind. All that can be 
expected of first-order calculations is an order-of-magnitude estimate of energy 
and a crude estimate as to structure. It would therefore be illusory to interpolate 
parameter values, to try to get an accurate minimum, for example. To get ac- 
curacy sufficient to warrant interpolation, one would need to use at least second- 
order terms, which involve three-center integrals and other functions not tabu- 
lated here ; such computations should probably be done on a high-speed computer. 

For what they are supposed to do, the tables seem adequate. They are hand- 
somely printed, with plenty of explanatory textual material and bibliographic 
references. The reviewer has not checked the accuracy of the entries. 


Puitrp M. Morse 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 


110[S ].—BatTH. VAN DER Pot, Editor, Atlas of Ground-Wave Propagation Curves 
for Frequencies between 30 Mc/s and 300 Mc/s, \’'Union Internationale des 
Telecommunications, Geneva, 1955, xxxv + 174 p., 29.5 cm., oblong. Price 
36.50 Swiss francs. 


Curves are presented for the field of a half-wave transmitting dipole as a 
function of distance for a wide range of the physical parameters. These are: 
frequencies of 30, 60, 100, 150, 200, and 300 megacycles per second, heights of 
transmitting antenna above the ground of 10, 20, 50, 100, 200, 500, and 1000 
meters, and heights of receiving antennas above the ground of 0, 2, 5, 10, 20, 50, 
100, 200, 500, and 1000 meters. The transmitting and receiving antennas are 
oriented, simultaneously, horizontally, or vertically. The electrical ground con- 
stants are taken corresponding to sea water (¢ = 4 mhos/m and ¢/e9 = 80) and 
to land (¢ = 10 mhos/n and ¢/e¢9 = 10). To account for “normal” atmospheric 
refraction, the earth’s radius is taken as 4/3 times the actual earth’s radius. 

The r.m.s. field strength in micro-volts/meter is plotted against distance D 
from 0 to 200 kilometers for a specified antenna height, frequency and ground 
constants. The strength of the transmitter is such that the hypothetical field, in 
its equatorial plane, would be 


2.22-10° 
7 eet) oe 


E Dun 


V/m 


if the latter was placed in free space. 

One of the underlying assumptions is: the earth is regarded as a smooth and 
homogeneous sphere. The inhomogeneous character of the atmosphere has not 
been considered. Effects such as “inversion layers’’ and turbulence are thus ex- 
cluded. Furthermore, the roughness and variability of the ground constants are 
not mentioned. It is unfortunate that a few calculations were not carried out for 




















. REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 257 


an earth conductivity of 10~* mhos per meter corresponding to poorly conducting 
ground which is more typical of arid regions of the world. 

As a result of the simplifying assumptions mentioned above, the problem is 
reduced to a scalar boundary value problem in electromagnetic theory. Essentially 
it is desired to obtain a Green’s function G exterior to a spherical region with an 
approximate boundary condition 8G/dr = yG at the surface. The constant y 
depends on the electrical constants of the ground, the polarization and the fre- 
quency. The validity of this boundary condition rests on the condition that the 
square of the refractive index N? of the earth is large compared to one. Some error 
might be introduced in employing values of y based on this assumption when 
| N?| is only of the order of 10 for over-land propagation at 30 Mc/s. 

For purposes of computation, G is represented as a sum of modes involving 
spherical Bessel Functions of large complex order with a large real argument. 
These Bessel functions are approximated by using their representation in terms 
of Hankel functions of order 1/3 or the equivalent modified Airey Integrals. This 
Hankel approximation has been discussed extensively in the early work of van 
der Pol and Bremmer referred to in the text. The statement on page XIII that 
“the introduction of the Hankel approximation is essentially equivalent to the 
use of an approximate boundary condition . . .’’ is, however, not true. These are 
two separate approximations with different regions of validity. 

When the transmitting and receiving antennas are at moderate or very large 
heights, geometrical optics provide a much simpler formula for computation. 
Most of the curves in the atlas are based both on the geometrical-optical and the 
spherical wave function formulas. The overlapping of the regions of usefulness 
of these two techniques provides an excellent check on the over-all accuracy. 

Most of the calculations were carried out under the supervision of Professor 
A. Van Wijngaarden on the electronic digital calculator (ARRA) of the Mathe- 
matical Centre in Amsterdam. Parts of the text were written by Dr. H. Bremmer. 

The curves provide a convenient reference for field strength estimates in the 
frequency range from 30 to 300 Mc/s when the atmosphere is undisturbed and 
inversion layers are not present. 

The list of errata below has been compiled since the publication of this book 
and furnished by Professor van der Pol. The corrections apply only to the intro- 
ductory text of the atlas and in no way influence the validity of the curves. 





page for read 
VI, lines —11 and —12 ona perfectly conducting plane in free space 
IX, line 9 1.47 1.047 
IX, line 13 having the same moment radiating 1 kilowatt 
X, line —11 azimuthal angular 
XIV, formula (7) 787 000 
XV, formula (10) 787 Ge Dim 


James R. Wait 


National Bureau of Standards 
Boulder, Colorado 
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111[S].—J. M. Kennepy & M. J. Curr, Transformation Coefficients between LS 
and jj Coupling, Atomic Energy of Canada, Limited, Chalk River, Ontario, 
A.E.C.L. Report No. 224, 1955, 39 p., 27 cm. Price $0.50. 


This report lists and discusses the transformation coefficients for converting 
wave functions expressed in the LS coupling representations to a jj representation 
or vice versa. The basic relation used is: 


hah 
V(LileL, $1525, JM) = Zz A le Se je V(Lisiji, lesoje, JM), 
ide a ee 


where /; and /; are orbital angular momenta whose resultant is L; s; and Sse are 
spins whose resultant is S; ji, the total angular momentum of particle 7, is the 
resultant of J; and s:, and J, the total angular momentum of the two particle 
system, is the resultant of L and S (and also of j; and j2). The expansion coeffi- 
cients A are independent of M the z component of the angular momentum J. 
In atomic and nuclear spectroscopy, the spins s; and sz are equal to }. In the 


report, the coefficients, 
heh 
All, 4 je 
Boda 


are given in three different forms: 


1. as algebraic functions of /;, /2, J, 

2. numerically (up to 8 decimal places) for all possible cases in which /; and 
lz do not exceed 5, 

3. in coded fractional form for all cases in which /; and J. do not exceed 5. 
The coded fractional form is useful since the squares of the A coefficients 
are rational fractions which have no prime factors exceeding the value 19 
for the ranges of parameters considered here. 


We give here an example of this method—it is discussed more fully in the report. 
18/275 = 2’xX 3? X 5? X7° X 117 — 1220, 1, 


the underlining denoting a negative exponent and a comma being placed after 
the fourth digit. Multiplication of fractions merely involves addition of their 
fractional representations. 

The tables are well presented (and legible) and should be useful to workers 
in atomic and nuclear spectroscopy. 

Related tables are: 


W. T. Suarp, J. M. Kennepy, B. J. Sears, & M. G. Hoyze, Tables of Coeffi- 
cients for Angular Distribution Analysis, CRT-556, A.E.C.L. Report No. 97 
[MTAC, Rev. 26, v. 10, 1956, p. 51-52]. 

J. M. Kennepy, B. J. Sears, & W. T. SHarp, Tables of X Coefficients, CRT- 
569, A.E.C.L. Report No. 106 [MTAC, Rev. 27, v. 10, 1956, p. 52]. 

S. A. MoszKowskI 


University of California 
Los Angeles, California 
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112[C, K, L, W, Z].—Rosert R. Bus & Freperick MosTeLier, Stochastic 
Models for Learning, John Wiley & Sons, Inc., New York, 1955. Price $9.00. 


This work is of interest to readers of Mathematical Tables and Other Aids to 
Computation both for the text and for tables which are included as incidental 
information to provide means for evaluation of the validity of hypotheses made 
or suggested in the text. 








Major tables include the functions: 4(a, 8) = 9 a’’*"g" and W(a, 8) 
nt log + — ; 
= Peo va? +) /28", a = .5(.02).99, 8 = .5(.02)1, 4D; T(a, 8) = eh" idliew 
on 
the range a, 6 = .7(.01).99 with a < 8, 3D; g,(a) =; a, a = .50(.01).99, 
—_ a’ 
a’B e va’B 





» = 1(1)10, 3D; F(a, 8,2) = r=07 


a, 8 = .5(.05).7(.02).98, @ = 4(4)16. 

In addition to these tables (which were considered by the authors to be 
sufficiently important to be presented separately as appendices) there are numer- 
ous tables scattered throughout the text. Some of these simply tabulate observed 
data from experiments reported, and others furnish functions which the authors 
need for their exposition but which they evidently consider to be of less general 
interest. An example of such a table is the function P(a, 8) = [] 0 (1 — a’8), 
a, 8 = 0(.1)1, 3D found on page 160. 

The major tables are used in maximum likelihood estimates and in other 
adjustment of parameters in the theory developed. 

The book also contains a general description of the mathematical and nu- 
merical methods used in making the estimates. These methods are not new, but 
the whole exposition seems to the reviewer to be a valuable contribution to litera- 
ture generally applicable to the rapidly expanding quantitative social sciences. 

The aim of the authors is to present a hypothesis with regard to the learning 
process. This is briefly that a reaction to a stimulus is probabilistic and that the 
governing probability distribution is dependent on the experience of the person 
subjected to the stimulus. A major portion of the book is devoted to an attempt 
to make quantitative hypotheses along these lines which accord with experience. 

The tables are auxiliary to the estimates required in these hypotheses. 

The theory itself is of plausible interest to students of automatic computing 
equipment. This interest might work in both directions—the equipment might 
give better simulations of the learning processes assumed and hence better insight 
into the theory, and the theory might be put to use in the design and application 
of some types of machinery, say, for language translation or for some similar 
application. However, the connection at present must certainly be considered to 
be nebulous. 

The tables are a minor part of the book; they are a little hard to read because 
of the space saving device of listing two functions with a single table of arguments. 
Acknowledgment of help in computation is given, but no indication of checks of 
accuracy or of means used. There are no aids to interpolation. 


and G(a, @,2) = > <0 [~~ 


— a’B — a’B 


Css. 
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113[Z ].—Martin H. Werk, PB 111996, A Survey of Domestic Electronic Digital 
Computing Systems (Reprint of Ballistic Research Laboratory Report No. 
971), United States Department of Commerce, Office of Technical Services, 
Washington, D. C., 1956, vii + 272 p., 28 cm. Price $4.75. 


It is certainly not easy to remain abreast of developments in the electronic 
computer field. This survey helps to bring up-to-date (end of 1955) a broad picture 
of the field in the United States. Previous surveys of a similar nature date back 
to 1953 and are now obsolescent. 

Eighty-four ‘Domestic Electronic Digital Computing Systems”’ are described 
in the report. Two foreign manufacturers, Ferranti and Olivetti, are represented, 
presumably because these companies have U. S. sales organizations. On the other 
hand at least some U. S. computers are missing, for example Harvard Mark IV, 
Unrvac 120, and Untvac File Computer. It would be highly desirable if all foreign 
computers could be added to the list together with annual supplements. 

The descriptive information on each system runs from 2 to 4 pages and 
generally includes photographs, manufacturer, user(s), arithmetic and logical 
organization, size, cooling and power, input-output and memory data, production 
record and price, and operating experience. 

Much of the information was obtained from manufacturers and users. For 
this reason, as the editor points out, care must be taken when comparing detailed 
minute points since the contributers may have used different frames of reference. 
Anyone considering the acquisition of a computer could use the information in 
the survey to select a small group of computers for further consideration. The 
survey would be very useful for familiarization purposes. 

A chapter on Analysis and Trends follows the system description. The 
machines are listed in a number of ordered tables according to the characteristics 
of word length, add time, memory capacity, memory access time, tube and diode 
quantities, power requirements and cost. In addition the editor gives his opinion 
of trends in the field. 

A brief bibliography and a glossary of computer terms are included. 


Dovucias L. HoGAN 
The George Washington University : 
Washington 6, D. C. 


TABLE ERRATA 


Reviews and papers in this issue mention errata in the following works: 


sin 20 
n? 


A. AsHour & A. Sasri, “Tabulation of the function ¥(@) = \?., 





Review 105, p. 250. 

RoBERT E. GREENWOOD, “‘Coupon collector’s test for random digits,” Review 
96, p. 243. 

H. Hasse, Vorlesungen tiber Zahlentheorie, p. 197. 

F. G. Tricomi, “Valori numerici di ortogonali di Laguerre,’’ Review 103, p. 248. 

BALTH. VAN DER PoL, Atlas of Ground-Wave Propagation Curves for Frequencies 
between 30 Mc/s and 300 Mc/s, Review 110, p. 256. 
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249.—W. MULLER, “Viscous flow within cylindrical boundaries,” Ann. d. Physik. 
[MTAC, v. 1, 1944, p. 263.] 
Items 7 and 9, for Jo(rj1,.) read J1(rj1, 2). 
H. D. Conway 
Cornell University 
Ithaca, New York 


250.—HERBERT E. Sauzer, “On numbers expressible as the sum of four tetra- 
hedral numbers,” London Math. Soc., Ju., v. 20, 1945, p. 3-4. 


The following erratum has been found: 
Page 3, line 25, the number 107 was omitted. It should be inserted between numbers 
103 and 137. 


HERBERT E. SALZER 
CONVAIR 


San Diego, California 


251.—GEORGE WELLINGTON SPENCELEY, RHEBA MURRAY SPENCELEY, & EUGENE 
RHODES EpPEerson, Smithsonian Logarithmic Tables to Base e and Base 10, 
The Smithsonian Institution, Washington, D. C., 1952. 


The following errata have been found. 


p. 336, log 6686, first figure, for 6 read 8 
p. 342, log 6997, third figure, for 9 read 4 
p. 349, log 7347, fourth figure, for 0 read 1 
p. 364, log 8100, first figure, for 8 read 9. 
ARNE OLANDER 


University of Stockholm 
Stockholm 6, Sweden 


NOTES 


Sponsorship by the Office of Naval Research, the National Bureau 
of Standards, the National Science Foundation, 
and various other organizations 


Much of the editorial and refereeing work for this and for earlier volumes of 
Mathematical Tables and Other Aids to Computation has been done under the 
sponsorship of the Office of Naval Research and the National Bureau of Standards 
as part of their program for the development of numerical analysis and research 
connected with numerical analysis. Considerable help has been given to the 
Editorial Committee by various committees sponsored by the National Science 
Foundation and several learned societies in connection with mathematical and 
statistical tables. Finally, the members of the Editorial Committee who are 
members of University faculties and the referees of papers submitted to the 
journal have carried out their work for Mathematical Tables and Other Aids 
to Computation as part of their University research assignment. The contribu- 
tions of all these organizations to the publication of this journal are gratefully 
acknowledged. 





CORRIGENDA 





Index 


The classified index of reviews in this volume continues along the general 
lines set in volume 9 (see MTAC, v. 9, 1955, p. 226-229). However, in accord 
with suggestions received from many different people, an attempt has been made 
to describe the work reviewed rather than to list the title. 

For tables in analysis, the description usually describes the functions tabulated, 
the increments and ranges of the arguments, the precision of the table. The 
description for other tables and works is less stereotyped, ranging from the fairly 
precise description above to vague statements about intended or presumed use 
of the work. 

In the formal description the symbol 6* is used to indicate that m-th differences 
of some kind—possibly modified with throw-back of a higher difference, for 
example-—are included. 

Tables from Number Theory are arranged more or less in the order of the 
classification scheme of D. H. Lehmer [1]. Tables from Analysis are arranged 
more or less according to the classification scheme of A. Fletcher, J. C. P. Miller, 
and L. Rosenhead [2], usually carried to 1D. Notation of these authors has been 
taken as standard wherever possible. Tables from Statistics are arranged more 
or less in the order of the classification scheme furnished last year by Professor 
H. O. Hartley. It has been convenient to add a few classes to these schemes. 

The Chairman of the Editorial Committee accepts full responsibility for the 
mistakes or bad judgment displayed in the classified and other indices. He 
invites suggestions for future volumes. He also acknowledges gratefully great 
aid rendered by Professor W. J. Dixon, particularly in connection with Tables 
from Statistics. 

Usually in Mathematical Tables and Other Aids to Cemputation ranges and 
increments of variables are expressed in the notation, 


@1(d1)a2(d2)as(ds)a4, 


etc., which means that the first value of the argument is a;, the increment is d; 
between argument values a; and az (inclusive), the increment is d, between da 
and as, the increment is dz; between a; and a4, and (in this example) the last 
value of the argument is a.. 

e ST. 


D. H. LEHMeER, erent to Tables in the Theory of Numbers, National Research Council, 
Washington, D. C., 194 


A. FLETCHER, J. c. P. Mitter, & L, ROSENHEAD, An Index of Mathematical Tables, Scien- 
tie Computing Service Limited, London, 1946. 


CORRIGENDA 


Anprew D. Bootu, Numerical Methods, MTAC, v. 10, 1956, Review 38, 
reference 3, p. 166, for 


on p. 53 we find + etc., ---, «++ etc. 
read 


on p. 53 we find + etc. ---, «++ ete. 
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Harvey Coan, “Stability configurations of electrons on a sphere,”” MTAC, 
v. 10, 1956, p. 117. 

Féppl’s results for n = 8 were misquoted. The configurations corresponding to 
the cube and dodecahedron are not stable. 

Mr. John Leech has been kind enough to point this out and suggest the further, 
rather comprehensive reference: L. L. WuyTE, “Unique arrangements of points 
on a sphere,”’ Amer. Math. Monthly, v. 59, 1952, p. 606-611. 

W. M@ ter, “Viscous flow within cylindrical boundaries,” Ann. d. Physik, 
MTAC, v. 1, 1944, p. 263, items 7 and 9 should be corrected in accordance with 
Table Errata 249 this issue. 

K. C. S. Prttar & K. V. RAMACHANDRAN, “On the distribution of the ratio of 
the ith observation in an ordered sample from a normal population to an inde- 
pendent estimate of the standard deviation,’”’ MTAC, v. 10, 1956, Review 13, 
p. 43, line 8 from bottom, for v. 25, 1955, p. 565-572, read v. 25, 1954, p. 565-572; 


p. 44, line 3, for 
eikat/6 (x a.x') 
1=0 


read 


i 
eikz*/6 (= a(t) : 
i=0 


THe RAND Corp., One Million Random Digits and 100,000 Normal Deviates, 
MTAC, v. 10, 1956, Review 11, p. 39, in the title, 


for One Million Random Digits . . . 

read A Million Random Digits .... 
In connection with the number of pages, 

for xxv + 200 p. 


read xxvii + 400 + i + 200 p. 
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F(0, ¢) and E(0, ); @ = 5°(5)90°, @ = 1°(1)90°; 4D 


Jacobian Elliptic Functions 
log (sn u/sinx), log (cnu/wsx), logdn«u with 6; q = 0(.01).55, 
z = —1(.05)1; 8D; also with 6 or &; q = 0(.01).5, z = —1(.1)1; 5D. 
Comparable tables for @ 


Elliptic Functions @,, 02, 0:, 3, 


Auxiliary functions, G(q‘,2), with 6, q and ®; q* = 0(.001).1, 
z = 1(.05)1; G to 9D, q to 7D, @ to 0’.0001. Also related functions to SD 
over reasonable grid. Auxiliary functions H(q*, z) for 33; similar to both 
sets above, 9D and 5D 


Mathieu Functions. Spheroidal Wave Functions. Lamé Functions 
Characteristic values of Mathieu's equation, Be,-(t) = be-(s) — (2r + 1)/t, 
Bo-(t) = bo-(s) — (2r — 1)/t, and mod &; here s = 1/#; r = 0(1)15, 
s > 10 in variable interpolable increments; 6 to 8D 

Solutions of V7W + k*W = 0 in oblate and prolate spheroidal coordinates 


Hypergeometric Functions, Confluent Hypergeometric Functions 
Laguerre polynomials L,(x) and mod &; n = 0(1)10, x = 0(.1)5; 5 
to 6D 

Weber parabolic cylinder functions W(a, x) and W(a, —x) and reduced 
derivatives for interpolation in x; a = —10(1)10, x = 0(.1)10; 8S. 
Much auxiliary material 

Laguerre functions e~*!2 L,(x); m = 0(1)10, x = .1(.1)1(.25)6, (1)14 
(2)34; 6S 

H(m, a,x) = e~** F(m +1 — ia, 2m + 2 = 2ix); x = 1(1)10, a = 1 
(1)10, m = 0(1)3; 7S. 

Zeros of L»(x); m = 1(1)15; 12D. Corresponding weights for quadrature 
12S 


Miscellaneous Functions 
Clauson's integral C\() = — f° log (2 sin 1/2t)dt = Zn sin nd; with 
6; @ = 0°(10’)180°; 6D 
Functions related to molecular structure for atoms up to neon 
Interpolation Coefficients and Formulas 
L(x), where M*, pM*, M*, 2M‘, 54M‘, pM*, ps*M?, 54M?, M5, pM, M3, 


8M? resp. = x ox L(jE~tdt; x = 0(.1)x0, where L(xo) = 0; exact 


Numerical Differentiation 
Wk? (af /axPay**); «a = 1(1)4, B = 0(1)a; in terms of differences of 
f with increments (h, k) 
Aj,viqs P = 1(1)n, 7, ¢ = O(1)m, m = 2(1)10; exact rational fractions. 
Here f® (xq) = 2 Ai, v:af (xs) 


Curve-Fitting, etc. 


Coefficients of bn, where T,(ax +b) = 2 bmTm, and Tm is a Chebyshev 
polynomial; n = 12, 24 





Page 


106-107 


246-248 


246-248 


108 


172-173 


174-175 


245-246 


248-249 
248-249 


248-249 


250-252 


255-256 


38-39 


171-172 


231 














Review Number 
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104 


104 
107 


77 


TABLES FROM STATISTICS 


Description 
Harmonic Analysis—Sines and Cosines 


sin ( i) and cos (» =) v = 0(1)25, m = 1(1)25; 4D. 
sin (355 ) and cos (on St ) v = 0(1)100, nm = 1(1)50; 4D 


Integral Transforms 
ff" edu and related functions; n = 0(1)20, x = 0(.01)2(.1)10; 
7-10D 
Jf 2 sin nvxdz and {” x*x* cosh xxdx;k = 0(1)10, » = 1(1)100; 10D 
Extensive tables of integral transforms 


TABLES FROM STATISTICS 


Normal Distribution. Area, log Area and Ordinates 
as Functions of the Deviate x 


p(x); x = 0(.001)2.499; 5D 


Normal Distribution. Area, log Area and Ordinates as Functions of x/V2 


74 


13 


77 


20 


77 


58 


Error function H(x); x = 0(.001)2(.01)3.09; 5D 


Normal Distribution. Ratios of Ordinate to Area, etc., 
as Functions of the Deviate x, etc. 


Ratios of ordinate to integral and reciprocals; x = 0(.01)3; 5D except 
ordinate to integral 5S 


Normal Distribution. Miscellaneous Functions with 
Deviate x as Argument 


a;™ such that [e"!*p(x)}* = © a,x; i = 0(1)30; k = 1(1)7; 8S 


Normal Distribution. Deviates as Functions of the 
Area p, Percentage Points, Probits 
x(p); P = 0(.001)1.000; 5D 


Normal Distribution. Tables Associated with Probit Analysis 


Minimum working probits, ta. = (y .. g log r) ,2! = p/Z,w = 23/p: 


Y = 1(.1)5; 4D or 5S 
2*/p(1 — p); p = 0(.001)1.000; 4D. Here z is the ordinate 


Higher Integrals of the Normal Area 


eric (x) = f~ a eric Ede; m = —1(1)11, + = 0(.05)1(.1)2.8; 
6-9D 
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Page 


35 


249-250 


249-250 
252-254 


233-234 


106-107 


232 


43-44 


48 


233-234 


176-177 








TABLES FROM STATISTICS 





Review Number 
(For author 
see p. 267-8) Description Page 
Rankings of Normal Samples or Order Statistics, 
Median and Percentiles 
13 (1) 5th and (2) 95th percentiles of qn = xn/x and (3) 95th percentile of 43-44 


ttn =|xn/s|, Xn = max of m samples from N(0, 1), s independently dis- 
tributed with v degrees of freedom; for (1) m = 1(1)8, » = 3(1)10(2)20, 
24, 30, 40, 60, 120, ©; for (2) m = 1(1)10, » = 1(1)5(5)20, 24, 30, 40, 
60, 120, ©; for (3) m = 1(1)8, » = 5(5)20, 24, 30, 40, 60, 120, ©; 2D 


Normal Distribution. Measures of Location and Dispersion 
other than the Standard Deviation 


16 (.95)?-! and (.99)?-! points for q(p, m2), where q(p, m2) = w/s, w being 45-46 
range of p independent normal variates with equal (m, o) and nzs* is from 
x? with no d.f.; p = 2(1)10(2)20, 50, 100, mz = 1(1)20(2)30, 40, 60, 


100, 

24 Tables relating to R = max;,; |x; — x;|, i, 7 = 1(1)4, x; @ sample from 50 
N(mijo); specialized m; 

88 Relative efficiency of additive adjustment method and standard ratio esti- 239-240 


mates of extrapolating population totals; r = .6(.1).9, .95, S = .4(.2)2, 
f = .6(.2)2.4; 4D. r is correlation, S ratio of standard deviations, f the 
ratio of means 
89 Percentage points of G; and Gz for levels a = .10, .05, .01; m = 2(1)15, 240 
m, m,, mz = 1(1)15; 2D. Samples are of size mn and (m, + m:2)n, 
Gi = |X ~ u|/R and G.= |X, as X.|/R 


Tests of Normality 
99 Fisher's ky and kyo in terms of power sums and sample moments 244 
Normal Bivariate Distribution - 
25 Proportions of success in high and low 27-percentiles, pu and pr; 50-51 


.5000 < p < .9713, r = .05(.05)1.00; 4D. Here p is item difficulty and 
r is item discrimination 
79 Tetrachoric r for median cuts; proportions 3D; 3D 234 
Normal Distribution. Miscellaneous. Quality Control 
87 Average lengths of uninterrupted runs in various inspection schemes;0D 238-239 


Ordinates of the x*-Distribution and Poisson Frequencies 
83 ea? /x! and 6, «++ 88; a = .5(1)15(1)100; 6D, and 236 
a = 0(.001).01(.01).4; 7D 
Factorials and the Gamma Function. Moments of x* and x 
41 I'(n) and log I'(m); m = 1.01(.01)2.50; 5D 106-107 


Auxiliary Tables for Testing with I" Variates 


15 Tables giving confidence limits of lot acceptance tests, exponential distri- 44-45 
bution of failures 


Miscellaneous Functions Based on the Poisson Distributions 


(N — 1)! yy; 
w-9i “=~ 


N¥/M1; i = 0, 1, N = 2(1)20(2)40(5)100 237 














Review Number 
(For author 
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12 


14 


17 


78 


112 


21 


22 


91 


23 


TABLES FROM STATISTICS 


Description 
Binomial Distribution. Auxiliary Tables and Charts 
for Testing Percentages 
Pwx (probability of kill by salvo of N bombs on unit square with normal 
circular distribution {t,,¢n} with (&,) normal circular {0, yo, ca} 
for Pix = .1(.3)1; yo = 0, 1, 2, 4, 7, 11, 16, 22; o4, on = 1, 2, 4, 7, 11, 
16, 22; N = 1, 5, 10, 25, 50, 100, 150, 200; 4D 
Critical values for a signifiance test of differences between small propor- 
portions, binomial variates 
Miscellaneous Functions Based on the Binomial Distribution 


E' (np, n) = z+(*) yo; np =1(1)10, 2 = 25, 50, ~, and 
np = 12(2)20, m = 25, 50, 100, and np = 25(5)45, n = 50, 100; 4-5D 


Tables of F-Distributions and Multivariate Variance Tests 


P(Qz < t), where Oy = D}_, aix#, x1 « N(O, 1) and independently k = 2, 
3, ¢ = .1(.1)1(.5)2(1)5; 4S 


Likelihood and Maximum Likelihood Criteria on 
Parameters in Poisson Parents 
Confidence levels for estimation of Poissen distribution parameter by 
timing M occurrences; M = 2(2)40, allowed error = 1, 5, 10, 20%; 4D 


Likelihood and Maximum Likelihood Criteria on 
Parameters in Miscellaneous Parents 


Various tables for maximum likelihood estimates in stochastic learning 
theory 


Distribution of Serial Correlations and Related Statistics 
Tables relating to use of & = X (xis, — x:)*/m — 1 as estimator of o* 


Rank Correlation and Associated Criteria 

For Kendall's rank correlation S, P,{S > So}; k = 3 and m = 2(1)5, 
k = 4 and m = 2(1)4, k = 5 and m = 2, 3, k = 6 and m = 2, and 
So = 0(2)---; 3-4S 
Rating estimates x;, a testing statistic B, and significance probability of B, 
for rank analysis; 4 treatments, 7 and 8 repetitions and 5 treatments, 1 to 5 
repetitions; x; to 2D, B, to 3D, significance to 4D 
Point biserial correlation coefficient for normal variate dichotomized at 
median; 1D chart 

nAV(r*); +p = 0(.10)1, p = .05(.05).95; 3D. Here AV(r*) is 
asymptotic variance of Pearson's biserial correlation coefficient 


Order Statistics Criteria from any Continuous Parent 


Rejection levels for assumption of symmetric distribution about 0 when n 
observations are negative; a = .01, .05 and n = 5(1)2S. Effective signifi- 
cance levels n = 5(1)30 


Non-parametric Two Sample and Analysis of 
Variance Randomization Tests 


Ao(u, m); m < 12 and u < 100; exact A2(u, m); m < 11 and u < 75; 
exact. Here Ax(u,m) = 2 Ax_i(u — rk, m) 
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43 


46-47 


234 


236-237 


259 


237-238 


48 


49 


240 


240-241 


49 


242 








Review Number 
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92 


93 


97 


98 


76 


100 


il 
52 


11 





TABLES FROM STATISTICS 


Description Page 
Runs 
Pn, probability that a run of n decimal digits will be repaired to include 243 
every digit; mn = 1(1)75; exact or 20D 
Miscellaneous Non-parametric Tests 


(n — 2k) (") Jo for use in a bivariate sign test; n = 1(1)30, 241 

k < n/3;5D 

Power of an unweighted sign test of a linear regression model; N= 15 241-242 
(15)135; 3D. Significance levels 

Upper and lower bounds of w/s; sample size n = 3(1)20(10)100, 150, 243-244 
200, 500, 1000; 3D. Percentage points 0, .5, 1, 2.5, 5, 10% for normal 

samples; n = 3; 5D 

A statistic Q nearly x? for testing equality of k distribution functions; 244 

k = 3, numbers of samples < 3; 4D. Also lower .10, .05, and .01 sig- 

nificance levels 4D 


Formulas for Moments of Moments and Cumulants 
E(un), E(un?), o(tn'); m = 3(1)25;4D. Here un = max {S:,52,°--,Sa}, 233 
Sp = Dj xs, x1 N(0, 1) independent 

Distributions of Means from Non-Normal Parents 


5, 1, .1% levels of s*/6; m = 20(5)100(50)1000; 4S. s* is lumped variance 235 
of sample from population of normal components with means differing 
by 2S 


Miscellaneous Tables Relating to Fitting of Frequency 
Distributions, Moments of Moments and Cumulants, 
Means from Non-Normal Parents, etc. 


Production efficiency, N machines; x = .005(.001).01(.01).12, for 244-245 
x. > .04, Nx = .1(.1)1.9 and for x < .04, Nx = .1(.1).8(.05)1.2(.1)1.9; 
3D 


Tables for Harmonic Analysis 
‘ 360 360 
sin ( a) and cos (» a) v = 0(1)25, m = 1(1)25, 4D 35 
: 360 360 
“s ( 33 — (m3 ) v = 0(1)100, x = 1(1)50, 4D 


Random Numbers 
10° random decimal digits 39-43 
10‘ random decimal digits 172 
Random Deviates from Mechanical Devices— 
Quasi-Random Numbers 
Five methods of generating quasi-randcm numbers 229 


Random Normal Deviates 
10* deviates from N(0, 1) 39-43 
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37 
18 
19 


81 


82 
95 


75 


108 


45 


47 


71 
72 


110 


TABLES FROM PHYSICAL SCIENCES 


Description 
Construction of Block Designs 
Maximal sets without isomorphs of Steiner triple systems of order n < 15 


Optimal Design Aid Tables 
F,(x); x = 0(.1)2(.5)3, & = 1(1)5; 3-5D 
AOQL contour curves; AOQL = .01(.01).06(.02).1, i = 1(1)50(2)100, 
f = .01(.001).05(.002).10(.01).5 
Tables of sample size and maximum defectives for single sample control, 
Poisson population 
Pilot designs of minimum cost giving 90% efficiency in two-stage sampling 
Relative efficiencies of methods for estimating proportion of Type III, 
Cauchy and normal distributions between given limits 
Partitions, Symmetric Functions, Differences of Zero, etc., 
Bernoulli and Euler Numbers. See Analysis 
S(r, 2), m = 1(1)50; exact 
Cross index of fundamental symmetric functions of weights < 12—one 
part and unitary functions 
A*%0"- 


1 
= io ;m = 10(1)75; exact for n < 35, 20D elsewhere 


Pa 


Bibliographic Material Relating to Statistical Tables 
Index of tables 
A guide to tables of the normal probability integral 
A list of unpublished tables 


TABLES FROM PHYSICAL SCIENCES 


General References 
Principles and techniques of applied mathematics 


Particle Mechanics 
Tables for rocket and comet orbits 


Fluid Mechanics 
Compressible fluids 


Exchanges, Reaction Kinetics, Reaction Rates, etc. 

Diffusion flames 
Computers, etc., in chemistry 

Electricity and Magnetism; Electric Waves; Electric Filters 
Weights for an n-signal Chebyshev filter 
Coefficients of bm, where T,(ax + b) = Z buTn and T,, is a Chebyshev 
polynomial; n = 12, 24 
Ground wave propagation; frequencies 30-300 Mels; curves 


Crystal Structure, X-ray Diffraction, Electron Microscopy 
Tables for conversion of X-ray Diffraction Angles to Interplanar spacing 
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Page 


104-105 


47 
47 


235-236 


236 
243 


37-38 
230 


243 


100-102 
231-232 
232-233 


254-255 


110-111 


112 


111 
177 


230-231 
231 


256-257 


177-178 








280 AUTOMATIC COMPUTATION 


Review Number 
(For author 
see p. 267-8) Description 
Atomic and Molecular Physics and Chemistry, Quantum Mechanics 
109 Functions related to molecular structure for atoms up to neon 
111 Transformation coefficients, LS and jj couplings 
Nuclear Physics and Scattering 
(Including Photon Scattering and Compton Effect) 
26 Tables of coefficients for angular distribution analysis 
27 Table of X coefficients 
61 Scattering of y-rays by K electrons 


WORKS CONCERNING SOCIAL SCIENCES 





Page 


255-256 
258 


51-52 
52 
178 


Note. Tables relating to experimental design, quality control, and statistical inference 


are listed under STATISTICS, not here 
Operations Research, Simulation 


68 Artificial simulation of traffic 
Econometrics 

15 Truncated life tests 

39 Conference on linear programming 

50 Symposium in linear programming 

63 Programming the transportation Problem 


Management Science 
(For Programming see Econometrics) 


100 A table for predicting the production from a group under the care of one 
operative 
Psychometrics 
25 Item analysis table for the high-low-27-per-cent method 
112 Stochastic learning theory 


Competition, Attrition, Warfare, etc. 
12 Salvo kill probabilities 


AUTOMATIC COMPUTATION 


General Works on Automatic Computation 


59 Computers, etc. in chemistry 
62 What is science? 


Radix Conversion—Decimal to Radix 2 and Radix 2* to Decimal 


6 Conversion of decimal x to periodic octonary fractions; 
x = .0001(.0001).9999 
65 nm X 8*; m = 0(1)511, & = 3(3)12; 10D and  X 10-*; n = 0(1)999, 
k = 3(3)9; 11 octal places 
Logic of Coding 


33 Binary multiplication procedure 





229 


44-45 


169-171 
179 


244-245 


50-51 
259 


43 


177 
178-179 


37 


227 


100 
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Review Number 
(For author 
see p. 267-8) Description Page 
Machine Descriptions—Specific Electronic Machines 
113 A survey of electronic computers made in U. S. 260 
Machine Programs and Codes for Specific Problems 
63 Programming the transportation Problem 179 
64 Stability configurations of electrons on a sphere 179 
Machine Design—Generalities 
112 Stochastic models for learning 259 
TABLE ERRATA 
No. Author Title Page 
247 GrusEPPE PaLama & “Tavola dei numeri primi dell’intervallo 12 012 000- 
L. PoLett1 oo 8 Sok ahh) are er 54 
248 Harorp T. Davis Tables of the Higher Mathematical Functions........... 180 
249 W. MULLER “Viscous flow within cylindrical boundaries”.......... 261 
250 HERBERT E. SALZER “On numbers expressible as the sum of four tetrahedral 
Es foe nak to nk awh atte oa ecaecveekaeussuaavaees 261 
251 GEoRGE WELLINGTON Smithsonian Logarithmic Tables to Base e and Base 10... 261 
SPENCELEY, RHEBA 
Murray SPENCELEY, & 
EUGENE RHODES EPPERSON 
H. Hasse Vorlesungen iiber Zahlentheorie, see remark by Emma 
A? caches Cite tabt ath eeeaaseeousdekeae 197 
Also see the following reviews : 
11 Tse RAND Corp. A Million Random Digits and 100,000 Normal Deviates 39-43 
15 BENJAMIN EpsTEIN “Truncated life tests in the exponential test’’........ 44-45 
22 R.A. BrapLtey & “Rank analysis of incomplete block designs. I. The 
M. E. TERRY method of paired comparisons”................+++++- 49 
24 S#Hozo SHIMADA eS a eR 50 
30 T. LarBLe “Hohenkarte des Fehler-integrals”................... 53 
56 Lucy J. SLATER “A short table of the Laguerre polynomials”. ......... 174 
96 Rosert E. GREENWOOD “Coupon collector’s test for random digits”........... 243 
103 _F. G. Tricom1 “Valori numerici di funzioni ortogonali di Laguerre” 248-249 
in 0 
105 A. AsHour & A. SABRI “Tabulation of the function ¥(@) = D2., ~ —" .. 250-252 
110 BALTH. VAN DER PoL Atlas of Ground-Wave Propagation Curves for Frequencies 
between 30 Mc/s and 300 Mc/s................44- 256-257 
NOTES 
Lin nas beebecesbevercceteebeoccevecesensa 112 
Electronic Computer for the Indian Statistical Institute..................60550eeeeeeees 113 
Vychislitel’naia Matematika i Vychislitel’naia Tekhnika....................00.0000: 113-115 
ns os sc ccheetedinneeturreresansscecesensvio’ 180-181 
Sponsorship by the Office of Naval Research, the National Bureau of Standards, the National 
Science Foundation, and various other organizations... ....... 2.2.2... ..ssseeeeeeee 261 
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